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The physics of Anderson transitions between localized and metallic phases in disordered systems is 
reviewed. The term "Anderson transition" is understood in a broad sense, including both metal- 
insulator transitions and quantum-Hall-type transitions between phases with localized states. The 
emphasis is put on recent developments, which include: multifractality of critical wave functions, 
criticality in the power-law random banded matrix model, symmetry classification of disordered 
electronic systems, mechanisms of criticality in quasi-one-dimensional and two-dimensional sys- 
tems and survey of corresponding critical theories, network models, and random Dirac Hamilto- 
nians. Analytical approaches are complemented by advanced numerical simulations. 
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developments and, in particular, on novel types of criti- 
cal systems. In this review we understand the term "An- 
derson transition" in a broad sense, including not only 
metal-insulator transitions but also critical points sep- 
arating phases with localized states (most prominently, 
quantum-Hall-type transitions). 

We now list the key developments in the field that took 
place during the last decade and constitute the main sub- 
ject of the review. 



A. Symmetry classification and universality classes 

Within the early classification scheme, three universal- 
ity classes for the Anderson transition were identified - 
orthogonal, unitary, and symplectic - in correspondence 
with the Wigner-Dyson classification of random matrix 
theory ensembles. Two basic symmetries of this scheme 
are the invariance of the Hamiltonian under time rever- 
sal and spin rotations. More recent research has shown, 
however, that this picture is in fact by far incomplete, 
for two reasons: (i) there exist more symmetry classes of 
disordered systems, and (ii) in many cases, the symmetry 
class does not uniquely determine the universality class 
of the transition. 



I. INTRODUCTION 



1. Additional symmetries. 



It is known since the seminal work o flAnd ersonI (I1958D 
that disorder can localize a quantum particle inspite of 
quantum tunneling processes and even if the particle is 
not localized classically. For a given energy and dis- 
order strength the quantum states are either all local- 
ized or all delocalized. This implies the existence of 
transitions between localized and metallic phases in dis- 
ordered electronic systems, known as Anderson transi- 
tions. A great progress in understanding of the cor- 
responding physics was achieved in the seventies and 
the eighties, due to the developments of scaling theory 
and field-theoretical approaches to localization, which 
demonstrated connections between the Anderson tran- 
sition and conventional second-order phase transitions. 
These results were summarized i n several review articles : 
(iKramer and MacKinnonl . ll993l : rLee and Ramakrishnanl 
I1985D . and, in the conte xt of the quantum Hall transi- 
tions , ( Huckesteinl . 1 1 995t) . as well as in the book I Efetovl . 
I1997D . 

During the last ten years considerable progress in the 
field has been made in several research directions. This 
has strongly advanced the understanding of the Ander- 
son localization phenomenon and the associated quan- 
tum phase transition physics in disordered electronic sys- 
tems and allows us to view it nowadays in a considerably 
broader and more general context. These important ad- 
vances have motivated us in writing the present review. 
While the article will also include a brief overview of basic 
earlier results, the main emphasis will be put on recent 



It has been understood that a complete set of ran- 
dom matrix theories includes, in addition to the three 
Wigner-Dyson classes, three chiral ensembles and four 
Bogoliubov-de Gennes ensembles. The additional en- 
sembles are characterized by one of the additional sym- 
metries - the chiral or the particle-hole one. The field 
theories (tr-models) associated with these new symmetry 
clas ses have in fact been consider ed already in the eight- 
ies (|Hikamil . 11983': 'Wegneil. Il989l ) . However, it was only 
after their physical significance had been better under- 
stood that the new symmetry classes were studied sys- 
tematically. For the chiral ensembles, important con- 
tributions in this di rectio n were made in Gadc (1S393) 

jSlevin and Nagaol (Il995 



, Gade and Wegneil (11991,. ^..^ 
Verbaarschot and Zahedl ( 19931 ). The particle-hole sym- 



metric ensembles were called i nto life several years later 
( Altland and Zirnbaueil Il997f ). Zirnbauer has also es- 
tablished a relation between random matrix theories, a- 
models and Cartan's classification of symmetric spaces 
fSec. lIVp. which provides the mathematical basis for the 
statement of completeness of the new random ma trix 
classification ( Heinzner et all 120051 : IZirnbaue"^ Il996l ). 

Amongst other things, these developments have led 
the theorists to predict two novel qua ntum Hall effects , 
the spin quantum H all effect (SQHE) (jKagalovskv et al] , 
119991: ISenthil et all Il998h . Sec. IVLDl and the ther- 
mal quantum Hall effect ( TQHE ) ([Chalker et all 120021: 
ISenthil et all . Il999[ ). Sec. lyTEl Both should occur in 
materials with paired fermions where the particle-hole 
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symmetry is realized. 



odd. These developments make a survey of disordered 
wires (Sec. |V]) a natural part of this review. 



2. From symmetry classes to universality classes 

The classification of fixed points governing the local- 
ization transitions in disordered metals has turned out 
to be much richer than that of symmetries of random 
matrix ensembles (or field theories) . The first prominent 
ex ample of this wa s in fact given more than 20 years ago 
by iPruiskenI ( 1984 ) who showed that the quantum Hall 
transition is described by a cr-model with an additional, 
topological, term. However, it is only recently that the 
variety of types of criticality - particularly rich in 2D 
systems - was fully appreciated: 

(i) In several symmetry classes, the field theory (a- 
model) allows for inclusion of the topological 0-term (re- 
sponsible for the quantum Hall criticality) or of the Wess- 
Zumino (WZ) term. 

(ii) The phase diagram may depend on the type of 
disorder. The class D represents a prominent example, 
with three different network-model realizations yielding 
vastly different phase diagrams, Sec. IVI.EI 

(iii) In some cases, the field theory may possess a line 
of fixed points, since the coupling constant corresponding 
to the conductivity is truly marginal. This situation is 
in particular realized in the chiral symmetry classes BDI, 
AIII, and CII, Sec. IVLFl 

(iv) In some cases, the symmetry of the cr-model may 
get enhanced under renormalization, so that the ultimate 
fixed-point theory may have a different form. A paradigm 
for this behavior is provided by the S^ sphere cr-model 
with 9 = n topological term (describing a spin-^ antifer- 
romagnet) which fiows into a SU(2) Wess-Zumino-Witten 
(WZW) model. It was conjectured that a similar mech- 
anism may be relevant to some cr-models of localization, 
including the critical theory of the integer quantum Hall 
effect (IQHE). 

(v) It is possible that the same critical theory is shared 
by systems belonging to different symmetry classes. This 
type of "super- universality" has been proposed to occur 
in disordered wires with critical states, Sec. IV. El IV. PI 

(vi) It was recently discovered that Griffiths effects can 
render the conventional RG analysis of a cr-model insuffi- 
cient. In the framework of the RG calculations the result 
can be recovered if infinitely many relevant couplings are 
kept, Sec. lVI.F.2l 



B. Multifractality of wave functions 

It was appreciated by the beginning of nineties that 
wavefunctions at the Anderson transition exhibit strong 
amplitude fluctuations that can be characterized as 
wave function multifractality. The corresponding re- 
sults are summarize d in the review papers (HuckcstcinI, 
Il995t iJanfienl Il994l ) published about a decade ago. In 
more recent years a considerable progress in the un- 
derstanding o f wav e function statistic s in metallic sam- 
ples (lMirhnl.l2000bh and at criticality (lEvers et aZ.l .l200ll: 
lEvers and Mirlinl . bood iMudrv et all Il996l) has been 
achieved. The multifractality implies the presence of in- 
finitely many relevant operators, which is a peculiarity of 
the Anderson transition critical point, and the spectrum 
of multifractal exponents constitutes a crucially impor- 
tant characteristics of the fixed point governing the tran- 
sition. The understanding of general properties of the 
statistics of critical wave functions and their multifractal- 
ity (see Sec. III. Cl) was complemented by a detailed study 
- analytical as well as numerical - for a number of lo- 
calization critical points, such as conventional Anderson 
transition in various dimensionalities, Dirac fermions in a 
random vector potential, IQHE, SQHE, and symplectic- 
class Anderson transition in 2D, as well as the power-law 
random banded matrix model. 

In several situations, the characterization of a criti- 
cal point by its multifractality spectrum has turned out 
to be particularly important. Specifically, a great deal 
of recent research activity has been devoted to confor- 
mal theories governing Anderson critical points in 2D 
systems. Further, in systems of the new symmetry 
classes, peculiar critical point s have been found that cor- 
respo i id to strong disorder ( Carpentier and LeDoussall . 
I2OOII : iMotrunich et all |2002[ ). such that critical wave- 
functions show at the same time some kind of localiza- 
tion. Entrance of a system into such a strong-coupling 
regime manifests itself as a phase transition in the mul- 
tifractality spectrum (the "freezing transition"). 

Finally, the notion of multifractality was very recently 
extended onto a boundary of a critical system, yielding a 
novel independent set of surface critical exponents. The 
importance of this notion has been confirmed by analyt- 
ical and numerical studies of the surface multifractality 
for several models at criticality fSec. III.C.8[) . 



3. Many-channel disordered wires 

Common wisdom has it that all states in one- 
dimensional disordered systems are localized. However, 
for several symmetry classes wires with critical states and 
even such with perfectly transmitting eigenchannels were 
identified recently. The emergence of criticality depends 
crucially on whether the number of channels is even or 



C. Quantitative understanding of critical behavior 

For several types of Anderson transitions, very detailed 
studies using both analytical and numerical tools have 
been performed during the last few years. As a result, 
a fairly comprehensive quantitative understanding of the 
localization critical phenomena has been achieved. The 



following developments played a particularly important 
role in this context: 



1. Power-law random banded matrix model 

An ensemble of power-law random banded matrices 
(PRBM), which can be viewed as a ID system with long- 
rang e hopping, has been ana lytically solved on its critical 
line (jMirlin and Ever^ I200Q) . This allowed, in particu- 
lar, a detailed study of the statistics of wave functions 
(in particular, multifractality) and energy levels at crit- 
icality. The PRBM model serves as a "toy model" for 
the Anderson criticality. This model possesses a truly 
marginal coupling, thus yielding a line of critical points 
and allowing to study the evolution of critical properties 
in the whole range from weak- to strong-coupling fixed 
points. 



2. Network models 

Formulations of quantum dyna mics in terms of net- 
work models, pioneered in Chalker and Coddingtoiil 
(|l988l ) in the IQHE context, have been developed and 
systematically exploited for both analytical studies and 
computer simulations. Such network models have played 
a key role in advancement of understanding of Quan- 
tum Hall critical points, including the conventional IQHE 
and the sys tems of unconventional symmetries - SQHE 
and TQHE dCh alker et aLV2002"; Cho an d Fisheil ll997al: 
iGruzberg et ah . [1999; : .Read and Ludwiei . l200l[ ) . In par- 
ticular, the investigation of the network model of SQHE 
has led to an analytical understanding of the critical be- 
havior of a number of most important physical observ- 
ables. Sec. IVLDl 



3. Progress in numerical simulations 

During the last ten years numerical mathematicians 
have developed highly efficient routines for diagonaliz- 
ing sparse matrices. Combined with the increase in com- 
puter power and an improved understanding of finite size 
effects, this development has recently paved the way for 
highly accurate numerical studies of critical behavior for 
a variety of Anderson critical points. 



4. Field theories: cr-models and Dirac fermions 

The development of the symmetry classification of dis- 
ordered systems has allowed to classify also the cor- 
responding field theories having a form of nonlinear 
cr-models defined on different symmetric spaces. The 
renormalization-group (RG) method was used to ana- 
lyze them at and near two dimensions. A complemen- 
tary approach is based on the analysis of 2D disor- 
dered Dirac fermions subjected to different types of dis- 



order (jLudwig et all Il994 iNersesvan et all . Il995l ). An- 
alytical methods have allowed to identify fixed points 
and determine the critical behavior for some types 
of disorder corresponding to unconventional symmetry 
classes. The interest to the random Dirac fermion mod- 
els has been largely motivated by thei r applications 
to disordered c?-wave superconductors, see lAltland et al] 
(2002) for review. Recent breakthrough in the fabrica- 
tion of monoatomic grap hene sheets and corresp o nding 
transport nieasure ments ( Novoselov et all l2005l . |2004 
IZhang et all l2005l ) has greatly boosted the theoretical 
activity in this field. 



II. ANDERSON TRANSITIONS IN CONVENTIONAL 
SYMMETRY CLASSES 

A. Scaling theory, observables, and critical behavior 

Quantum interference can completely suppress the dif- 
fusion of a particle in random potential, a phen omenon 
known as Anderson localization ( Andersonl . ll958[ ). When 
the energy or the disorder strength is varied, the sys- 
tem can thus undergo a transition from the metal- 
lic phase with delocalized eigenstates to the insulating 
phase, where eigenfunctions are exponentially localized. 



|^2(r)|^exp(-|r-ro|/0, 



(2.1) 



and ^ is the localization length. The character of this 
transition remained, however, unclear for roughly 20 
ye ars, until Wegne r conjectured, developing earlier ideas 
of iThoulessI ( 1974[ ). a close connection between the An- 
derson tr ansition and t he scaling theory of critical phe- 
nomena (|Wegneil . fl97i) . Three years later, Abrahams, 
Anderson, Licciardello, and Ramakrishn an form ulated 
a scaling theory of localization (Abraham s et al\ .\l979). 
which describes the flow of the dimensionless conduc- 
tance g with the system size L, 



dlng/dlnL ^ P{g). 



(2.2) 



This phenomenological theory was put on a solid ba- 
sis after Wegner's discovery of the field-theoretical de- 
scription of the lo calization prob lem in terms of a 
nonlinear cr-model (jWegned . 1 19791) . Sec. ITLB] This 
paved the way for the resummation of singulari- 
ties i n perturbation theory at or near two dimen- 
sions (jCor'kov et all \197^ IVollhardt and Wolfld . ll980D 
and allowed to cast the scaling in the systematic 
form of a field-theoretical RG. A microscopic deriva- 
tion of the cr-mode l wor k ed out in a number of 
pape r s jEfetov et all 19801 : Jiingling and Qppermannl . 
gneii I: 



119801 : ISchaefer and Wegneii ll980D has completed a case 
for it as the field theory of the Anderson localization. 

To analyze the transition, one starts from the Hamil- 
tonian H consisting of the free part Hq and the disorder 
potential U (r): 



H = Ho + C/(r) 



Ho = p2/2m. 



(2.3) 
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The disorder is defined by the correlation function 
{U{r)U{r')); we can assume it to be of the white-noise 
type for definiteness, 



(?7(r)C/(r')) = (27rpr)-M(r-r'). 



(2.4) 



Here, p is the density of states, r the mean free time 
and (. . .) denote the disorder average. It may be shown 
that models with finite-range and/or anisotropic disorder 
correlations are equivalent with respect to the long-time 
and long-distance behavior (hydrodynamics) to the white 
noise model with re normalized parameters ( tensor of dif- 
fusion coefhcients) (|W61fie and Bhattl . flQSl . 

More convenient for numerical simulations is the lat- 
tice version of (|2.3p . (|2.4p known as the Anderson tight- 
binding model. 



(2.5) 



where the sum (ij) goes over nearest neighbor sites and 
the random site energies Ui are chosen from some distri- 
bution V{u); the standard choice is the uniform distribu- 
tion over an interval [~W/2; W/2] ("box distribution"). 

The physical observables whose scaling at the transi- 
tion point is of primary importance is the localization 
length ^ on the insulating side (say, E < Ec) and the DC 
conductivity a on the metallic side [E > Ec), 



crcx {E~ EcY. 



(2.6) 
(2.7) 



The corresponding critical indices ly and s satisfy the scal- 
ing relation s = v{d — 2), first derived in , Wegner (1971). 

On a more technical level, the localization transition 
manifests itself in a change of the behavior of the diffusion 
propagator, 

n(ri,r2;c^) = (Gi+^/2(ri, r2)G^_/2(r2, ri)), (2.8) 

where (...) denotes the disorder averaging and G^, 
are retarded and advanced Green functions. 



G 



l^^{v,v')^{v\{E-H±^r^)-'\v'), 



j;^+0. (2.9) 



In the delocalized regime 11 has the familiar diffusion 
form (in the momentum space). 



Ii{ci,uj) = 2TTp{E)/{Dq^ -iuj), 



(2.10) 



where p is the density of states and D is the diffusion 
constant, related to the conductivity via the Einstein re- 
lation a — e^pD. In the insulating phase, the propagator 
ceases to have the Goldstone form (12.101) and becomes 
massive, 



n(ri,r2;t^) 



2ttp 



It is worth emphasizing that the localization length ^ 
obtained from the averaged correlation function 11 — 
(G^G'^), Eq. (|2.8p . is in general different from the one 
governing the exponential decay of the typical value 
Iltyp = exp(ln G^'^G"^). For example, in quasi-lD systems 
the two lengths differ by a factor of 4. However, this is 
usually not important for the definition of the critical in- 
dex We will return to observables that are related 
to critical fluctuations of wave functions and discuss the 
corresponding family of critical exponents in Sec. III. CI 



B. Field-theoretical description 

1. Effective field theory: Non-linear cr-model 

In the o rigina l derivation of the g - model 



( EfetoveiaZj. 



, - ■ 19801 I Jiingl ng and OppermannL Il980l : 

ISchaefer and Weeneil . Il980l: IWeeneil . Il979^ ■ the replica 
trick was used to perform the disorder averaging. 
Within this approach, n copies of the system are 
considered, with fields (j)a, ct = l,...,n describing 
the particles, and the replica limit n — > is taken 
in the end. The resulting cr-model is defined on the 
ri — > limit of either non-compact or compact sym- 
metric space, depending on whether the fields 0q are 
considered as bosonic or fermionic. As an example, 
for the unitary symmetry class (A), which corresponds 
to a system with broken time-reversal invariance, the 
CT-niodel target manifold is U(n, n) /U(n) x U(n) in the 
first case and U(2n)/U(n) x U(n) in the second case, 
with n 0. A supersymmetric formulation given by 
lEfetovl (|198,'^ combines fermionic and bosonic degrees 
of freedom, with the field $ becoming a supervector. 
The resulting a-model is defined on a supersymmetric 
coset space, e.g. U(l, 1|2)/U(1|1) x U(l|l) for the 
unitary class. This manifold combines compact and 
non-compact features and represents a product of the 
hyperboloid = U(l,l)/U(l) x U(l) and the sphere 
= U(2)/U(l) X U(l) "dressed" by anticommuting 
(Grassmannian) variables. For a detailed presentation 
of the supersymmetry fo rmalism and its app li cations to 



mesos co pic systems, see Efetovl (1 1983 ^ 1997h : Fvodorov 



1997I)- 



1995 ): iFvodorov and Sommeri 
19981): [Mir lin (2000a, b); Verbaarschot tt al 



Guhr et aL. 
(1985!); 



IZirnbaueilT2004[ ). While be ing equivalent to the replica 



version on the level of the perturbation theory (including 
its RG resummation) , the supersymmetry formalism al- 
lows also for a non-perturbative treatment of the theory, 
which is particularly important for the analysis of the 
energy level and eigenfunction statistics, properties of 
quasi-lD systems, topological effects, etc. 

We briefiy sketch the key steps in the conventional 



H\ri-r2\/0, 



(2.11) 



with the function J^{r) decaying exponentially on the 
scale of the localization length, !F{r/£) ^ exp(— r/^). 



^ A remarkable exception is the behavior of ^ in disordered wires 
of the chiral symmetry, see Eqs. US. 21( 1. US. 22( 1. 



6 



derivation of the cr-model; to be specific, we consider 
the unitary symmetry class. One begins by expressing 
the product of the retarded and advanced Green func- 
tions in terms of the integral over a supervector field 

= {Sl,Xl,S2,X2)- 

G's+„/2(ri, i-2)G^_^/2(r2, ri) 
= y"l?$ i?$t^i(ri)5i*(r2)52(r2)52*(ri) 



X exp 



|i J dr$t(r)[(£;-iJ)A+|+i77]$(r)|,(2.12) 



where A = diag{l, 1,— 1,— 1}. After disorder averaging, 
the resulting quartic term is decoupled via the Hubbard- 
Stratonovich transformation, by introducing a 4 x 4 su- 
permatrix variable 7?.p^(r) conjugate to the tensor prod- 
uct $^(r)$J,(r). Integrating out the <i> fields, one gets 
the action in terms of the TZ fields, 



S[TZ] = npr J d^rStiTZ^ + Stilii[E+{- + iri)A~ Ha-U], 

(2.13) 

where Str denotes the supertrace. The next step is to 
use the saddle-point approximation, which leads to the 
following equation for TZ: 

7^(r) = (2^pr)"i(r|(i; -Ho- Uy^lr) (2.14) 

The relevant set of the solutions (the saddle-point mani- 
fold) has the form: 



7^ = E • / - (V2r)Q 



(2.15) 



where / is the unity matrix, S is certain constant, and 
the 4x4 supermatrix Q — T^^AT satisfies the condi- 
tion = 1 and belongs to the cr-model target space 
described above. Finally, one performs the gradient ex- 
pansion of the second term in (|2.13p . for TZ having the 
form (|2.15p with a slowly varying Q(r). The expression 
for the propagator 11, Eq. (|2.8p . then reads. 



(2.16) 



n(ri,r2;c^) = J DQ Qll{v^)Q'Jiir^)e-'^Q\ 
where S[Q] is the cr-model action 

S[Q] = ^ J d'^r Str [-D{\/Q)^ - 2iujAQ], 



(2.17) 



The size 4 of the matrix is due to (i) two types of the 
Green functions (advanced and retarded) , and (ii) neces- 
sity to introduce bosonic and fermionic degrees of free- 
dom to represent these Green's function in terms of a 
functional integral. The matrix Q consists thus of four 
2x2 blocks according to its advanced-retarded structure, 
each of them being a supermatrix in the boson-fermion 
space. In particular, is the boson-boson element 
of the RA block, and so on. One can also consider an 
average of the product of n retarded and n advanced 
Green functions, which will generate a cr-model defined 



on a larger manifold, with the base being a product of 
U(n,n)/U(n) x U(n) and U(2n)/U(n) x U(n) (these are 
the same structures as in the replica formalism, but now 
without the n — limit). 

For other symmetry classes, the symmetry of the cr- 
model is different but the general picture is the same. 
For example, for the orthogonal class (AI) the 8 x 8 Q- 
matrices span the manifold whose base is the product 
of the non-compact space 0(2,2)/0(2) x 0(2) and the 
compact space Sp(4)/Sp(2) x Sp(2). The cr-model sym- 
metric spaces for all the classes (Wigner-Dyson as well 
as unconventional) are listed in Sec. IIVI 



2. RG \n 2 + e dimensions; e-expansion 

The cr-model is the effective low-momentum, low- 
frequency theory of the problem, describing the dynamics 
of interacting soft modes - diffusons and cooperons. Its 
RG treatment yields a flow equation of the form (|2.2p . 
thus justifying the scaling theory of localization. The 
/3-function P{t) = —dt/d\nL can be calculated pertur- 
batively in the coupling constant t inversely proportional 
to the dimensional conductance, t — l/2TTg? This al- 
lows one to get the e-expansion for the critical exponents 
in 2 + e dimensions, where the transition takes place at 

<C 1. In part icular, for the orthogonal symmetry class 
(AI) one finds (|Wegneij . Il989t) 



(i{t) ^et- 2t^ - 12C(3)t'5 -t- 0{t^). 



(2.18) 



The transition point is given by the zero of the /3- 
function. 



i* = |-^C(3)6^ + 0(e^)- 



(2.19) 



The localization length exponent v is determined by the 
derivative 

= -l//3'(i.) = - ^C(3)e' + 0{e^), (2.20) 

and the conductivity exponent s is 

s = j/e = 1 - ^C(3)e^ + 0(e4). (2.21) 

Numerical simulations of localization on fractals with 
dimensionality slightly above 2 give the behavior 
of v that is in good agreement with Eq. (|2.20[) 
( Schreiber and Grussbaclil [l9 96) . For the unitary sym- 
metry class (A), the corresponding results read 



I3{t) et-2t^ -6t^ + 0{t'^); 



(2.22) 
(2.23) 



s^^-^e + Oie'). (2.24) 



^ For spinful systems, g here does not include summation over spin 
projections. 
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In 2D (e = 0) the fixed point in both cases becomes 
zero: the /3-function is negative for any t > 0, implying 
that all states are localized. The situation is qualitatively 
different for the third Wigner-Dyson class - the symplec- 
tic one. The corresponding /3-function is related to that 
for the orthogonal class via /3sp(t) — — 2/3o(— yicld- 



mg'' 



l3{t) =et + t^ 



(2.25) 



In 2D the /3-function (|2.25p is positive at sufficiently small 
t, implying the existence of a truly metallic phase at 
t < tir, with an Anderson transition at certain ~ 1, 
Sec. IVI.BI This peculiarity of the symplectic class rep- 
resents one of possible mechanisms of the emergence of 
criticality in 2D, see Sec. IVI.AI The results for the (3- 
functions in all the symmetry classes will be given in 
Sec HVn 



3. Additional comments 

A few interrelated comments are in order here. 

(i) e-expansion vs 3D exponents. The e-expansion is of 
asymptotic character, yielding numerically accurate val- 
ues of the critical exponents only in the limit of small e. 
It is thus not surprising that, if Eqs. ()2.20p . ()2.2ip are 
used to estimate the indices in 3D (e = 1), the best thing 
to do is to keep just the leading (one-loop) term, yielding 
a quite substantial error (v = 1 instead ofv c± 1.57 ±0.02 
( Slevin and Ohtstikil . Il999[ ) known from numerical simu- 
lations). Rather unexpectedly, the agreement turns out 
to be remarkably good for the multifractal exponents, see 
Sec. III. C. 61 Independently of the accuracy of these num- 
bers, the e-expansion plays a major role in understanding 
of qualitative properties of the transition. 

(ii) Composite operators. The RG can also be used to 
calculate scaling dimensions of composite operators. In 
particular, the operators of the type (QA)" determine 
multifractal fluctuations of wave functions at criticality 
that will be discussed in Sec. III. Cl Another class of opera- 
tors - those with high de r ivatives, (VQ)^ " - were studied 



in lKravtsov et all 



iLerner and Wegneil ( 1990f) and 



found to have the scaling dimensions 

yn^d-2n + 2t^n{n - 1) + 0{tl), 



(2.26) 



where t^ is given by Eq. (|2.19p for the orthogonal class 
and Eq. (|2.23p for the unitary class. The fact that 
the one-loop result (|2.26p becomes positive for a suf- 
ficiently large number of gradients, n > t~^ (suggest- 
ing that the corresponding operators are relevant and 
might drive the system into an unknown fixed point) has 



launched a debate about the stability of the cr-model and 
the one-parameter scaling. This question is not specific 
to the localization problem but is equally applicable to 
a broader class of cr-models, including t he 0(n) model 
of a Heisenber g ferro magnet dCastilla and Chakravartvi 



119931: IWegneil. ll990D It was, however, pointed out 



^ Here t = I/ttj, where g is the total conductance of the spinful 
system. 



i n iBrezin and Hikamil ( 19971 ): iDerkachov and Manashovl 
(119971 ) that since the expansion (|2.26p is asymptotic and 
its true parameter is t*n, the behavior of the one-loop 
result at i:,,n > 1 does not allow to make any reliable 
conclusions. 

(iii) Order parameter. In view of the analogy with con- 
tinuous thermodynamic phase transitions, it is natural to 
ask what is the order parameter for the Anderson tran- 
sition. While naively Eq. (|2.17p suggests that it is the 
expectation value of Q, the latter is in fact uncritical. 
To describe the transition in terms of symmetry break- 
ing, one has to introduce an order parameter function 
(OPF) F{Q) resulting from integrating out Q-fiel ds at all 
points except one (rp). with given g(ro) = Q (|Efetovl . 
119871 : IZirnbaueil Il986allbl ). One can introduce an OPF 
i^($) with similar properties also within the supervector 
formalism (iMirlin and Fvodorovi Il9^li ). It was shown 
(jMirlin and Fvodorovi Il994aibl ) that the OPF is closely 
related to the distribution of one-site Green functions, in 
particular local density of states (LDOS), and wave func- 
tion amplitudes. In the framework of scattering theory, 
this suggests an inerpretation of the Anderson transition 
as a phenomenon of spontaneous breakdown of S*— matrix 
unitarity (Fyodorgv, 2003). 

(iv) Upper critical dimension. For conventional criti- 
cal phenomena, there exists an upper critical dimension 
dc above which the transition is governed by a Gaus- 
sian fixed point, with exponents being d-independent and 
given by their mean-field values. As a consequence, an 
e-expansion near dc (in the most standard case, in 4 — e 
dimensions) exists, alternative to 2 -t- e expansion. One 
can ask whether this is also the case for the Anderson 
localization transition. The answer is negative: there 
is no conventional mean field theory for the Anderson 
transition, and it was ar gued that the the upper criti cal 
dimension is dc = co (|Mirlin and Fvodorovi Il994al ibl). 
The closest existing analog of the mean-field theory is 
the model on the Bethe lattice corresponding to d = oo; 
we will discuss it in more detail in Sec. III.D] 



C. Critical wave functions: Multifractality 

1. Scaling of inverse participation ratios 

Multifractality of wave functions, describing their 
strong fluctuations at criti cality, is a striking feature 
of the Anders on transitions ( Castellani and Pelitil . [l986l : 
I Wegneil . 119801). Multifractality as a concept has been in- 
troduced by Mandelbrod ( 1974i ). Multifractal structures 
are characterized by an infinite set of critical exponents 
describing the scaling of the moments of some distribu- 
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tion. This feature has been observed in various complex 
objects, such as the energy dissipating set in turbulence, 
strange attractors in chaotic dynamical systems, and the 
growth probability distribution in diffusion-limited ag- 
gregation. For the present problem, the underlying nor- 
malized measure is just |V'^(r)| and the corresponding 
moments are the inverse participation ratios (IPR) * 



(2.27) 



At criticality, Pq show an anomalous scaling with the 
system size L, 



(2.28) 



governed by a continuous set of exponents Tq. One often 
introduces fractal dimensions Dq via Tq — Dq{q — 1). In a 
metal Dq = d, in an insulator Dq =^ 0, while at a critical 
point Dq is a non-trivial function of q, implying wave 
function multifractality. Splitting off the normal part, 
one defines the anomalous dimensions Aq, 



d{q-l) + Aq, 



(2.29) 



which distinguish the critical point from the metallic 
phase and determine the scale dependence of the wave 
function correlations. Among them, A2 = —rj plays the 
most prominent role, governing the spatial correlations 
of the "intensity" 

L2'i(|V,2(r)V'2(r')|)~(|r-r'|/i)-''. (2.30) 



Equation ((2?30)) can be obtained from (|2.28[) by us- 
ing the fact that the wave function amplitudes be- 
come essentially uncorrelated at |r — r'| ^ L. 
Scaling behavior of higher order spatial correlations, 
(|?/'^«i(ri)V'^'"(r2) . . . V^''"(rn)|) can be found in a sim- 
ilar way, e.g. 



i'i(9i+92)(|^2gi^j.^)^292(j.2)|) 



i-^"-^-(|ri-r2|/L)^ 



'91+52 ^91 ^92 



Correlations of two different (but close in energy) eigen- 
functions possess the same scaling properties. 



i2'i(V,,(r)^*(r)V*(r')^,(r')> 



|r — r 



(2.32) 

where uj = ei — ej, ^ {puj)^'^/'^^ p is the density 
of states, and |r — r'| < L^. For conventional classes. 



* Strictly speaking, Pq as defined by Eq. 12.271 1. diverges for suffi- 
ciently negative q {q < —1/2 for real ip and q < —3/2 for complex 
?/)), because of zeros of wave functions related to their oscillations 
on the scale of the wave length. To find Tq for such negative q, 
one should first smooth \tf>'^\ by averaging over some microscopic 
volume (block of several neighboring sites in the discrete version). 



where the DOS is uncritical, the diffusion propagator 
(|2.8p scales in the same way. 

In the field-theoretical language (Sec lILBjl . are the 
leading anomalous dimensions of the operators Tr{QA)'^ 
(or, more generally, Tr(Q A)i^ . . . Tr(gA)«'" with qi+. . .+ 
q) (|Wegneil . 11980) . The strong multifractal fluctu- 
ations of wave functions at criticality are related to the 
fact that Ag < for q > 1, so that the corresponding op- 
erators increase under RG. In this formalism, the scaling 
of correlation functions [Eq. (|2.30p and its higher-order 
gener alizations] results from an operator product e xpan- 
sion (iDuplantier and Ludwid . Il99ll; iMudrv et cdl . 1 19961: 
IWegneil . ll985l^ 



2. Singularity spectrum f{a) 

The average IPR {Pq) are (up to the normalization fac- 
tor L'^) the moments of the distribution function T'dV'P) 
of the eigenfunction intensities. The behavior (|2.28p of 
the moments corresponds to the intensity distribution 
function of the form 

^ (2.33) 

Indeed, calculating the moments (IV'^'^I) with the distri- 
bution function (|2.33p . one finds 



(Pq) = L'^(|i/;29|) ^ / daL-9"+/(°) 



(2.34) 



where we have introduced a = — In In L. Evaluation 
of the integral by the saddle-point method (justified in 
the limit of large L) reproduces the result p.28p . with 
the exponent Tq related to the singularity spectrum /(a) 
via the Legendre transformation. 



T{q) ^ qa~ f{a) , q = ,f'{a) 



(2.35) 



The meaning of the function /(a) is as follows: it is the 
fractal dimension of the set of those points r where the 
(2.31) eigenfunction intensity is |V'^(r)| ~ In other words. 



in a lattice versi on of the model the number of such points 
scales as L^^") (|Halsev a/.l . ll986f) . 

General properties of the functions Tq and /(a) follow 
from their definitions and the wave function normaliza- 
tion: 

(i) Tq is a non-decreasing, convex function (t^ > 0, 
T^' < ), with To = ~d, n = 0; 

(ii) f{a) is a convex function (/"(a) < 0) defined on 
the semiaxis a > with a maximum at some point uq 
(corresponding to g = under the Legendre transforma- 
tion) and f{ao) = d. Further, for the point ai (corre- 
sponding to q — 1) we have /(ai) — a\ and j'(pi\) = 1. 

If one formally defines f{a) for a metal, it will be con- 
centrated in a single point a — d, with f(d) = d and 
/(a) = —00 otherwise. On the other hand, at critical- 
ity this "needle" broadens and the maximum shifts to a 
position ao > d, see Fig. [TJ 



d a 

FIG. 1 Schematic plot of the multifractal spectrum, f{a). 
A metal is represented by a "needle", i.e. f{a) having zero 
width, at Q = d. At criticality /(a) acquires a finite width 
and the apex, qq, shifts to a value larger than d. The negative 
parts of /(a) (grey area) correspond to rare events - values 
of the wavefunction amplitude that typically do not occur in 
a single sample. 

3. Weak multifractality: approximately parabolic spectrum 

One of the situations in which the Tq spectrum can be 
evaluated analytically is the regime of weak multifractal- 
ity, when the critical point is, in a sense, close to a metal. 
This happens, in particular, for the Anderson transition 
in 2 + e dimensions with e <C 1, see Sec. lILCTBl and in the 
PRBM model with 6 » 1, Sec. IIILBI In this situation, 
one finds generically a spectrum of the form 

~d(g-l)-7(7(g-l), 7 « 1, (2.36) 

i.e. the anomalous dimension ~ 7'?(1 — q)- (We re- 
mind that Ao = Ai = by definition.) The approxima- 
tion (|2.36p is valid in general as long as the second term 
(Aq) is small compared to the first one, i.e. for q <C d/j. 
After the Legendre transformation Eq. (|2.36p yields 

/W^rf- I," ao-d + 7. (2.37) 
4(q;o - d) 

In some specific cases, the parabolic form of the spec- 
trum (|2.36p . (|2.37p is not just an approximation but 
rather an exact result. This happens, in particular, 
for the random vector potential model, see Sec. IVI.G.3I 
There is a conjecture corroborated by numerical simula- 
tions that this is also the case for the IQHE transition, 
see Sec. I VI. C. 71 Note that exact parabolicity cannot ex- 
tend to all q: at qc — {d+'y)/2^ the derivative becomes 
zero (i.e. the corresponding a = 0) , so that Tq should stay 
constant for larger q. We will discuss this issue, known as 
"termination" of the multifractal spectra, in Sec. III. C. 71 



FIG. 2 Multifractal exponents A, for the PRBM model with 
6 = 4, 1, 0.3, 0.1. The symmetry (fOS)) with respect to 
the point q = 1/2 is evident. A small difference between 
Aq (full line) and Ai-g (da shed) is due to numerical errors. 
(jMildenberger et fflU[2007bl ). 
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FIG. 3 The data of Fig.[2]in terms of the singularity spectrum 
/(q). Dashed lines represe nt /(2 — a) -|- a — 1, demons trating 
the validity of Eq. (|Mildenberger et aU UoOTbl ) . 

an exact symmetry relation^ 

A, = Ai_, , (2.38) 

connecting exponents with g < 1/2 (in particular, with 
negative q) to those with q > 1/2. In terms of the singu- 
larity spectrum, this implies 

/(2d - a) = /(a) -I- d - a. (2.39) 

The analytical derivation of Eqs. (|2.38p . (|2.39p is based 
on the supersymmetric cr-model; it has been confirmed 
by numerical simul ations on the P RBM model at crit- 
icality (Mildcnbcrg er tt all l2007bf) . see Figs. [2 [3] and 
Sec. Illli and the 2D Anderson tran si tion o f the symplec- 
tic cla ss f|Mildenberger and Eversl l2007t lObuse tt al\ . 
l2007bD . Secini 

5. Role of ensemble averaging 

a. Average vs. typical spectra. It should be stressed that 
the definition (I2.28P of Tq is based on the ensemble- 



4. Symmetry of the multifractal spectra 



As was recently shown ( Mirlin et all |2006|) , the mul- 
tifractal exponents for the Wigner-Dyson classes satisfy 



^ If the multifractal spectrum possesses a termination (non- 
analyticity) point qc, Sec. III. C. 71 the status of the relation ||2.38| I 
beyond this point is not clear. 
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averaged IPRs, (Pq). On the other hand, until recently 
most numerical studies of multifractality were dealing 
with properties of a single (representative) wave func- 
tion. Formally, this corresponds to an analysis of the 
typical IPR, 

P*yP = exp(lnP,). (2.40) 
Similarly to Eq. (|2.28p . one can define the exponents t^^^, 



ptyp 
g 



L- 



(2.41) 



and introduce the spectrum /*yp(a) as the Legendre 
transform of r^^P. The relation between Tq, f{a), on one 
side, and t^^ p, P^Pja), on the other s ide, was analyzed 
i n det ail in lEvers and Mirlii] ( 2000l ): iMirlin and Ever j 



(I200GD .6 The function r^^P has the form 




q<q^ 

q- <q <q+ 

q> q+ ■, 



(2.42) 



where a± are determined by the condition /(a) = 0, and 
q± are the corresponding values of q, with q^ < q+7 
The singularity spectrum /*yp(a) is defined on the in- 
terval [Q:_|_,a_], where it is equal to f{a). The informa- 
tion about the negative part of /(a) (on a < and 
a > a-), or, equivalently, about the part of Tq with q 
outside the range g+J, gets lost when one considers a 
single wave function. This is because the average number 
of points with such a singularity a for a single eigenfunc- 
tion is L-^^") <C 1, so that the ensemble averaging is of 
crucial importance for determination of this part of the 
multifractal spectrum, see Fig. [1] 



The upper cutoff of this tail, {Pq/ Pq^^)niax depends on 
the system size L, moving to infinity with L —f oo. It 
is clear that the relation between t^^p and Tq depends 
crucially on the power-law exponent Xq. li Xq > 1, 
the two definitions of the fractal exponents are identi- 
cal, T{q) — T'yP(g). On the other hand, if Xq < 1, the 
average (Pq) is determined by the upper cut-off of the 
power-law tail, which depends on L. As a result, (Pq) 
shows scaling with an exponent r*^P different from Tq . In 
this situation the average value (Pq) is not representative 
and is determined by rare realizations of disorder. Thus, 
Xq = 1 ioT q = q±, Xq > 1 for q- < q < q+, and Xg < 1 
other wise. Furthermore, it was found ( Mirlin and Everi . 

that the power-law-tail index Xq is related to the 
fractal exponents as follows: 



rr. T-typ 

■^q I q 



iqxq 



(2.44) 



More precisely, Eq. (|2.44p was proven for the case when 
Xq is an integer. Also, it was shown to hold for the small- 
b limit of the PRBM model, at q > 1/2. The generic 
validity of this formula remains a conjecture. In the range 
of non-self-averaging IPR it yields 



Xq=q+/q, q>q+: 



(2.45) 



and similarly for q < q^. As to the range q- < q < q-^., 
the behavior of Xq depends on the specific form of Tq. In 
the particular case of the weak multifractality, Eq. (|2.36p . 
the solution of Eq. (|2.44p reads 



Xq ^ iq+/qf 



q- <q <q+ , 



with 



(2.46) 



(2.47) 



b. IPR distribution and tail exponents. A closely related 
issue is that of the distribution function of the IPR 



Pq. It was c onjectured in Fvodorov and Mirliiil (Il995 ) 



and shown in lEvers and MirlinI (|200o( )7 iMirlin and Everg 
(l2000h that the distribution function of the IPR normal- 
ized to its typical value P^^p has a scale invariant form 
V{Pq/Pq^^) at criticality. In other words, the distribu- 
tion function of the IPR logarithm, VihiPq), preserves 
its form and only shifts along the x axis with increas- 
ing L. On the large- side, this distribution develops a 
power-law tail, 

V{Pq/Pq"') OC (Pg/Pf P)-l-^', Pq » Pf P . (2.43) 



6 In lEvers and MirlinI ||2000D : IMirlin and Ever^ l l2000D different no- 
tations were used: Tq, /(a) for the typical spectra, and Tq, f{a) 
for the averaged spectra. 

It is tacidly assumed here, that q± and Q- actually exist, i.e. 
are not infinite. To the best of our knowledge, an example to 
the opposite has never been encountered in the context of the 
Anderson transitions. 



6. Dimensionality dependence of the wave function statistics at 
the Anderson transition 

In this subsect i on, which is largely based on 
iMildenberger et oD ( 20021 ). we summarize the results for 
the Anderson transition in d dimensions, obtained by an- 
alytical and numerical means. This allows us to ana- 
lyze the evolution of the critical statistics from the weak- 
multifractality regime in d = 2 -I- e dimensions to the 
strong multifractality at 1. 

In 2-f e dimensions with e ^ 1 the multifractality expo- 
nents can be obtained within the e-expansion, Sec. III.BT2] 
The 4-loop results for the orthog onal and unitary sym- 
metry classes read ( Wegneil . Il987 ) 



Af) = 9(1 - q)e + ^q{q - l){q^ - q + l)e' 



+0(e5); (2.48) 



Al^^^q{l-q)ie/2Y/'-lq'iq-ini3)e 



+0(e5/2-)_ ^2.49) 
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FIG. 4 IPR distribution at the Anderson transition (GOE) 
(a) 3D (system sizes L = 8,11,16,22,32,44,64, 80) and (b) 
4D (L = 8, 10, 12, 14, 16). l|Mildenberger et all 1200^ 



Keeping only the leading (one-loop) term on the r.h.s. of 
Eqs. (|2.48|) and (|2.49p . we get the parabolic approxima- 
tion for T,, Eq. (12311), and f{a), Eq. (EST]), with 7 = e 
for the orthogonal class and 7 = (e/2)^/'^ for the unitary 
class. The IPR fluctuations (Sec. III.C.5|) can be studied 
analytically as well. In particular, in the orthogonal sym- 
metry class, the variance of the distribution P(lnPg) 
is given, to the leading order in e <C 1, by 



aq = Sn^ae^q^iq-l) 



\q\ < q+ , 



(2.50) 



where a ~ 0.00387 for the periodic boundary conditions. 
The values q± of q beyond which the typical and the 
average IPR scale differently are given by Eq. (|2.47p with 
d ~ 2 and the above values of 7. The power-law exponent 
Xq of the IPR distribution is given by Eqs. ^J5\), ^M\). 
At g ^ g-f the variance aq is governed by the slowly 
decaying power-law tail, yielding 

aq ~ Xq'^ ~ q/q+ . (2.51) 

Results of numerical simulations of the wave function 
statistics in 3D and 4D for the orthogonal symmetry 
class are shown in Figs. HHTJ in comparison with the 
one-loop analytical results of the 2 + e expansion. Fig- 
ure m demonstrates that the critical IPR distribution 
P(l nPn) acquires the s cale- invariant form (as also found 
m ^ (|Cuevas et all |2002[ )). The corresponding variance is 
shown in Fig. [SJ in 3D it is described very well by the 
analytical formulas with e — I. The evolution from the 
weak- to strong-multifractality with increasing d is nicely 
seen in Fig. [5] for /(a). As is demonstrated in the inset, 
the one-loop result of the 2 + e expansion with e = 1 
describes the 3D singularity spectrum with a remarkable 
accuracy (though with detectable deviations). In partic- 
ular, the position of the maximum, ao — 4.03 ± 0.05, is 
very close to its value ao — d + e implied by Eq. (|2.37p . 
As expected, in 4D the deviations from parabolic shape 
are much more pronounced and ao = 6.5 ± 0.2 differs 
noticeably from 6. 



Evolution of the fractal dimension 



rj{q - 1) 




FIG. 5 The rms deviation Oq of InP, extrapolated to L ^ 00 
in 3D (x) and 4D (*). The dotted line is the analytical result 
(f230| for e = 0.2; the full lines represent Eqs. (fl^SOjl . ([ZSTjl 
with e = 1. Inset: evolution of cr, with L in 3D for values of 
q = 0.5, 1.5, 2, 2.5, 3, 4, 5, 6. The leading finite size correction 
of all data has the fo rm with y = 0.25 .5 for 3D and 
y = 0.1 0.4 in 4D. (jMildenberger ef all l2002t ) 




FIG. 6 Singularity spectrum /(a) in 3D (dashed) and 4D (full 
line). To illustrate the evolution of the spectrum from d — 2 
to d — 4, analytical results for d = 2 + e are shown for e = 
0.2 (dotted) and e = 0.01 (dot-dashed). Inset: comparison 
between f{a) for 3D and th e one-loop result of the 2 -I- e 
expansion with e = 1 (solid). (|Mildenberger et all [200I ) 



dimensions Dq with g ^ 1 decrease with increasing d. 
As an example, for g = 2 we have D2 ~ 2 — 2e in 2 -I- e 
dimensions, D2 = 1.3±0.05 in 3D, and D2 = 0.9±0.15 in 
4D. This confirms the expectation based on the Bethe- 
lattice results (Sec. III.D[) that ^ at d ^ 00 for 
q > 1/2. Such a behavior of the multifractal exponents 
is a manifestation of a very sparse character of critical 
eigenstates at d ^ 1, formed by rare resonance spikes. 
In combination with the symmetry relation (j2.38p this 
implies the limiting form of the multifractal spectrum at 
d 00, 



with d is shown in Fig. [71 It is seen that the fractal 



, 
2d{q 



1/2) 



9>l/2 
q<l/2 



(2.52) 
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FIG. 7 Fractal dimensions Dq/q{q — 1) in 3D (dashed) 
and 4D (full line). Analytical results for d = 2 + e with 
e = 0.2 (dotted) and e = 0.01 (dot-dashed) are also shown. 
(jMildenbereer et all [200^ ) 



This corresponds to /(a) of the form 

f{a) = a/2 , < a < 2rf , (2.53) 

dr opping to — oo at th e boun daries of the interval [0, 2d]. 
In iMildenberger et al\ ( 20021 ) arguments were given that 
the way the multifractahty spectrum approaches this lim- 
iting form with increasing d is analogous to the behavior 
found in the PRBM model with 6 < 1, Sec. UlLCl 



7. Possible singularities in multifractal spectra: termination 
and freezing 

In this subsection, we discuss what kinds of singular- 
ities may be typically encountered in the multifractahty 
spectra /(a) and Tq. First of all, we recall that the spec- 
trum r*''P of a typical eigenfunction has non-analyticity 
points at q±, corresponding to the termination of /'^p(q;) 
at its zero a±, see Sec. III. C. 51 However, the ensemble- 
averaged spectra Tq and f{a) (that we are considering 
throughout) do not have any singularity there. 

Singularities in Tq and /(a) may arise, depending on 
the behavior of /(a) at a = in the particular critical 
system under investigation. One possibility is that f{a) 
approaches the a = axis continuously, with f{a) 
—oo as a ^ (Fig. [5^). Then Tq increases monotonically 
with q, without any non-analyticities. Such a situation 
is realized e.g. in the PRBM model, see Sec. IIII.CI An 
alternative option is that /(O) is finite, see Fig. [SJj. This 
generically implies that Tq has a discontinuity in the sec- 
ond derivative at certain Qc = f'{a)\a~^o and is strictly 
constant, Tq = — /(O) at q > qc- Such a behavior of 
the multifractahty spectrum at q — qc is called "termi- 
nation" . In particular, it takes place unavoidably if the 
spectrum is exactly parabolic, as is the case, e.g., for the 
random vector potential problem, Sec. IVI.G.3I From the 
point of view of the underlying field theory, termination 
implies that there is a qualitative change in properties of 
the operators Oq describing the moments of the density 
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FIG. 8 Possible behavior of the singularity spectrum /(a) at 
a ^ 0: (a) no singularity, (b) termination, (c) freezing. The 
corresponding behavior of is shown as well. 



of states. An explicit example of how this may happen is 
provided by the 2D Liouville field theorv llKogan et al. , 
1996 : [Se"iberg| . ll990l : IZamolodchikov and Zamolodchikov . 
1996| ) (closely related to the random vector potential 
problem), where the operators Oq cease to be local for 
q > Qc- 

A spectrum with termination may show another pe- 
culiarity. While normally /(O) is negative, one can also 
imagine a situation with /(O) = 0, Fig.[5J; (corresponding 
to <7c < !)■ In fact, this is exactly what happens in the 
random vector potential problem. Sec. IVI.G.31 when the 
disorder strength exceeds a certain critical value. The 
transition into this phase is termed "freezing transition" . 
In the "frozen" phase the wave functions combine proper- 
ties of localized and critical states: while the wave func- 
tion normalization is governed by a vicinity of one or few 
(of order unity) points, the tails away from these points 
show multifractal fluctuations and correlations. 



8. Surface vs. bulk multifractality 

Recently, th e concept of wave function multifractality 
was extended ( Subramaniam et all . [20061 ) to the surface 
of a system at the critical point of an Anderson transi- 
tion. It was shown that the fluctuations of critical wave 
functions at the surface are characterized by a new set 
of exponent (or, equivalently, anomalous exponents 
Ap, which are in general independent from their bulk 
counterparts. 



T'q=d{q-l) + qfl+l + Al. 



(2.54) 
(2.55) 



Here /z is introduced for generality, in order to account 
for a possibility of non-trivial scaling of the average value. 
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(|'0(r)P) oc L~'^~'^, at the boundary in unconventional 
symmetry classes. For the Wigner-Dyson classes, /i — 
0. The normalization factor L"^^^ is chosen such that 
Eq. (|2.54p yields the contribution of the surface to the 
inverse participation ratios (Pq) = {J d'^r\tp{r)\'^'^). The 
exponents as defined in Eq. ()2.55|) vanish in a metal 
and govern statistical fluctuations of wave functions at 
the boundary, (|V'(r)|2'3)/(|V'(r)|2)9 ~ L'^l, as weh as 
their spatial correlations, e.g. L^(''+^)(|'0^(r)-0^(r')|) ^ 

(|r-r'|/i)^'=. 

Wave function fluctuations are much stronger at the 
edge than in the bulk. As a result, surface exponents are 
important even if one performs a multifractal analysis 
for the whole sample, without separating it into "bulk" 
and "surface" , despite the fact that the weight of sur- 
face points is down by a factor 1/L. This was an alyt- 
ically demonstrated in ISubramaniam et all ( 2006t ). us- 
ing a model of a 2D weakly localized metallic system 
(large dimensionless conductance g ^ 1), which shows 
weak multifractality on length scales below the localiza- 
tion length ^ ~ e^'^^-' , where /3 = 1 (2) for systems with 
preserved (resp. broken) time-reversal symmetry. With 
minor modifications, the formulas below describe also the 
Anderson transition in 2 -I- e dimensions. 

For the bulk multifractal spectr um one gets the result 
(lOel) with 7 = ( Pttq)'} < f 1 (lAltshuler et all . Il986t 



iFal'ko and Efetovl . il995allbl : |Wegner|, llSSOD; generaliza- 
tion of this result to the surface case reads: 

r^^ = 2(g-l) + l + 279(1 -g). (2.56) 

The corresponding /(a)-spectra have the form: 

/^a) = 2 - (a - 2 - 7) V47, (2.57) 
r{a) = 1- (a- 2 -27)2/87. (2.58) 



These results are illustrated in Fig. [51 When the mul- 
tifractality in the whole sample is analyzed, the low- 
est of the Tq exponents "wins" . The surface effects be- 
come dominant outside the range < q < q^, where 
qhs ^ ±^-1/2 are the roots of the equation — t^- The 
lower panel of Fig. [5] shows how this is translated into the 
f{a) representation. The total singularity spectrum is 
given by the bulk function f^{a) only for < a < a'i, 
where a\ — 2 ~ +27^/^. Outside this range the sur- 
face effects are important. Specifically, /(a) is equal to 
the surface spectrum P{a) for a < and a > aL, 
where a!^ — 2 ~ +47^/^, while in the intermediate inter- 
vals q;^ < a < and < a < at_ its dependence 
on a becomes linear (shown by dashed lines). The latter 
behavior is governed by intermediate (between "bulk" 
and "surface") points with a distance from the surface 
r ~ , < /3 < 1; their /(a) spectrum is found to 
be fpia) = Pf^ia) + (1 - P)r{a). Note that in this 
case the surface effects modify f{a) in the whole range 
below /(a) ~ 1. Therefore, the surface exponents affect 
the multifractal spectrum of the sample not only for rare 
realizations of disorder (governing the negative part of 
/(a)) but already in a typical sample. 
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FIG. 9 Surface and bulk multifractal spectra Tq and /(a) 
for a 2D metal with 7 = 0.01. For details see text. 
(|Subramaniam et aZ.I . |2006| ). 

The boundary multifractality was also explicitly stud- 
ied, analytically as well as numerically, for several other 
syste ms at criticality: the 2D s pin quantum Hall transi- 
tion ( Subramaniam tt aLl . [200^ . the Ander son transition 
in a 2 D system with spin-orbit co upling (lObuse et al. 
2007bl) . and the PRBM model (!Mil denberger et al. 



2007bl), see Sec. IVI D.5I IVLB:! andHlLEl respectively. 

In lObuse et all ( 2007bf ) the notion of surface multi- 
fractality was further generalized to a corner of a critical 
system. It was shown that for a 2D system at criticality 
conformal invariance leads to the following dependence of 
the corresponding anomalous exponent on the open- 
ing angle Q, 



(2.59) 



More carefully, for 9 < n the spectrum will terminate 
at some qg (see Sec. III.C.71 Fig. [SId), even if the surface 
spectrum showed no singularity, as in Fig. [S^. Equation 
(|2.59p holds then for q < qg; for larger q the exponent Tq 
is constant. 



9. Manifestations of multifractality in other observables. 

The multifractal structure of wave functions at crit- 
icality manifests itself also in other physical character- 
istics of the system. In particular, one can open the 
system by attaching a local lead at some point r. The 
system can then be characterized by the Wigner de- 
lay time tw (energy derivative of the scattering phase 
shift), whose statistical properties in chaotic and dis- 
ordered systems have b ecome a subject of resea r ch ac- 
tivity i n rec ent years ( Fvodorov and Sommersl Il997t 
iKottod . |2005[ ). At criticality, the mo ments of the in 



verse delay time show a scaling behavior (iFvodoroyl 



Mendez-Bermudez and Kottosll2005l:lMirhn et all . 
Ossipov and Fvodorov . 2005t) 



2003; 



20061: 



it-' 



w / 



(2.60) 
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It was shown that for all Wigner-Dyson classes the 
exponents 7^ are linked to the wave functi o n ex- 
ponents by an exact re lation ( Mirlin et al\ . 120061 : 
lOssipov and Fvodorovl [20051) . 



Iq = Tl+q- 



(2.61) 



If a second local lead is attached to the system, the 
statistics of the two-point conductance (7(r,r') can be 
studied. Specifically, one can analyze the scaling of the 
moments (q'^(r^r)) with the distance |r — r^| between the 
contacts (|Janfien et ali Il999l : IZirnbaueii Il994 Il999l ) , 



(5'(r,r')) 



(2.62) 



For the case of the unitary symmetry class (A), a re- 
lation linking the exponents Xq to the wave function 
anomalous dimensions [and based on a result of 



Klesse and Zirnbaued ( 20oil )] was obtained ( Evers et al\ . 
200lh. 



Aq + Al_ 

2A1/2 , 



g< 1/2 
q > 1/2. 



(2.63) 



In view of l|2.38p . the first line of Eq. (|2.63p can be equiva- 
len tly written as = 2 A^, which has also been proposed 
in (jJanBen et a/.l. Il999l) . 

A relation between the exponents Xq and A^ was also 
derived for the case of the spin quantum Hall transition 
fSec. lVI.D.Sp . belonging to the unconventional symmetry 
class C. In contrast to critical points of the Wigner-Dyson 
classes, the density of states p in this case is critical, i.e. 
it has a non-trivial scalin g dimensions p oc L ~^p with 



Xn > 



0. It was shown in [Mirlin et al. I (|2003[) (see also 



[Bernard and LeClairl ([2002bl )V that in this case 



Xq = 



2qxp + 2Aq 
X. 



qo J 



<i<qo 



(2.64) 



2Aq reaches its 



where go is the point at which 2qxp 
maximum. It is plausible that the relation (|2.64p (which 
reduces to (I2.63P for the Wigner-Dyson classes) holds in 
fact for critical points of all symmetry classes 



D. Anderson transition in d = 00: Bethe lattice 

The Bethe lattice (BL) is a tree-like lattice with a fixed 
coordination number. Since the number of sites at a dis- 
tance r increases exponentially with r on the BL, it ef- 
fectively corresponds to the limit of high dimensionality 
d. As has been already mentioned in Sec. III.BTsI the BL 
models are the closest existing analogs of the mean-field 
theory for the case of the Anderson transition. 

The Anderson tight-binding model (lattice version of 
Eqs. (1^. (1^) on the BL was studied for the first 



time in lAbou-Chacra et al\ ( 1973() . where the existence 
of the metal- insulator transition was proven and the po- 
sition of the mobility edge was determined. The ana- 
lytical results were confirmed by numerical simulations 



([Abou-Chacra and Thoulesi . [1974 [Girvin and Jonson[ . 
I1980D . In later works the BL versions of the cr-model 
(I^TTl) (lEfetovLll985L ll987HZirnbaueilll986a [ bl) and of the 

tight-binding model ( Mirlin and Fvodorovl [l99lf) were 
studied within the supersymmetry formalism, which al- 
lowed to determine the critical behavior. It was found 
that the localization length diverges in the way usual 
for BL models, ^ oc \E — Ec\~^, where E' is a micro- 
scopic parameter driving the transition. When rein- 
terpreted within the effective-med ium approximation 
(|Efetovl[T990l : [Fvodorov et aU[l992D . this yields the con- 
ventional mean-field value of the localization length expo- 
nent, ly = 1/2. On the other hand, the critical behavior of 
other observables is very peculiar. The inverse participa- 
tion ratios Pq with q > 1/2 have a finite limit at E Ec 
when the critical point is approached from the localized 
phase and then jump to zero. By comparison with the 



scaling formula, Pq 



this can be interpreted as 



Tq — for all q > 1/2. Further, in the delocalized phase 
the diffusion coefficient vanishes exponentially when the 
critical point is approached. 



D cx n-^ W n ; 

n ^ expjconst \E ~ E^ 



(2.65) 
(2.66) 



which can be thought as corresponding to the infinite 
value, s = 00, of the critical index s. The distribution 
function of the LDOS v = p{r)/{p) (normalized to its 
average value for convenience) was found to be of the 
form 



(2.67) 



and exponentially small outside this range. Equation 
(I2.67P implies the following behavior of the LDOS mo- 
ments: 



oc 



^|g-l/2|-l/2^ 



(2.68) 



The physical reason for the unconventional critical be- 
havior was unravelled in [Mirlin and Fvodorovl ( 1994a[ [b[) . 
It was shown that the exponential largeness of the fac- 
tor n reflects the spatial structure of the BL: the "cor- 
relation volume" (number of sites within a distance 
^ from the given one) on such a lattice is exponentially 
large. On the other hand, for any finite dimensional- 
ity d the correlation volume has a power-law behavior, 
Vd{0 (X oc \E - Ecl""^, where 1/ ~ 1/2 at large d. 
Thus, the scale cannot appear for finite d and, assum- 
ing some matching between the BL and large-d results, 
will be replaced by Vd{0- Then Eq. ((2:68)) yields the fol- 
lowing high-d behavior of the anomalous exponents Aq 
governing the scaling of the LDOS moments fSec. III.C)) . 



A„ 



d{l/2 



1/2|) , 



(2.69) 



or, equivalently, the results (|2.52p . (I2.53P for the mul- 
tifractal spectra Tq, f{a). These formulas describe the 
strongest possible multifractality. 
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in the Poisson limit {SN"^) = {N). At criticality, {5N'^) 
shows an intermediate Unear behavior 



FIG. 10 

symplectic 



2 3 4 

Level spmjti^ 8 

Critical level spacing distribution of 2D 
systems for periodic (PBC) and Dirichlet 
(DEC) boundary conditions and various system sizes 
L. PBC: L =20 (o); DEC: L=40 80, 120 (+,n,x). 
Adapted from ISchweitzer and Potempal (jl999i ): see also 
ISchweitzer and ZharekeshevI l|l997f r 



The critical behavior of the conductivity, Eq. (|2.65p . 
is governed by the same exponentially large factor fl. 
When it is replaced by the correlation volume Vd{^), the 
power-law behavior at finite d ^ 1 is recovered, a cx 
\E — Ecl'^ with s ~ d/2. The result for the exponent s 
agrees (within its accuracy, i.e. to the leading order in d) 
with the scaling relation s — v{d — 2). 



E. Level statistics at criticality 

We restrict ourselves here to a brief account of key 
results on the critical level statistics; more detailed ex- 

Jiosition and list of references can be found in the review 
Mirlin . ,2000b) . The primary quantity is the two-level 
correlation function (the superscript "c" standing for the 
connected part) 

R['\s)^{p)-^{p{E-Lo/2)p{E + uj/2))~l , (2.70) 

where p{E) ^ V-^Tt5{E-H) is the fluctuating DOS, V 
the system volume, s = uj/A, and A = l/{p)V the mean 

level spacing. In a metallic system i?2"^'' (s) is given, as a 
first ap proximation, by the RMT (Wigner-Dyson statis- 
tics; see lMirlinI (j2000bl ) for analysis of deviation and lim- 
its of applicability) , while in the insulating limit the levels 

are uncorrelated, R^\s) = 5{s) (Poisson statistics). In 
a critical point the level statistics takes an intermediate 
scale-invariant form. Specifically, R^\s) (and higher- 
order correlation functions) is independent under rescal- 
ing of the sample, although it does depend on the sam- 
ple shape and the boundary conditions, see Fig. [TOl A 
closely related quantity is the variance {SN{E)'^) of the 
number of levels N{E) in a spectral window of the width 
E, 



{dN{Er) = / ds{{N{E)) - \s\)Ri'\s) . (2.71) 



(NiE)) 
-{N{E)) 



In the RMT limit the level number variance increases log- 
arithmically, ((57V2) = (2/7r2/?)ln(7V) at {N) > 1, while 



{6N^) = x{N) 



(2.72) 



with a coefficient < % < 1 called "spectral compressibil- 
ity" . The parameter x is a universal characteristics of the 
critical theory, i.e. it has a status analogous to critical in- 
dices. Evolution of x from the "quasi-metallic" (x ^ 1) 
to the "quasi-insulating" (x close to 1) criticality can be 
explicitly analyzed for the family of critical PRBM the- 
ories. Sec. mil It is expected that similar evolution takes 
place for the Anderson transition in d dimensions when 
d is changed from c? = 2-|-etod^l, see a related dis- 
cussion of the wave function statistics in Sec. III.CT6l The 
"quasi-metallic" d = 2 + e limit can be studied analyti- 
cally with the result x = t* /P, where the critical coupling 
t* is given by Eqs. ()2.19p or ()2.23|) . depending on the 
symmetry class. The approach of the critical statistics 
to the Poisson limit at large d was clearly demonstrated 
in the recent numerical work (Garcia-Garcia and Cuevai, 
l2006f ) , where systems of dimensionality up to d = 6 were 
studied (with x reaching the value ~ 0.8 in 6D). 

In systems of unconventional symmetry classes already 
the one-point correlation function (average DOS) is non- 
trivial and acquires, in analogy with the two-level corre- 
lation function discussed above, a scale-invariant form at 
criticality. In particular, this will be shown in Sec. lVI.D.6l 
for the case of the SQH transition (class C). 



III. CRITICALITY IN THE POWER-LAW RANDOM 
BANDED MATRIX (PRBM) MODEL 



A. Definition and generalities 



The PRBM model is defined (jMirlin et all Il996l ) as 
the ensemble of random L x L Hermitean matrices H 
(real for /? = 1 or complex for (3 = 2). The matrix ele- 
ments Hij are independently distributed Gaussian vari- 
ables with zero mean (i/,, ) = and with variance 



\H,, 



|2\ _ 



where a(r) is given by 

a2(r) = [l + (r/6)2"] 



(3.1) 



(3.2) 



At a = 1 the model undergoes an Anderson transition 
from the localized (a > 1) to the delocalized {a < 1) 
phase. Below, we concentrate on the critical value a = 1, 
when a{r) falls down as a(r) cx l/r at r ^ b. 

In a straightforward interpretation, the PRBM model 
describes a ID sample with random long-range hop- 
ping, the hopping amplitude decaying as 1/r" with the 
distance. Also, such an ensemble arises as an effec- 
tive description i n a n umber of physical contexts, see 
lEvers and MirlinI (|2000l ) for relevant references. At a = 1 
the PRBM model is critical for arbitrary value of b and 
shows all the key features of the Anderson critical point. 
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including multifractality of eigenfunctions and non-trivial 
spectral compressibility. The existence of the parameter 
b which labels the critical point is a distinct feature of 
the PRBM model: Eq. p.ip defines a whole family of 
critical theories parametrized by b. The limit 6 3> 1 
represents a regime of weak multifractality, analogous 
to the conventional Anderson transition in d — 2 + e 
with e <C 1. This limit allows for a systematic analyti- 
cal treatment via a mapping ont o a supermatrix g-model 
and a w eak-cou p ling expansion (pavers and Mirlin j_ 200C ; 
Mirlinl . l2000bl : iMirlin and Eversl l2000l : iMirlin aL . 
19961 ). The opposite limit 6 <C 1 is characterized by 



very strongly fluctuating eigenfunctions, similarly to the 
Anderson transition in d ^ 1, where the transition 
takes place in the strong disorder ("strong coupling" 
in the field-theoretical language) regime. It is also ac- 
cessible to an analytical treatment using a real-space 
renormalization-group (RG) method ( Mirlin and Ever^ 
[2 000') introduced earlier for related models bv iLevitovl 

In addition to the feasibility of the systematic ana- 
lytical treatment of both the weak-coupling and strong- 
coupling regimes, the PRBM model is very well suited 
for numerical simulations in a broad range of couplings. 
For these reasons, it has attracted a considerable inter- 
est in the last few years as a model for the investigation 
of various properties of the Anderson critical point, see 
Sec. UlLEl HIlFI 



B. Weak multifractality, 6 > 1 

The q uasi-metallic regi m e b ^ 1 c an be 

studied (lEve rs and Mirliii'. '2000|; iMirlinl . l2000bl : 
iMirlin and Evera. ,200Q : .Mirlin et all 119961 1 via mapping 
onto the supermatrix a-model, cf. Sec. III. B. 11 



S[Q] = ^Str 

L rr' r 

(3-3) 

In momentum (k) space and in the low-fc limit, the action 
takes the form 



5[Q] =/3Str 



(3.4) 



where Qk = J2r e*'"'(5(?') and Q{r) is a 4 x 4 (/3 = 2) or 
8x8 {P — 1) supermatrix field constrained by Q^{r) — 1, 
see Sec. III.B.li p is the density of states given by the 
Wigner semicircle law 

p{E) = (l/2^2^)(4^fe - £;2) 1/2 ^ < , (3.5) 

and t ^ 1 is the coupling constant, 

l/t = (7r/4)(7rp)262 = (6/4)(l - E^/Awb) . (3.6) 



The main difference between the action (|3.4I) and that of 
the diffusive cr-model, Eq. (|2.17p is in the replacement of 
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FIG. 11 Fractal dimension D2 as a function of the parameter 
b of the PRBM ensemble. The data points are the results of 
the numerical simulations, while the lines represent the b ^ 1 
and 6 <^ 1 analytical asym ptotics, D2 = 1 — 1/nb [Eq. (|3.7|l ] 
and D2 = 26 [Eq. JS^DI]- (|Mirlin and Eversl I2OOOI ). 



the diffusion operator ^Dk^ by \\k\. Consequently, all 
calculations within the weak coupling expansion of the 
cr-model are generalized to the PRBM case by substitut- 
ing n(fc) = i/8|fc| for the diffusion propagator 11 (fc) = 
l/npDk^. The l/|fc| behavior of n(fc) implies that the 
kinetics is superdiffusive, also known as Levy flights, an d 
leads to criticality in ID ( Boucha ud and Geor ged . l l99Cl| ). 



In particular, calculating the average IPR (Pq), one flnds 
the following result for the fractal dimensions 

T, ~ (g- 1)(1 -^VStt^) , g<87r/3/i, (3.7) 

i.e. the weak- multifractality- results (|2.36p . (|2.37p for Tq 
and /(a) with d — 1 and 7 = i/87r/3. For dcfiniteness, 
we focus below on the band center, E = 0, where 7 = 
l/27r/35. 

These results are in good agreement with numerical 
simulations. Figs. [Til EH The deviations from the 
asymptotic (parabolic) form in Fig. 1121 which are par- 
ticularly pronounced at 6 = 1, are a precursor of the 
crossover to the small-6 regime (Sec. IIII.Cp . where the 
parabolic approximation breaks down comple tely. 

The IPR fl uctuations arc also found (Evcrs and Mirliiil . 
I2OOCII : IMirlin and Eversl |2&^ by generalizing the results 
obtained for metaflic samples (iFvodorov and Mirlinl . 
I1995I : IMirlinl . l2000bl: iPrigodin and Altshuled . I1998D . In 
particular, the IPR variance is given for q <^ Q+ib) = 
(2/37r6)i/2 by 



var(P,)/(P, 



(3.8) 



cf . Eq. (|2.50p . Calculating higher cumulants, one can re- 
store the corresponding scale- invariant distribution func- 
tion. 



P(P) = e-^-^exp(-e-^-^) 



where 



P 



(Pn 



27r/36 



(3.9) 
(3.10) 
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FIG. 12 Multifractal spectrum f{a) for 6 = 1 (□) and 
b — 4 (o). Solid line indicates the parabolic approximation 
Eg. (lOTI). Inset: expo nent Tq (□) and t^^p (■) for 6=1. 
(jMirlin and Eversl [2OO0I ') . 




FIG. 13 Evolution of the distribution P(lnP2) for b = 
1 with the system size L (from left to right: L = 
4096, 2048, 1024, 512, (256)). The scale invariance of the IPR 
distribution is clearly seen. (jMirlin and Everd I2OOOI ). 
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FIG. 14 Distribution function P(Pg) at g = 2 (o), 4 (□), and 
6 (o) at 6 = 4 for systems of size L — 4096. The solid line 
represents the analytical result Eq. (|3.9[l . Inset: Asymptotic 
of V{Pa)- Dashed line ind i cates power law with exponent 
Xi = 1.7. (jMirlin and Eversl. l200dh . 



The asymptotic behavior (|3.11|) is outside the reach of 
the numerical simulations for q = 2, however, since the 
condition of its vahdity P ^ 271(3 /q{q — 1) ~ 12.5 corre- 
sponds to very small values of the distribution function 
'P{P) <^ 10~^, and its clear resolution would require a 
much larger statistical ensemble. The situation changes, 
however, with increasing q (see the data for g = 4, 6 
in Fig. [141). Equation p.9p becomes inapphcable (since 
the condition of its validity g <C 9+ is no longer met), 
and the power-law asymptotic behavior p. lip becomes 
clearly seen. In particular, the inset of Fig. [HI shows the 
tail for q = A\ the extracted value of the index Xi ~ 1.7 
is in good agreement with the prediction of the 3> 1 
theory, 2:4 = 7r/2. 



and C ~ 0.5772 is the Euler constant. Equation ([3^ C. Strong multifractality, 6 < 1 

is valid for Pq/{Pq) - 1 < 1. At Pq/{Pq) - 1 - 1 the 
exponential fallofF (|3.9p crosses over to a power-law tail 



ViPq) ^ {Pq/{Pq))- 



with 



Xq = 2TTl3b/q^ 



q^ < 2-K(3b 



(3.11) 



(3.12) 



see the discussion of general properties of the IPR distri- 
bution in Sec. III. C. 51 

The analytical results on the IPR distribution are con- 
firmed by numerical simulations. Figure [TSj demonstrates 
the scale invariance of the distribution. Figure [HI shows 
results for the distribution of the IPRs Pq with q = 2, 
4, and 6 at 6 = 4 (the corresponding value of q+ being 
g+ = (Stt)-^/^ ~ 5). It is seen that at g = 2 the an- 
alytical formula (|3.9p nicely describes the "main body" 
of the distribution, with the upward deviations at large 
P indicating the crossover to the power-law tail (|3.1ip . 



In the q uasi-insulating case b <C 1 the problem can 
be studied (jMirlin and Eversl l2000f) via t he RG method 
earlier developed for related problems bv iLevitovl ( 1990l 
I1999D . The idea of the method is as follows. One starts 
from the diagonal part of the matrix H , each eigenstate 
with an energy Ei — Ha being localized on a single site 
i. Then one includes into consideration non-diagonal ma- 
trix elements Hij with \i—j\ = 1. Most of these matrix el- 
ements are essentially irrelevant, since their typical value 
is ^ b, while the energy difference \Ei — Ej\ is typically 
of order unity. Only with a small probability (~ b) is 
\Ei — Ej\ also of the order of b, so that the matrix ele- 
ment mixes strongly the two states, which are then said 
to be in resonance. In this case one is led to consider a 
two-level problem 



two— level 



E, V 
V E, 



V = H,, 



(3.13) 
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The corresponding eigenfunctions and eigenenergies are 



(+) 



COSf 

sin t 



(-) - 



sin6' 



= {E, + Ej)l2 ± \VW\ + T^ , 



(3.14) 
(3.15) 



where tan 6* = -r + \/l + and t = {Ei - Ej)/2V. 

In the next RG step the matrix elements Hij with \i — 
j\ — 2 are taken into account, then those with — = 3, 
and so forth. Each time a resonance is encountered, the 
Hamihonian is re-expressed in terms of the new states. 
Since the probabihty of a resonance at a distance r is 
~ b/r, the typical scale r2 at which a resonance state 
formed at a scale ri will be again in resonance satisfies 



ln(r2/ri) ~ 1/6 , 



(3.16) 



so that r2 is much larger than ri. Therefore, when con- 
sidering the resonant two- level system at the scale r2 , one 
can treat the ri-resonance state as point-like. 

This leads to the following evolution equation for the 
IPR distribution (for /? = 1): 



d 



2b 

TT 



Tr/2 



de 



[-f{P. 



q, ' 



sin 6 cos^ 9 



XS{Pg 



P^^' COS 



29, 



Pj^^' sm 



29, 



(3.17) 



Eq. (|3.17p is a kind of kinetic equation (in the fictitious 
time t — bhur), with the two terms in the square brackets 
describing the scattering-out and scattering-in processes, 
respectively. 

Multiplying Eq. p.lTp by Pq and then integrating over 
Pq, we get the evolution equation for (Pq) 



d{Pq)/d\nr = -2bT{q){Pq] 



(3.18) 



with 



T{q) 



tt/2 



de 



sin 9 cos2 9 



,c29, 



n29, 



r(2q-l) 



2 r(g-l/2) _ ] 

r(g-l) 229-3 r(g)r(g-l) 



.(3.19) 



Integrating p.l8p from i 
multifractal behavior (P„ 



: 1 to r ^ L, we find the 
with the exponents 



Tq = 2bT{q) 



(3.20) 



It is assumed here that q > 1/2, which is the condition 
of the existence of the integral in Eq. (|3.19p . For q < 1/2 
the resonance approximation breaks down; the exponents 
can be found then from the symmetry relation (j2.38p . 
The function T(q) is shown in Fig. 1151 Its asymptotics 
are 



T{q)^^l/[n{q-l/2)] 
T{q) ^ (2/V^)gi/2 , 



9->l/2 
g > 1 . 



(3.21) 
(3.22) 




2 4 6 8 

q 

FIG. 15 Universal function T(q) characterizing the expo- 
nents r(g) via r(g) = 2bT{q) at fo <C 1. Dashed line in- 
dicates the pole position. Inset: Legendre transform F{A) 
describing the mu l tifrac tal spectrum via f{a) = 2bF{a/2b). 
(|Mirlin and Eversll2000l ). 



We see that the fractal exponents are proportional to 
the small parameter b. This is characteristic of wave 
functions that are very small, typically, with rare and 
strong peaks (resonances) . In the limit 6^0 the fractal 
exponents tend to their insulator value Tq ~ for all 
g> 1/2. 



Legendre transformation of (|3.20p produces the /(a)- 
spectrum of the form 



/(a) - 2bF{A) 



A = a/2b 



(3.23) 



where F{A) is the Legendre transform of T{q). The func- 
tion F{A) is shown in the inset of Fig. [151 its asymptotics 



F{A)c, 
F{A) 



-I/ttA 
A/2, 



A- 
A 



(3.24) 
(3.25) 



Furthermore, it changes sign at A_ ~ 0.5104, corre- 
sponding to qc ~ 2.4056. These analytical findings are 
fully supported by numerical simulations. Fig. 1161 The 
above results, Eqs. are valid for g > 1/2, 

which corresponds to a < 1. The other part of the spec- 
trum can be obtained via the symmetry relation (j2.38p . 
(|2.39p . which is also confirmed by numerical results, see 
Figs. [21 El [21 [Ml 



We return now to the IPR distribution function. Fig- 
ure [T7| shows the results of the numerical integration of 
Eq. (|3.17p for q = 2 with the initial condition /(-P2) = 
6{P2 — 1) a,t t — Q. It is seen that at sufficiently large t 
the distribution of lnP2 acquires a limiting form, shift- 
ing with t without changing its shape. This conclusion of 
scale-invariance of the IPR distribution is also supported 
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FIG. 16 Multifractal spectrum f{a) for b = 0.25 (O) and 
6 = 0.1 (A). Inset: exponent r(g). Dashed and dotted 
lines indicate the analy tical results Eqs. (|3.23p and (|3.20|) . 
(jMirlin and Eversll2000l ) . 




FIG. 17 Flow of the distribution of lnP2 calculated from the 
kinetic equation (|3.17p at t — blur = 1.2. ..1.7 (from right 
to left). The oscillations near l nP2 ~ —1.5 are numer ical 
artifacts due to rounding errors. IIMirlin and Eversl . l2000D . 



by analytical arguments: it is found that the fixed-point 
distribution has the form 



(3.26) 



with r*^P as defined in Sec. III. C. 51 In agreement with a 
general discussion in Sec. lII.(J?5] the distribution is found 
to have a power-law tail, /o(-fg) ~ p-^q-i ^ with the 
index Xq given by Eqs. (I^Tii)) . (PTiS)) . 

All the formulas of this subsection remain valid for the 
case /? = 2, with a replacement 6 — > (7r/2\/2)6. 



D. Levels statistics 



the form^ (iKravtsov and Muttalibl . Il997t iMirlhl l2000bl: 
iMirlin and Eversll2000D 



R'^\s) = 5{s) 



sin^(7rs) {irs/AhY 
(7rs)2 sinh^(7rs/46) 



(3.27) 



The correlation function (|3.27p follows the RMT result 

rI^\s) = 5{s) — sin^(7rs)/(7rs)2 up to the scale s ^ b 
(playing the role of the Thouless energy here), and then 
begins to decay exponentially. The spectral compress- 
ibility at 6 3> 1 is given by 



X ^ 1/2^/36 



6> 1 



(3.28) 



In the opposite limit, 5^1, the evolution equation for 
R^juj^r) can be written down in analogy with Eq. ()3.17p 
(jMirlin and Eversl 120001 ). with the result (for (3=1) 



4 («) = S[s) - erfc(|s|/2V^6) , (3.29) 

where erfc(a;) = (2/y^) J^°° exp(— i^)^^ is the error func- 
tion. This yields the spectral compressibility 



X 



1 - 46 



5< 1 



(3.30) 



The results for the (3 — 2 case are qualitatively similar, 
(s) = (5(s) - exp(-s V27r6^) (3.31) 



X 



1 -7rV2fe 



6< 1 



(3.32) 



Thus, in the limit of small b the level repulsion is efh- 
cient in a narrow region |s| < b only, and the spectral 
compressibility tends to the Poisson value x = 1- The 
physical reason for the reduced range of the level repul- 
sion is quite transparent. Consider two nearby in energy 
states separated by a typical distance r ^ L \a. the coor- 
dinate space. If their energy difference s b, such two 
states will form a resonance pair, so that their levels will 
repel. On the other hand, if s ^ &, these two states will 
not be in resonance, their wave functions remain weakly 
overlapping, and the level repulsion between them will 
be inefRcient. 

These results are fully supported by the numerical 
data. Figure [15] represents the level correlation function 
i?2(s) at 6 = 0.1, in agreement with Eq. ([X^ . The level 
number variance var[n(5)] is plotted versus the average 
{n{£)) in Fig. [HI the data show an extended plateau 
region in var[n(£)]/(n(£)), determining The upper 
bound for this region is set by the matrix size L, while 
the lower bound is ~ b. The numerically obtained spec- 
tral compressibility in a broad range of b is shown in 
Fig. [201 in the large-6 and small-6 regions it agrees well 
with the corresponding analytical asymptotics. 



In the 6^1 regime the two-level correlation func- 
tion (|2.70p is obtained by an appropriate generalization 
of the earlier findings for the diffusive samples. In par- 
ticular, considering for simplicity the [3 = 2 ensemble 
at the band center, the level correlation function has 



The precise form of the level correlation function R2 depends 
on the boundary conditions, which are chosen to be periodic in 
this subsection. The value of the spectral compressibility r) is 
independent on the boundary conditions. 
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FIG. 18 Two-level correlation function R2{s) for two system 
sizes L = 256 (o) and L = 512 (□) at b = 0.1. The solid 
line i ndicates the theoretical result (|3.29p . (jMirlin and Everi . 
|2000| 1. 



0.8 



0.6 



!§ 0.4 



0.2 



o o 



fl i i ■ ■ 



fl □ 



• • • • . 8 



10" 



10 



<n) 



10- 



10^ 



FIG. 19 Variance of the number of levels {5n^) as a function 
of the energy width of the interval parametrized by the mean 
level number (n). Traces correspond to fe = 1 (open o; L = 
4096, fiUed: L = 2048), b = 0.25 (open □: L = 4096, filled: 
L = 2048) and b = 0.05 (open O: L = 4096, O with dot: 
L = 1024, filled: L = 512). Th e dashed line ind i cates the 
analytical prediction Eq. (fTSUI) . (|Mirlin and Eversl . [2OO0I) . 
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E. Boundary criticality 

The presentati o n in this subsection follows 
iMildenberger et al\ (|2007bl) . In spirit of Sec. III.C.81 
we consider the critical PRBM model model with a 
boundary at i = 0, which means that the matrix element 
Hij is zero whenever one of the indices is negative. 
The implementation of the boundary is, however, not 
unique. The important point is that this degree of 
freedom affects the boundary criticality. Specifically, 
one should specify what happens with a particle which 
"attempts to hop" from a site i > to a site j < 0, 
which is outside of the Hilbert space. One possibility 
is that such hops are simply discarded, so that the 
matrix element variance P) = Jij is simply given by 
J,j = [1 + (i _ for i,j > 0. More generally, the 

particle may be reflected by the boundary with certain 
probability p and "land" on the site — j > 0. This leads 
to the following formulation of the model. 



J, 



1 



+ 



P 



(3.33) 



While the above probability interpretation restricts p to 
the interval [0,1], the model is defined for all p in the 
range —1 < p < 00. The newly introduced parameter p 
is immaterial in the bulk, where i,j 3> |i — i| and the 
second term in Eq. ()3.33p can be neglected. Therefore, 
the bulk exponents depend on b only (and not on p), 
and their analysis in Sec. 1111.51 IIII.CI remains applicable 
without changes. On the other hand, as discussed below, 
the surface exponents are function of two parameters, 
b and p. 

Equation (|3.33p describes a semi-infinite system with 
one boundary at i = 0. For a finite system of a length 
L (implying that the relevant coordinates are restricted 
to < i,j < L) another boundary term, p'/[l + {i + 
j — 2L)^/&^], is to be included on the right-hand side of 
Eq. (j3.33|) . In general, the parameter p' of this term may 
be different from p. This term, however, does not affect 
the boundary criticality at the i = boundary, so it is 
discarded below. 

The regime of weak criticality, 6 3> 1, can be studied 
via a mapping onto the supermatrix u-model as in the 
bulk case. Sec. IIII.BI This results again in an approxi- 
mately parabolic spectrum, which, however, differs by a 
constant factor Rp from its bulk counterpart. 



A' 



[<Z(1 - q)/27Tpb]Rp 



(3.34) 



This factor is determined from the solution of the classi- 
cal Levy flight problem with boundary, 



FIG. 20 Spectral compressibility x a function of b: 
crossover from the "quasi-metallic" (6 3> 1) to the "quasi- 
insulating" (6 ^ 1) behavior. The lines indicate the analyt- 
ical results for b > 1 a nd 6 < 1, Eqs. ((3?28ll and ((3^ . 
(jMirlin and Eversll2000l ) . 



dt 



6ij Jg 



Wj{t)=0, (3.35) 



where J, 



(i) 



Tli'kLoJik, with the initial condition 
Wi(0) = Sir, where r is near the boundary. Specifically, 
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FIG. 21 The ratio Rp = Al{b,p)/ Ag(b) of the surface and 
bulk anomalous exponents for large 6, as a function of the 
reflection parameter p. Diamonds represent the results of 
the (T-model analysis with a numerical solution of the cor- 
responding classical evolution equation (|3.35p . Circles rep- 
resent a direct computer simulation of the PRBM model 
with & = 8. The ratio Rp has been evaluated for the range 
< g < 1, where the numerical accuracy of the anomalous 
exponents is the best. Within this interval the obtained Rp 
is q-indepen dent (within numerical erro rs) in agreement with 
Eq. (ITM)) . l|Mildenberger et am2007bl ). 



the return probability Wr{t) decays with time as 1/t, and 
the corresponding prcfactor yields the fractal exponents, 

Ayq{l-q)^ r'tWr{t)\t^co- (3.36) 

The evolution equation p.35p was solved numerically; 
the obtained results for Rp are shown in Fig. [21] For p = 
1 the evolution equation can be obtained from its bulk 
counterpart (defined on the whole axis, — oo < i < oo) 
by "folding the system" on the semiaxis i > according 
to Wi{t) + W-i{t) — > Wi{t), leading to the analytical 
result 

Ri = 2. (3.37) 

As in the bulk case, the regime of small b can be studied 
via the real-space RG method (Sec. IIII.Cp . The multi- 
fractal exponents are found to be 

T^^ = {l+pf/%T{q) = [(1 +p)i/V2]r^ (3.38) 

with T{q) given by Eq. (|3.19p . i.e. they differ from the 
bulk exponents by the factor (1 -|-p)^/^. In full analogy 
with the bulk formula p.20p . the result (|3.38p is valid 
for q > 1/2, where the multifractal exponent Tq is small. 
The results can, however be extended to the range of 
q < 1/2 by using the symmetry relation (|2.38p . For all 
q the obtained relation between the surface and the bulk 
multifractal spectra can be formulated in the following 
way: 

r^^(6,p) = r,^6^ 6(1+^)^2/2). (3.39) 

These predictions were corroborated by numerical sim- 
ulations of the PRBM model. Figure [221 shows the de- 
pendence of the anomalous dimension A2 = D2 — 1 on 6 
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FIG. 22 Upper panel: Anomalous exponent A2 = D2 — 1 as 
a function of b from numerical simulations in the bulk and at 
the boundary for the reflection parameter p = and 1. The 
inset shows data for p = 3 compared to the p = bulk values. 
Lower panel: Surface and bulk data for p = Q compared with 
analytical results for small and large b (using Rg = 2.78) , 
Eqs. ^J^, Jajll), ([3:2011 . ([338|l . ([Mildenberger et aLl . [2007bl ). 



in the bulk and at the boundary, for three different values 
of the reflection parameter p. In Fig. [221 the whole mul- 
tifractal spectra are shown for fixed large values of 
6. Specifically, the anomalous dimensions and A^ are 
presented for 6 = 2, 4 and 8, with the reflection parame- 
ter chosen to be p = 1 . For all curves the q dependence is 
approximately parabolic, as predicted by the large-6 the- 
ory, Eq. (|3.34p , with the prefactor inversely proportional 
to h. It is also seen in Fig. [23 that the bulk multifractal- 
ity spectrum for 6 = 4 and the surface spectrum for 6 = 8 
are almost identical, in agreement with Eq. (|3.37p . The 
same is true for the relation between the bulk spectrum 
for 6=2 and the surface spectrum for 6 = 4. The ra- 
tio of the large-6 surface and bulk anomalous dimensions, 
Rp — A^/A^, is in good agreement with the cr-model pre- 
dictions for Rp, as shown in Fig.[2U In Fig.[211the surface 
and bulk multifractal spectra are shown for the case of 
small 6. While the spectra are strongly non-parabolic in 
this limit, they clearly exhibit the symmetry q —t 1 — q, 
Eq. (|2.38p . The data are in good agreement with the 
real-space RG results. 



F. Further related activities 

As a model of the Anderson critical point, the 
PRBM ensemble has attracted a lot of interest 
during recent years. In view of space limita- 
tions, we can only briefly mention some of these 
works. The corresponding research directions in- 
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FIG. 23 Upper panel: Boundary and bulk multifractal spec- 
tra, Aq and A^, at 6 = 2, 4, and 8 for the reflection parameter 
p = 1. In accordance with Eq. (|3.37p . the surface multifrac- 
tality spectrum is enhanced by a factor close to two compared 
to the bulk. Middle panel: Surface spectrum divided by the 
analytical 6^1 result, Eq. (I3.34|l . With increasing 6, the nu- 
merical data converges towards the analytical result (dashed 
line). Lower panel: Analogous plot for the bulk spectrum. 
(jMildenberger et ali^mnh . 




q 

FIG. 24 Main panel: Numerically determined boundary and 
bulk anomalous dimensions Ag at 6 = 0.1 for p = 0, 1, and 3. 
As expected, the bulk anomalous dimension is independent 
of the value of p. In accordance with Eq. (|3.39|l . for p = 3 
surface and bulk dimensions have the same values. Inset: The 
p = data compared to the analytical results, surface [solid 
line, Eq. ((331}] and bulk [dashed Hue, Eq. (fT20) l]. Analytical 
data have been calculated for q > 0.6 and mirrored for q < 0.4 
by using the symmetry relation A, — Ai-,. In the vicinity of 
q = 1/2, at \q - 1/2| < l/lnb~\ the ana lytical result 
breaks down. (jMildenberger et aLll2007bl ). 



elude: relation to the sy stem of interacting: fermion s 
in ID (Luttinger liquid) (|Kravtsov and Tsvelikl |2000D : 
a generalization of the mod el to two dimensions 
( Potempa and Schweitzeii |2002|): wave f unctio n statis- 



pr oof of the comp l etene ss of the classification was given 



tics ( Cuevasl . ~ 2003allbi r Cuevas et al. 2002t Varga . 



|S|;l^corS^s hZ:^^}:^^^ 
20061 ): virial expansion for level statistics of al 



most diagonal random ma. t rices ( Kravtsov et al\ . 120061 : 
lYevtushenko and Kravtsovl [200^ 12004D; chiral PRBM 
and possible applications to qua i itum chromo dynam- 
ics ( Garcia-Garcia and Takahashil . |2004[ ). and manifes- 
tations of multifr actality in scattering characteri s tics of 
an open system (iMendez-Bermudez and Kottos'. '2005'; 
iMendez-Bermudez and Varg al l2006l ). Rotationally- 
invariant random matrix ensembles with level statistics 
largely similar to t hat in the critical PRBM ensem- 
bles were studied in lMoshe et al\ (|l994f) : iMuttalib et al\ 
(|200lL 119931 ): a relation between these ensembles and the 
PRBM model is di scussed in ( Kravtsov and Muttalibl . 
119971: lMirlinl . l2000bl) . 



IV. SYMMETRIES OF DISORDERED SYSTEMS 

In this section we briefly review the symmetry classi- 
fication of disordered systems based on the relation to 
the classical symmetric space s, which was estab lished in 
lAltland and Zirnbaued (|l997l ): IZirnbaueil (| 19961 ). For a 
detailed presentation of the scheme and the underlying 
mathemat ical structures the reader is referred to a recent 
review bv lCaselle and Magneal (|2004l) . A mathematical 



r oot ot the comp l etene s 
i lHeinzner eiM (|2005[ ) 



A. Wigner- Dyson classes 

The ran dom matr i x theo ry (RMT) was introduced into 
physics by Wigneil ( 195ll ). Developing Wigner's ideas, 
iDvsonI (|l96a) " put forward a classification scheme of en- 



sembles of random Hamiltonians. This scheme takes into 
account the invariance of the system under time rever- 
sal and spin rotations, yielding three symmetry classes: 
unitary, orthogonal and symplectic. 

If the time-reversal invariance (T) is broken, the 
Hamiltonians are just arbitrary Hermitian matrices. 



H = H^ 



(4.1) 



with no further constraints. This set of matrices is in- 
variant with respect to rotations by unitary matrices; 
hence the name "unitary ensemble" . In this situation, 
the presence or absence of spin rotation invariance {S) 
is not essential: if the spin is conserved, H is simply 
a spinless unitary-symmetry Hamiltonian times the unit 
matrix in the spin space. In the RMT one considers most 
frequently an ensemble of matrices with independent, 
Gaussian-distributed random entries, which is called the 
Gaussian unitary ensemble (GUE). While disordered sys- 
tems have, of course, much richer physics than the Gaus- 
sian ensembles, their symmetry classification is inherited 
from the RMT. 
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Let us now turn to the systems with preserved time- 
reversal invariance. The latter is represented by an an- 
tiunitary operator, T = KC, where C is the operator of 
complex conjugation and K is unitary. The time-reversal 
invariance thus implies H — KH^^K^^ (we used the Her- 
miticity, H* = H"^). Since acting twice with T should 
leave the physics unchanged, one infers that K*K = p, 
where p = ±1. As was shown by Wigner, the two cases 
correspond to systems with integer (p ~ +1) and half- 
integer [p = —1) angular momentum, li p — 1, a repre- 
sentation can be chosen where K — 1, so that 



(4.2) 



The set of Hamiltonians thus spans the space of real sym- 
metric matrices in this case. This is the orthogonal sym- 
metry class; its representative is the Gaussian orthogonal 
ensemble (GOE). For disordered electronic systems this 
class is realized when spin is conserved, as the Hamil- 
tonian then reduces to that for spinless particles (times 
unit matrix in the spin space). 

If T is preserved but S is broken, we have p = —1. In 
the standard representation, K is then realized by the 
second Pauli matrix, K — iuy, so that the Hamiltonian 
satisfies 



H 



(4.3) 



It is convenient to split the 2N x 2N Hamiltonian in 2 x 2 
blocks (quaternions) in spin space. Each of them then is 
of the form q — qoao + iqiCTx + iq2'Jy + iqaaz (where ctq is 
the unit matrix and cr^.y^z the Pauli matrices), with real 
q^, which defines a real quaternion. This set of Hamil- 
tonians is invariant with respect to the group of unitary 
transformations conserving ay, UdyU'^ — (Jy, which is 
the symplectic group Sp(2iV). The corresponding sym- 
metry class is thus called symplectic, and its RMT rep- 
resentative is the Gaussian symplectic ensemble (GSE). 



B. Relation to symmetric spaces 

Before discussing the relation to the families of sym- 
metric spaces, we brie fly remind the reader how the latter 
are c onstructed ( Caselle and Magneal . |2004 iHelgasoil 
Il978l ). Let G be one of the compact Lie groups SU(A^), 
S0(7V), Sp(2iV), and g the corresponding Lie algebra. 
Further, let 9 be an involutive automorphism — > 
such that 9^ = 1 but 9 is not identically equal to unity. 
It is clear that 9 splits g in two complementary sub- 
spaces, = J? e Cp, such that 9{X) ^ X ior X e ^ 
and 9{X) = —X for X € *p. It is easy to see that the 
following Lie algebra multiplication relations holds: 

[Si,Si]cK, [V,V]cSi. (4.4) 

This implies, in particular, that ^ is a subalgebra, 
whereas *P is not. The coset space G/K (where K is the 
Lie group corresponding to A) is then a compact symmet- 
ric space. The tangent space to G/K is Cp. One can also 



construct an associated non-compact space. For this pur- 
pose, one first defines the Lie algebra g* = R(Bi^, which 
differs from g in that the elements in *p are multiplied 
by i. Going to the corresponding group and dividing K 
out, one gets a non-compact symmetric space G* / K . 

The groups G themselves are also symmetric spaces 
and can be viewed as coset spaces G x G/G. The corre- 
sponding non-compact space is G^/G, where G^ is the 
complexification of G (which is obtained by taking the 
Lie algebra g, promoting it to the algebra over the field 
of complex numbers, and then exponentiating). 

The connection with sy mmetric spaces is now estab- 
lished in the following way ( Altland and Zirnbaueil[l997l : 
Zirnbauer, 1996). Consider first the unitary symmetry 
class. Multiplying a Hamiltonian matrix by i, we get 
an antihermitean matrix X = iH. Such matrices form 
the Lie algebra u{N). Exponentiating it, one gets the 
Lie group U(A'^), which is the compact symmetric space 
of class A in Cartan's classification. For the orthogonal 
class, X ~ iH is purely imaginary and symmetric. The 
set of such matrices is a linear complement 'P of the alge- 
bra ^ = (N) of imaginary antisymmetric matrices in the 
algebra g — u{N) of antihermitean matrices. The corre- 
sponding symmetric space is G/K = U(iV)/0(iV), which 
is termed AI in Cartan's classification. For the symplec- 
tic ensemble the same consideration leads to the symmet- 
ric space U(A^)/Sp(A^), which is the compact space of the 
class AIL If we don't multiply H hy i but instead proceed 
with H in the analogous way, we end up with associated 
non-compact spaces G*/K. To summarize, the linear 
space *P of Hamiltonians can be considered as a tangent 
space to the compact G/K and non-compact G* /K sym- 
metric spaces of the appropriate symmetry class. 

This correspondence is summarized in Table HI where 
the first three rows correspond to the Wigner-Dyson 
classes, the next three to the chiral classes (Sec. IIV.CP 
and last four to the Bogoliubov-de Gennes classes 
(Sec. lIV.Pp . The last two columns of the table specify 
the symmetry of the corresponding cr-model. In the su- 
persymmetric formulation, the base of the tr-model tar- 
get space Mb x Mj^ is the product of a non-compact 
symmetric space Mb corresponding to the bosonic sec- 
tor and a compact ( "fermionic" ) symmetric space My^. 
(In the replica formulation, the space is Mb for bosonic 
replicas or My^ for fermionic replicas, supplemented with 
the replica limit n ^ 0.) The Cartan symbols for these 
symmetric spaces are given in the sixth column, and the 
compact components Myr are listed in the last column. 
It should be stressed that the symmetry classes of A^g 
and Myr are different from the symmetry class of the en- 
semble (i.e. of the Hamiltonian) and in most cases are 
also different from each other. Following the common 
convention, when we refer to a system as belonging to 
a particular class, we mean the symmetry class of the 
Hamiltonian. 

It is also worth emphasizing that the orthogonal groups 
appearing in the expressions for Myr are 0{N) rather 
than SO{N). This difference (which was irrelevant when 
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TABLE I Symmetry classification of disordered systems. First column: symbol for the symmetry class of the Hamiltonian. 
Second column: names of the corresponding RMT. Third column: presence (+) or absence ( — ) of the time-reversal (T) and 
spin-rotation (S) invariance. Fourth and fifth columns: families of the compact and non-compact symmetric spaces of the 
corresponding symmetry class. The Hamiltonians span the tangent space to these symmetric spaces. Sixth column: symmetry 
class of the cr-model; the first symbol corresponds to the non-compact ( "bosonic" ) and the second to the compact ( "fermionic" ) 
sector of the base of the cr-model manifold. The compact component M f (which is particularly important for theories with 
non-trivial topological properties) is explicitly given in the last column. 



we were discussing the symmetry of the Hamiltonians, as 
it does not affect the tangent space) is important here, 
since it influences topological properties of the manifold. 
As we will detail in Sec. IVIVIi the topology of the cr- 
model target space often affects the localization proper- 
ties of the theory in a crucial way. 



C. Chiral classes 

The Wigner-Dyson classes are the only allowed if one 
looks for a symmetry that is translationally invariant in 
energy, i.e. is not spoiled by adding a constant to the 
Hamiltonian. However, additional discrete symmetries 
may arise at some particular value of energy (which can 
be chosen to be zero without loss of generality), leading to 
novel symmetry classes. As the vicinity of a special point 
in the energy space governs the physics in many cases (i.e. 
the band center in lattice models at half filling, or zero 
energy in non-s-wave superconductors), these ensembles 
are of large interest. They can be subdivided into two 
groups - chiral and Bogoliubov - de Gennes ensembles - 
considered here and in Sec. lIV.Di respectively. 

The chiral ensembles appeared in both contexts of par- 
ticle physics and physics of dis o rdere d electronic syst ems 
abou t fifteen years ago (iGadel 1199,4 iGade and Wegne^ . 



1991 



ISlevin and Nagad 19931 : IVerbaarschot and Zahed 



19931 ). The corresponding Hamiltonians have the form 



H = 




(4.5) 



i.e. they possess the symmetry 



t,Ht, = -H 



(4.6) 



where is the third Pauli matrix in a certain "isospin" 
space. In the condensed matter context, such ensembles 
arise, in particular, when one considers a tight-binding 
model on a bipartite lattice with randomness in hopping 
matrix elements only. In this case, H has the block struc- 
ture (|4.5p in the sublattice space. 

In addition to the chiral symmetry, a system may pos- 
sess time reversal symmetry and/or spin- rotation invari- 
ance. In full analogy with the Wigner-Dyson classes, 
IIV.AI one gets therefore three chiral classes (unitary, or- 
thogonal, and symplectic) . The corresponding symmetric 
spaces, the Cartan notations for symmetry classes, and 
the cr-model manifolds are given in the rows 4-6 of the 
Table H 



D. Bogoliubov - de Gennes classes 



As was found in lAltland and Zirnbaued (|l997l ). the 
Wigner-Dyson and chiral classes do not exhaust all pos- 
sible symmetries of disordered electronic systems. The 
remaining four classes arise most naturally in supercon- 
ducting systems. The quasiparticle dynamics in such 
systems can be described by the Bogoliubov-de Gennes 
Hamiltonian of the form 

N ^ N 

H = Y^ Kpcicp + 2 ^ - A^^CaC/s) (4.7) 

a/3 a/3 
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where and c are fermionic creation and annihilation 
operators, and the N x N matrices h, A satisfy h — 
and A"^ — — A, in view of hermiticity. Combining cj,, 
in a spinor -0^ = (c^ ,Ca), one gets a matrix representa- 
tion of the Hamiltonian, H — 'ip^ Hip, where H has the 
structure 



H = 




h^h'^ , A = -A^ . (4.. 



The minus signs in the definition of H result form the 
fermionic commutation relations between and c. The 
Hamiltonian structure (|4.8p corresponds to the condition 



(4.9) 



H = -T^H^ 



(in addition to the Hermiticity H — H^), where Tx is the 
Pauli matrix in the particle-hole space. Alternatively, 
one can perform a unitary rotation of the basis, defining 
H = g'' Hg with g = {l+iTx)/^/2. In this basis, the defin- 
ing condition of class D becomes H = —H^, so that H 
is pure imaginary. The matrices X = iH thus form the 
Lie algebra so (27V), corresponding to the Cartan class D. 
This symmetry class described disordered superconduct- 
ing systems in the absence of other symmetries. 

Again, the symmetry class will be changed if the time 
reversal and/or spin rotation invariance are present. The 
difference with respect to the Wigner-Dyson and chiral 
classes is that now one gets four different classes rather 
than three. This is because the spin-rotation invariance 
has an impact even in the absence of time-reversal invari- 
ance, since it combines with the particle-hole symmetry 
in a non-trivial way. Indeed, if the spin is conserved, the 
Hamiltonian has the form 



E. Additional comments 

i) In addition to Table [H where the symmetry classes 
are ordered according to the discrete symmetries of the 
underlying physical systems, we include Table |TT] with 
a more mathematical ordering. There, the first row is 
formed by compact symmetric spaces that are groups and 
the rest is a "matrix" of symmetric spaces G/K arranged 
according to the type (U, Sp, or O) of the groups G and 
K. This ordering illustrates the completeness of the clas- 
sification scheme: all entries in the "matrix" are filled, ex- 
cept for two, as there is no symmetric spaces of the 0/Sp 
and Sp/0 types. Also, we will see in Sec. |VT]that this 
ordering is relevant to types of 2D critical behavior that 
the corresponding systems may show. In particular, the 
first row are the classes where different types of the QHE 
(IQHE, SQHE, and TQHE) take place. The second row 
is formed by the classes allowing for the Wess-Zumino 
type of criticality, while the third row are the systems 
which allow for a Z2 topological term. The diagonal of 
the "matrix" is occupied by three chiral classes. 

ii) Strictly speaking, one should distinguish between 
the orthogonal group SO(iV) with even and odd N, which 
form different Cartan classes: S0(2A^) belongs to class D, 
while S0(2iV-|-l) to class B. In the conventional situation 
of a disordered superconductor, the matrix size is even 
due to the particle-hole space doubling, see Sec. IIV.DI 
It was found, however^ that the class B can arise in p- 
wave vortices ( Ivanovl . 2002al lbh. In the same sense, the 
class Dili should be split in Dlll-even and Dlll-odd; the 
last one represented by the symmetric space S0(4Af -|- 
2)/U(2iV -I- 1) can appear in vortices in the presence of 
time-reversal symmetry. 



N N 

(4.10) 

where h and A are NxN matrices satisfying h = W and 
A = A"^. Similar to (l48l) . we can introduce the spinors 
V'l = (c]^|5Ci|) and obtain the following matrix form of 
the Hamiltonian 



H = 



A 



h = h) , A = A^ 



It exhibits a symmetry property 



H = -TyH' Ty. 



(4.11) 



(4.12) 



The matrices H = iX now form the Lie algebra sp(2iV), 
which is the symmetry class C. 

If the time reversal invariance is present, one gets two 
more classes (CI and Dili). The symmetric spaces for 
the Hamiltonians and the u-models corresponding to the 
Bogoliubov-de Gennes classes are given in the last four 
rows of the Table HI 



F. Perturbative /^-functions for cr-models of different 
symmetry classes 

Perturbative /^-functions for cr-models on all t he types 
of symmetric space s were in fact calculated ( Hikamil . 
I198II : IWegneij . Il989l ) long before the physical significance 
of the chiral and Bogoliubov-de Gennes classes has been 
fully appreciated. In Table Hill we have collected the re- 
sults for all /3-functions, /3{t) — —dt/dlnL, in d = 2 -I- e 
dimensions up to four-loop order. Here t is the coupling 
constant inverse proportional to the dimensionless con- 
ductance g, and the Anderson localization problem cor- 
responds to the replica limit N = p = 0. The corre- 
sponding results for the Wigner-Dyson classes have al- 
ready been quoted in Sec. III. B. 21 For each symmetry 
class of disordered electronic systems, the perturbative 
/3-function can be equivalently obtained from either com- 
pact or non-compact cr-models on the appropriate man- 
ifolds. As an example, the /^-function for the Wigner- 
Dyson orthogonal class can be found by using the replica 
limit of the compact cr-model of the type CII or of the 
non-compact cr-model of the type BDI. 

It is seen that for the classes A, AI, C, CI the f3- 
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U(A^) A 


Sp(2A') C 


SO(A') BD 


U(p + g)/U(p)xU(g) AIII 


Sp(2iV)/U(iV) CI 


SO(2Af)/U(iV) Dili 


U(2iV)/Sp(2iV) All 


Sp(2p + 2g)/Sp(2p)xSp(2g) CII 




U(iV)/0(iV) AI 




SO(p + g)/SO(p)xSO(g) BDI 



TABLE II Compact symmetric spaces arranged in a form of a "matrix", with the corresponding Cartan symbols. First row: 
U, Sp, and O groups. Second (third, fourth) row: symmetric spaces G/K (which are not groups themselves) with K being the 
unitary (resp. symplectic, orthogonal) group. 



function is negative in 2D in the replica limit N = p = 
(at least, for small t). This indicates that normally all 
states are localized in such systems in 2D. (This con- 
clusion can in fact be changed in the presence of topo- 
logical or Wess-Zumino terms, Sec. IVLAI ) Above 2D, 
these systems undergo the Anderson transition that can 
be studied within the 2 + e expansion, Sec. III. B. 21 For 
the classes AIII, BDI, and CII (chiral unitary, orthogo- 
nal and symplectic classes, respectively) the /3-function 
is exactly zero in 2D, implying a line of fixed points, 
Sec. IVLFl Finally, in the classes All, D, and Dili the 
/3-function is positive at small t, implying the existence 
of a metal- insulator transition at strong coupling in 2D. 
The 2D Anderson transitions in the Wigner-Dyson sym- 
plectic class All and in the Bogoliubov-de Gennes class 
D will be discussed in more detail in Sec . I VI. B1 and I VI . El 
respectively. 



V. QUASI-ID SYSTEMS: DISORDERED WIRES 

Under usual conditions, electronic states in disordered 
wires are localized, with localization length ^ ~ Nl, 
where / is the mea n free path and N the number of con- 
ducti ng channels ( Berezinskvl Il974l : lEfetov and Larkiiil . 
I1983D . This is, however, not the full story, and that is 
why we include a section about quasi-lD systems in this 
review devoted to localization-delocalization transitions 
and criticality. In fact, one route to delocalization and 
criticality in systems of one ID geometry - long-range 
l/r hopping - has already been discussed in Sec. IIIII In 
this section we consider possible types of delocalization 
in disordered wires which are related to the symmetries 
of the problem. We will see later that in many cases 
there are close connections between such phenomena in 
quasi-lD and 2D systems fSec. IVip. 

Two approaches to quasi- ID disordered electronic sys- 
tems have been developed. The first one is the super- 
symmetric (T-model approach. Sec. III. B. II For the wire 
geometry, the cr-model field Q depends on the longitudi- 
nal coordinate only. As a result, the problem becomes 
a kind of imaginary-time quantum mechanics on the a- 
model manifold, with the longitudinal coordinate playing 
the role of the (imaginary) time. This has allowed re- 
searchers to get a number of exact results for systems of 
Wigner-Dyson symmetry classes, including the asymp- 
totics of the density-density correl ation function and 
the value of the localization length ( Efetov and Larkinl . 



Il983l), a detailed description of the wave function statis- 
tics ()Fvodorov a nd MirliUj 199 4), average conductanc e 
and its variance (iMirlin aZ.I . [1994 IZirnbaueii Il992[ ). 
For reviews of results for statistical properties of vari- 
ous quantities in disordered wires obtaine d within the c r- 
model approa ch the reader is referred to lEfetovl ( 19971 ): 
iMirlinI (|2000bD . 

The second approach is based on the description of 
a wire in terms of its transfer-matrix M. The evolu- 
tion of the corresponding distribution V{M) with the 
wire length is described by the Dorokhov-Mello-Pereyra- 
Kumar (DMPK) equation (' Dorokhovill982HMello et all 
.1988) . At variance with the cr-model approach, which 
allows to address any observable, the DMPK approach 
is designed to study the transport properties. On the 
other hand, the advantage of the DMPK approach is that 
it allows one to study wires with arbitrary number N 
of channels (not necessarily iV ^ 1 as required for the 
derivation of the diffusive cr-model) and provides a very 
detailed information on the whole distribution of trans- 
mission eigenvalues. The two approaches are thus com- 
plementary; their equivalence for transport properties of 
many-channel wires (wh en both of them are applicable) 
was explicitly shown in iBrouwer and Frahml U996). A 
review of the results of the DMPK method for Wigner - 
Dyson symmetry classes was given in iBeenakkeil (|l997l ). 

Peculiar transport properties of disordered wires of un- 
conventional symmetry classes have been mainly studied 
within the DMPK approach. Below we briefly describe 
this method and present the key results. Whenever ap- 
propriate, we will also make contact to the results ob- 
tained within the cr-model formalism. 



A. Transfer matrix and DMPK equations 

In the transfer matrix approach, one imagines the wire 
attached to two clean electrodes, where one can define 
asymptotic states. This allows to formulate a scattering 
problem. The transfer matrix M relates the amplitudes 
in N incoming and N outgoing channels to the right of 
the wire to the corresponding amplitudes on the left side 
of the wire. 




(5.1) 
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TABLE III Perturbative /3-functions, f3{t) = —dt/d\n L, up to the four- lo op order for the cr-models in d = 2 + e dimensions 
with symmetric spaces as target manifolds (Hikam|, llOSlI : IWegned . Il989l ). First column: Cartan symbol for the cr-model 
symmetric space. Second column: compact tj-model manifold. (The associated non-compact spaces can be found in Table |l]) 
Third column: symmetry class of random Hamltonians described by the replica version of the compact cr-model (c) and of 
its non-compact counterpart (n-c). The last five columns give the coefficients of the terms from t to f in the /3-functions 
for compact cr-models. The /3-functions for the corresponding non-compact cr-models are obtained by the substitution j3{t) 
—(5{~t), i.e. by flipping the sign of the terms with even powers of t. Following notations are used: a = ^(^(3), C5{N,p) = 
-{§N\N -p) + 5p\N - pf - ^iV^ - (f + 8a)pN{N - p) + (-§ + 16a)p{N - p) + H + 16a)N^ + (i - 64a)Ar + 64a); 
C4(iV,p) = ^(lpN\N^)^p\N^)^^N\llp{N^}iGy, c3(iV) = -{f^N^ + mN' + mN^ + ^N+Sl^)- c^iN) = -(f iV^- 
(f - 8a)Af^ + (9 + 8a)Af - 1) ; ci{N) = -~[{§+a){N ~2f^{N^3){N^){N+2)] 



(In fact, for symmetry classes A, C, and BD with broken 
time-reversal symmetry the number of incoming and out- 
going channels can differ. We will discuss this peculiar 
situation in Sec. IV. Dl ) The requirement of the current 
conservation 



|^out|2_ |^i„|2 ^ 



(5.2) 



implies that M is an element of the pseudounitary group 
G = \J{N,N). The transfer matrix can be presented in 
the form (Cartan decomposition) 



M = 




cosh X sinh x 
sinhi; coshx 




(5.3) 



where x = diag(xi, . . . , xat) are "radial coordinates", 
while the left and right matrices ( "angular coordinates" ) 
are the elements of if = l]{N)xlJ{N). If the Hamil- 
tonian of the system possesses some additional (time- 
reversal, spin-rotation, chiral, particle-hole) symmetries, 
the group G of the transfer matrices will change corre- 
spondingly, with K being the maximal compact subgroup 
of G. The coset spaces G/K (which play the central role 
for the DMPK equations, see below) are non-compact 
symmetric spaces; they are listed for all the symmetry 
classes of the Hamiltonian in Table IIVI The dinien- 
sionless conductance of the wire is expressed in terms 
of the radial coordinates a;„ as G = d^^^^Tn, where 
T!„ = 1/ cosh^ Xn are transmission eigenvalues and d is 
the degeneracy (1, 2, or 4) depending on the symmetry 
class. 



When an additional thin slice with a transfer matrix 
T is attached to the wire, the new transfer matrix is ob- 
tained by simple multiplication M' = TM . For a given 
distribution function of the elementary transfer-matrices 
T, one then gets a stochastic process on the space of 
transfer matrices. It is assumed in the derivation of the 
DMPK equations that the distribution of T is fully in- 
variant with respect to the subgroup K (isotropy assump- 
tion), which corresponds to a complete mixture of chan- 
nels at each elementary step. (This model assumption 
is justified by the fact that the mixing of channels in a 
realistic wire takes place on a scale much shorter than 
the localization length.) As a result, the stochastic pro- 
cess develops in fact on the coset space G/K, describing 
a Brownian motion on this manifold. In view of the ro- 
tational symmetry, the distribution function depends on 
radial coordinates Xi only. The corresponding evolution 
equation — which is the DMPK equation — has the form 



dL ~ 2ej ^ dxi ^'^' dxi ^"^^ 

2—1 



(5.4) 



where the Jacobian J{x) of transformation to the radial 
coordinates and the parameter 7 are specified below. 

For the Wigner-Dyson and Bogoliubov-de Gennes 
classes the radial coordinates Xi satisfy^ < a;i < a;2 < 



® For the classes D and Dili the proper variation domain of Xj is 
slightly different, |a:i| < a'2 < . . . < xpf llGruzberg et a^.ll2005l ). 
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. < and the Jacobian J{x) is 

N N 

j{x) = Y[ nisinh(x,±xj)r° ni^i"^2.Tfcr 

i<j lb k 



N 



Yl I sinh 



Xl\ 



Here mo, rUg, and mi are the multipUcities of the ordi- 
nary, short , and long roots for th e corresponding s y mmet - 
ric space (|Caselle and Magneal 12004 iHelgasonl . Il978[) . 
The root multipUcities for all the symmetry classes are 
listed in Table |lVl The first factor in Eq. ((5?5|) is re- 
sponsible for the repulsion between eigenvalues x^, anal- 
ogous to the energy level repulsion in RMT. The last 
two factors govern (when the corresponding multiplici- 
ties are non-zero) the repulsion between an eigenvalue Xi 
and its mirror —Xi, and the repulsion between Xi and 
zero, respectively. For the chiral classes, the variation 
domain of the coordinates Xi does not restrict their sign, 
xi < X2 < ■ ■ . < xn , and the Jacobian is 



N 



— W I sinh(a;i — Xj)Y 



(5.6) 



(For these classes = mi =0.) For the Wigner-Dyson 
and the chiral classes the multiplicity mo of the ordinary 
roots is the familiar parameter (3 of the RMT, equal to 
1,2, and 4 for the orthogonal, unitary, and symplectic 
symmetry classes, respectively. Further, the parameter 
7 in the DMPK equation is given by 



mo{N - \) + mi 4 
i[2 + m„(A^-l)] 



1 



WD and BdG 
chiral . 



(5.7) 



In the short-wire regime, L <^ 7^, where the conduc- 
tance is large, g ^ 1, the DMPK equation yields Ohm's 
law for the average conductance, the quasi-lD universal 
conductance fluctuations, and the weak- localization cor- 
rections in the form of a \/ g series. Our interest here is 
in the opposite, long-wire limit, L ^ 7^, where localiza- 
tion or critical properties of the problem manifest them- 
selves. This will be the subject of the remaining part 
of Sec. |Vl We also briefly mention that at the crossover 
scale, L ^ 7^, an analytical treatment is most compli- 
cated. An approximate sch e me has, however, been deve l- 
oped inlGopar al\ j2002l ): lMuttalib and Wolfli (|l999l ): 
iMuttalib et al. ( 20031 ) that allows one to understand key 
properties of the conductance distribution in this regime. 



B. Conventional localization in ID geometry 

We begin by considering the standard quasi-lD lo- 
calization in Wigner-Dyson classes (m; — 1, mo — f3, 
ms =0). In the long- wire limit the transmission eigen- 
values satisfy the hierarchy 1 ^ Ti ^ T2 ^ . . ., i.e. 



Ham. 


transfer matrix 


tr.matr. 




mi 


ms 


class 


symmetric space 


class 









Wigner-Dyson classes 





A 


V{p + q)/\J(p)x\J{q) 


AIII 


2 


1 


2|p-q| 




AI 


Sp(2iV, IR)/U(iV) 


CI 


1 


1 





(5.5) 


All 


SO*(2Ar)/U(iV) '^^^^'^'^ 
N odd 


Dlll-e 
DIII-o 


4 


1 




4 



chiral classes 



AIII 


GL(iV, C)/U(iV) 


A 


2 








BDI 


GL(iV, IR)/0(iV) 


AI 


1 








CII 


U*(2Af)/Sp(2iV) 


All 


4 








Bogoliubov - de Gennes classes 


C 


Sp(2p,2g)/Sp(2p)xSp(2g) 


CII 


4 


3 


4|p-g| 


CI 


Sp(2iV, C)/Sp(2Af) 


c 


2 


2 





BD 


0{p,q)/0{p)xO{q) 


BDI 


1 







Dili 


SO(iV,C)/SO(iV) ^ '''''''' 
N odd 


D 
B 


2 






2 



TABLE IV Transfer matrix spaces. First column: symmetry 
class of the Hamiltonian. Second and third columns: sym- 
metric space for the transfer matrix and the corresponding 
Cartan symmetry class. Last three columns: multiplicities of 
the ordinary (mo), long (mi) , and short (ms) roots. 



the consecutive Xi are separated by intervals much larger 
than unity. One can thus approximate the hyperbolic 
sine functions in Eq. (|5.5p by the exponentials. As a re- 
sult, all the variables in the DMPK equation (|5.4[) fully 
decouple; the resulting Fokker-Planck equation for each 
of Xk reads 



dVjxk) 
dL 



1 d^V 
2^94 



1 dp 

S.k dxk 



(5.8) 



where ^^T"^ = [1 -I- (3{k — l)]/7^. This is an equation of 
the advection-diffusion type, with playing the role 

of the diffusion constant and l/^i of the drift velocity. 
The solution has a Gaussian form with {xk) — L/^k and 
var(xfc) = L/^i. This immediately implies that the log- 
arithm of the conductance g ^ dj cosh^ x\ has a Gaus- 
sian distribution w ith the following average and variance 
(lBeenakkeilll997D : 



- (In.g) = 2L/7£ ; var(ln.9) 4i/7^ . (5.9) 

On the side of atypically large g this distribution is cut 
at 5 ~ 1. The average conductance is straightforwardly 
found to be determined by this cutoff (i.e. by rare events 
of g ^ 1), — ln(g) = L/27L Defining the typical and the 
average localization length via 



- (ln.g) =2L/Ctyp ; 
we thus find 

^typ = 7^ ; 



-ln(5) = 2L/Ca 



Cav = 47^ 



(5.10) 



(5.11) 
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These results are in full aRree ment with tho se obtained 
within the cr-model formalism (lMirlinl . [2000bh . which cor- 
responds to 3> 1, so that 7 = (3N. 

A similar behavior is als o found in the Bog'oliubov- 
de Gennes classes C and CI ( Brouwer et aLl . l2000 a). The 
same derivation yields the equation ()5.8p with ^^T^ = [3 + 
4(fc - l)]hl for class C and ^ = [2 + 2(fc - l)]ht for 
class CI. Thus, 

-(Ing) = 6L/7£ ; var(ln.g) 4L/7£ (C) ; (5.12) 
-(Ing) = 4L/7£ ; var(lng) = 4L/7£ (CI) . (5.13) 

Calculating the average conductance, one gets — ln(5) — 
AL/^l (class C) and —\n{g) — 2L/^l (class CI), so that 



6yp - 

^typ 



: 7^/3 ; 
= 7^/2 



= 7^/2 
V =7^ 



(C) ; 
(CI) 



(5.14) 
(5.15) 



A quasi-lD model of c lass C was also studied w ithin the 
(T-model formalism in iBundschuh et al\ (|l998t ). It was 
found there that ^av is 8 times shorter than for class A. 
This agrees with the result of the DMPK approach, as is 
seen by comparison of ^av in Eq. (|5.14p and Eq. (|5.1ip . 



C. Types of delocalization in disordered wires 



Having discussed in Sec. IV. Bl how conventional local- 
ization in disordered quasi-lD systems takes place, we 
now analyze how electrons in such a system can escape 
exponential localization. Table IIVI is very useful in this 
respect, showing that there are two distinct mechanisms. 

First, we see that in five classes (in all three chiral 
classes AIII, BDI, and CII, as well as in the supercon- 
ducting classes with broken spin-rotation invariance, BD 
and DIII) the multiplicity mi of long roots in zero. This 
implies that there is no repulsion between the eigenvalue 
Xi and its mirror — x^, so that the smallest (by absolute 
value) Xi can be close to zero. We will analyze the crit- 
ical behavior that emerges in these classes in Sec. IV. El 
andlVJ] 

The second type of delocalization is related to the mul- 
tiplicity of short roots wis in the Table IIVI In the con- 
ventional situation it is equal to zero. However, there are 
five classes (A, C, BD, AH, and DIII), where it can be 
non-zero. A non-zero value of implies a repulsion of 
Xi from zero, which is an indicator of the existence of ex- 
actly zero eigenvalues of the transfer-matrix. These zero 
eigenvalues imply perfectly transmitting channels, yield- 
ing a finite conductance in the limit of infinite system 
length. We will discuss this class of systems in Sec. IV. Dl 



D. IVIodels with perfectly conducting channels 

In this subsection, we consider the models with per- 
fectly transmitting channels, in which case ^ 0, see 
Table IIVI These models are in turn subdivided in two 
types. 



In the classes A, C, and BD, the transfer-matrix spaces 
given in the Table HVl correspond to a model with p chan- 
nels propagating to the left and q channels propagating 
to the right. While in the conventional situation p — q, 
symmetric spaces with p ^ q are allowed as well. It is 
not difficult to understand what is the physical realiza- 
tion for these models, if one recalls that these are exactly 
those classes whose 2D representatives show the quan- 
tum Hall effects, Sec. I VII When such a 2D system is on 
the quantum Hall plateau of order p, there are p edge 
channels that propagate on its boundary. These chan- 
nels are chiral in the sense that they can propagate in 
one direction only. The edge of such a system represents 
thus a wire of the corresponding symmetry class with p 
modes propagating in one direction and zero in the oppo- 
site. Since there is no backscattering, the conductance of 
such a wire is identically equal to p. Let us consider now 
parallel edges of two quantum Hall systems with coun- 
terpropagating modes separated by a potential barrier. 
Assuming that there are p modes in one edge and q in 
the other and that they are coupled by tunneling, we get 
a wire of the {q,p) type from the corresponding class. 
In the long length limit, \p — q\ modes will then remain 
perfectly conducting, while the rest will be localized. 

The situation with the classes AH and DIII is much 
more intricate. The systems of these classes may possess 
a single perfectly conducting channel. As we explain be- 
low, this reflects an underlying Z2 topological structure. 
A delocalization in the symplectic Wigner-Dyson class 
AH was obta ined for the first ti m e within t h e a-m odel 
formahsm in iMirhn et all (|1994 : IZirnbaued (|1992D . It 
was found in these works that the average conductance 
and its variance remain finite in the long- wire limit, 
(g) 1/2, var(g) -^1/4 due to a zero mode of the corre- 
sponding transfer-operator. The physical significance of 
these results was not understood at this s tage. Further, it 
was shown in lBrouwer and Frahml (|l996D that the above 
zero mode is double- valued on the cr-model manifold and 
thus does not contribute in the case of a conventional wire 
with spin-orbit interaction. More recently, it was under- 
stood, however, that a model of symplectic symmetry 
with a perfectly conducting channel arises if one consid- 
ers tr ansport i n carbon nanotub es ( Ando and Suzuural . 
l2002t feuzuura and Andol . 12002^ . The problem is de- 
scribed by two species of Dirac fermions corresponding to 
two valleys in the graphene spectrum. If the scatterers 
are of long-range character and the intervalley scatter- 
ing can be neglected, the problem acquires the symplec- 
tic (AH) symmetry, with the sublattice space taking the 
role of isospin (Sec. IVI.G.2p . Furthermore, in contrast to 
conventional wires of All symmetry, where the number 
of channels is even, there is just a single channel here. 
(More precisely, its "partner" belongs to the other node, 
and they "do not talk to each other".) As a result, the 
channel remains perfectly transmitting independently of 
the length of the wire. 

In subsequent works (ICase lle and Magneal l2006t 
ISakai and Takand . l2005l : iTakand . l2004a.b ) quasi-lD sys- 
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terns of the All symmetry class with odd number N of 
channels were studied within the DMPK approach. It 
was found that for any odd N a single perfectly trans- 
mitting channel remains in the limit L ^ ji. As empha- 
sized in iTakana ()2004ciV these results are i n full agree- 
ment with the earlier g-model results of iMirlin et all 
(|l994D : IZirnbauerj(|l992D . if the latter are complemented 
by the follo wing interpretation. The Fourier expa nsion 
employed in iMirlin a/] (| 19941 ): IZirnbauerj (|l992[) con- 
tained eigenfunctions of the Laplace operator in the a- 
model manifold of two types - with even and odd parity 
(the zero mode is of the odd-parity type). For systems 
with even (odd) number of channels one should keep only 
even-parity (resp. odd-parity) eigenmodes ar id include an 
overal l factor of two. T he original result of iMirlin et all 
(|l994l ): IZirnbaued(fl992) with (a), (g^) ^ 1/2 at L ^ oo 
corresponds thus to an average over wires with even and 
odd number of channels. 

In the problem with perfectly transmitting channels, 
one can determine the localization length for the remain- 
ing modes by considering the deviation Sg of the conduc- 
tance from its L — + cxD limit. A straightforward general- 
i zation of the consideration sketched in Sec. IV. Bl yields 
( Caselle and Magneal [20061 ) 



(InSg) = (2m; -I- ms)L/j£ 



(5.16) 



so that ^typ = 7^/(wi + ins/2). This implies for both the 
All and Dili models with an odd number of channels 



typ 



7^/3; 



j£/2 



(5.17) 



The qualitative different behavior of class-All (and 
Dili) wires with even and odd number of channels is 
intimately connected with a non-trivial topology of the 
corresponding cr-model manifold (or, more specifically, of 
its compact component Mf)'- the first homotopy group 
■7Ti{Mf) is equal to Z2. T his enables a topologica l 6- 
term with 6 equal to or tt ( Ostrovskv et all . l2007al) , in 



analogy with the 2D situation, see Sec. IVI.A.Si IVI.B.5I 
for these symmetry classes. The topological term with 
6 ^ TT is present if the number of channels is odd. 

Another realization of a symplectic-symmetry wire 
with an odd number of channels has recently emerged 
in the context of the quantum spin Hall (QSH) ef- 
fect in systems of Di r ac fermions with spin - orbit c ou- 
pling (jBernevig et all . 120061: iKane and Meli . l2005al lbl). 
Such systems were found to possess two distinct insu- 
lating phases (with a transition between them driven 
by Rashba spin-orbit coupling strength), both having a 
gap in the bulk electron spectrum but differing by the 
edge properties. The topological distinction between the 
two insulating phases ret ains its validity in the pres- 
ence of disorder (.Essin and Moore ^ 20071: Obuse et all . 



l2007al : IOnoda et aZ.1 . 120071: ISheng et a^.l ■ l2005^ . While the 
normal insulating phase has no edge states, the QSH 
insulator is characterized by a pair of mutually time- 
reversed spin-polarized edge states penetrating the bulk 
gap. These edge states in the QSH phase, which do not 



get localized by disorder, represent a realization of the 
class- AH wire with a non-trivial Z2 topology. 



E. Chiral classes 

For the chiral classes , the cons i deratio n analogous to 
that in Sec . ITbI (Brouwer et all . l2000bl : iMudrv et all . 
I1999L I2OOOD yields {xk) = (2fc - 1 - N)PL/2je. If the 
number of channels is even, then the smallest by abso- 
lute value eigenvalues Xk are separated by a large gap 
from zero: —{xN/2) = {x N/2+1) = L/2j£. Therefore, the 
exponential localization is preserved, —{hug) = f3L/"f£, 
vai{\ng) = AL/^£, yielding the localization lengths 



6yp = 2-f£/l3 



167£//32 



(5.18) 



On the other hand, if TV is odd, one of the eigenvalues is 
close to zero, {x(^m+i)/2) = 0. This leads to a completely 
different behavior of the conductance. 



{g) = (27£/7rL)V2 



1/2 



var(lng)=(^4-^j A;(5.i9) 
var(g) = {8j£/9TTLy^^ . (5.20) 



It is seen from Eq. (|5.19p that while the typical conduc- 
tance decays in a stretched-exponential way, its fluctu- 
ations are very strong, so that the probability to have 
g ~ 1 is small as L^^l"^ only. This determines the very 
slow decay (|5.20p of the average conductance, which is 
even slower than the in classical Ohm's law. 

The delocalization takes place for arbitrary odd A^, in- 
cluding TV = 1. The single-channel model with chiral- 
class disorder has been studied, in its various incarna- 
tions, in a large n umber o f work s, starting form the pi- 
oneering paper by iDvsonI 



1951 . We list most salient 



features characterizing (in addition to the above results 
for the statistical properties of the conductance) this crit- 
ical point. 

i) Localization length. If the energy E deviates from 
zero, the chiral symmetry is broken, and the exponential 
localization establishes. One can thus ask how the cor- 
responding localization length diverges at ^ 0. It has 
been found that one should distinguish between average 
and typical observables (e.g. conductance) whose spatial 
depende nce is governed by two parametrically different 
lengths (|Balents and Fishej . Il997l: iFisherj . Il995l ) , 



typ 



|ln£;| 



Ca 



ln£;| 



(5.21) 



ii) Staggering. An alternative way to drive the system 
out of criticality is to introduce a staggering M in the 
hopping strength which opens a gap around zero energy 
in the spectrum of a clean system. The co rresponding lo- 
calization lengths behave as follows (fBalents and Fisherl . 
,1997, : .Fisher . 1995; Mathur, 1997'), 



Ctyp ^ M- 



M' 



(5.22) 
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Hi) Wave function at criticality. The Hamiltonian of 
a single-chain problem ca n be written in a Dirac form 
(|Balents and Fished . [iQQTh . 



H 



-iozdx + m{x)ay 



(5.23) 



where m(x) = M -\- Th{x) and 771(2;) is the disorder (e.g. 
of the white-noise type) . The zero-energy eigenfunction 
can then be found explicitly, 




(5.24) 



The properties of t his w ave function were analyzed in 
iBalents and Fisherl (1993). The following scaling of the 
spatial correlation function of the moments of ipix) at 
criticality {M = 0) was found, 



(|V(x)v(0)|^ 



L-'\x\-'/' 



(5.25) 



for all g > 0. This result can be interpreted in terms 
of the following pic t ure o f wave functions at criticality 
(IBalents and Fished . Il997t) . The wave function is typi- 
cally quasi-localized, showing a stretched-exponential de- 
cay with respect to its principal maximum. However, 
with a probability ~ x~^^^ it shows a secondary maxi- 
mum of a magnitude close to the primary one and sepa- 
rated from it by a distance x. 

iv) Density of states. The DOS s h ows at cri t icality 
the Dyson singularity ( Dvsonl . Il953l : iMcKenzi^ . Il996t 
iTitov a/.. . ,200 1). 



p{E) - l/\E\n^E\ 



(5.26) 



When the system is driven away from criticality by a non- 
zero staggering parameter M, the singularity weakens 
and becomes non-universal, p{E) ~ with 6 > 0. 

Finally, it has been shown that a sufficiently strong 
staggering M c an also drive a system with even N into 
a critical state ( Brouwer et aZ.l .ll9981. More specifically, 
the staggering shifts all the variables Xk by a constant. 
With increasing M, the average values (xk) consecutively 
cross zero; whenever this happens, the system is at criti- 
cality. Therefore, whether N is odd or even, changing M 
will drive the system through N transition points. 



F. Bogoliubov-de Gennes classes with broken spin-rotation 
invariance 

Analysis of the DM PK equation for the classes ED and 
Dili (Brouw er et Ml, [2000a) leads to results identical to 
those obtained for chiral classes, Eq. (|5.19p . (|5.20p . Fur- 
thermore, the DOS was found to show the Dyson singu- 
larity (|5.26p . again in fu ll analogy with the chira l classes. 

On the other hand, iMotrunich et al. 1 (|200lh studied 
certain single-channel models of the classes D and Dili 



via a strong-disorder real-space RG. It was found that 
generically these systems are in localized phases and the 
DOS diverges in a power-law fashion with a non-universal 
exponent, p{E) ~ E~^^^ with 5 > 0. Only at phase 
boundaries the system is critical and the DOS takes the 
Dyson form (|5.26p . 

An apparent contradi ction between th e resul ts of both 
papers was resolved in iGruzberg et al. I dlooi)- 

It was 

shown there that, generically, quasi- ID systems of the 
class es BD and Dili ar e in a l ocalized phase, in agreement 
with IMotrunich et ai\ (|200ll ). The terms that drive the 
system towards localization are usually neglected within 
the DMPK approach, as they are irrelevant at the short- 
distance (diffusive) fix ed point. It was found, however, in 
iGruzberg et al. I (l2005h that these terms become relevant 
at the long-distance (critical) fixed point and drive the 
system away from it, into the localization fixed point. 
Only if the disorder is fine tuned, the system is at the 
critical point. On the other hand, the length at which the 
crossover from criticality to localization happens becomes 
exponentially large with increasing number of channels 
N . Therefore, in the thick-wire limit, 1, the system 

is essentially at criticality. An analo gous conclusion was 
also reached in lBrouwer et ^ (|2003D . 

It was also argu ed in IGruzberg et al\ (|2005'); 
IMotrunich et al\ ( 200lh that critical points with Dyson 
singularities of all five symmetry classes (AIII, BDI, 
CII, BD, and Dili) belong to the same universality 
class. To estab l ish th is remarkable "superuniversality" , 
IGruzberg et al\ (|2005l) pointed out that all the universal 
properties can be obtained from A'^ = 1 models and then 
constructed mappings between single-channel models of 
all the five classes. 



VI. CRITICALITY IN 2D 



A. Mechanisms of criticality in 2D 



As was discused in Sec. III.B.2[ conventional Ander- 
son transitions in the orthogonal and unitary symmetry 
classes take place only if the dimensionality is d > 2, 
whereas in 2D all states are localized. It is, however, 
well understood by now that there is a rich variety of 
mechanisms that lead to emergence of criticality in 2D 
disordered systems. Such 2D critical points have been 
found to exist for 9 out of 10 symmetry classes, namely, 
in all classes except for the orthogonal class AL A nice 
su mmary of possi ble types of 2D criticality was given 
in ( Fendlevl . I2OOOI ): we closely follow this work in our 
presentation in this subsection. We now list and briefly 
describe the mechanisms for the emergence of criticality; 
a detailed discussion of the corresponding critical points 
will be given in the remaining subsections of Sec. IVII 
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1. Broken spin-rotation invariance: IVletallic phase 

We begin with the mechanism that has been aheady 
mentioned in Sec. III. B. 21 in the context of the Wigner- 
Dyson symplectic class (All). In this case the /3-function 
[ (|2.25p with e = 0] is positive for not too large t (i.e. suffi- 
ciently large conductance) , so that the system is metallic 
(t scales to zero under RG). On the other hand, for strong 
disorder (low t) the system is an insulator, as usual, i.e. 
P{t) < 0. Thus, /3-function crosses zero at some t*, which 
is a point of the Anderson transition. Properties of this 
critical point will be discussed in detail in Sec. IVI.BI 

This mechanism (positive /3-function and, thus, metal- 
lic phase at small t, with a transition at some t*) is also 
realized in two of Bogoliubov-de Gennes classes - D and 
Dili, see Table IIIII All these classes correspond to sys- 
tems with broken spin-rotation invariance. The uncon- 
ventional sign of the /3-function in these classes, indi- 
cating weak antilocalization (rather then localization), is 
physically related to destructive interference of time re- 
versed paths for particles with spin s = 1/2. 



2. Chiral classes: Vanishing /3-function 

Another peculiarity of the perturbative /3-function 
takes place for three chiral classes - AIII, BDI, ad GIL 
Specifically, for these classes /3(t) = to all orders of 
the perturbat ion theory, as was first discovered by Gade 
and Wegner (iGadd . Il993l : iGade and Wegned . Il99lh . As 
a result, the conductance is not renormalized at all, serv- 
ing as an exactly marginal coupling. There is thus a 
line of critical points for these models, labeled by the 
value of the conductance. In fact, the cr-mod e ls for 
these classes contain an additional term (iGadel . 119931 : 
iGade and Wegneil . Il99ll ) that does not affect the absence 
of renormalization of the conductance but is crucial for 
the analysis of the behavior of the DOS. A detailed dis- 
cussion of the chiral classes will be given in Sec. IVI.FI 



3. Broken time-reversal invariance: Topological ^-term and 
quantum Hall criticality 

For several classes, the cr-model action allows for inclu- 
sion of a topological term, which is invisible to any order 
of the perturbation theory. This is the case when the sec- 
ond homotopy group tt2 of the cr-model manifold M (a 
group of homotopy classes of maps of the sphere into 
A^) is non-trivial. From this point of view, only the 
compact sector A4 p (originating from the fermionic part 
of the supervector field) of the manifold base matters. 



There are five classes, for which 7r2(A^F) is non-trivial, 
namely A, G, D, All, and GIL 

For the classes A, G, D the homotopy group tt2 {M f) = 
Z. Therefore, the action S[Q] may include the (imagi- 
nary) 0-term, 



iStop[Q]^iON[Q] 



(6.1) 



where an integer N[Q] is the winding number of the field 
configuration Q{r). Without loss of generality, 6 can be 
restricted to the interval [0, 27r], since the theory is peri- 
odic in with the period 27r. 

The topological term (|6.ip breaks the time reversal 
invariance, so it may only arise in the corresponding 
symmetry classes. The by far most famous case is the 
Wigner-Dyso n unitary class (A). A s was first understood 
by Pruisken (jPruiskenl 1 1984 Il987t) . the a-model of this 
class with the topological term (|6.ip describes the inte- 
ger quantum Hall effect (IQHE) , with the critical point of 
the plateau transition corresponding to 9 = n. More re- 
cently, it was understood that counterparts of the IQHE 
exist also in the Bogoliubov-de Gennes classes with bro- 
ken time- reversal invariance - classes G and D. They 
were called spin and i/iermaZ quantum Hall effects (SQHE 
and TQHE), respectively. The criticahty at the IQHE, 
SQHE, and TQHE transitions will be discussed in detail 
in Sec. lVLCllVLDl and lVLEl respectively. 



4. Z2 topological term 

For two classes, AH and GII, the second homotopy 
group is tt2{.Mf) = Z2. This allows for the 6'-term but 
9 can only take the value 9 = and 9 = n. It was 
shown very recently ( Ostrovskv et aLl . l2007al ) that the cr- 
model of the Wigner-Dyson symplectic class (AH) with 
a, 9 = TT topological angle arises from a model of Dirac 
fermions with random scalar potential, which describes, 
in particular, graphene with long-range disorder. Like in 
the case of quantum-Hall systems, this topological term 
inhibits localization. Whether the model then flows un- 
avoidably into the ideal-metal fixed point or, else, there 
is also a novel attractive fixed point is a matter of ongo- 
ing research. The reader is referred to Sec. lVI.Bj for more 
detail. 



A pedagogical introduction of topological concepts in the context 
of conde nsed matter theory can be found in the recent mono- 



graph bv lAltland and Simon 



5. Wess-Zumino term 

Finally, one more mechanism of emergence of criticalty 
is the Wess-Zumino (WZ) term. It is known that this 
term may appear in cr-models of the classes AIII, GI, 
and DHL For these classes, the compact component A4 p 
of the manifold is the group H x H/H = H, where H is 
U(ri), Sp(2n), and 0(2n), respectively. The correspond- 
ing theories are called "principal chiral models" . The WZ 
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term has the following form: 

ik 



iSwzig) 



247r 



r I ds e 





xStiig~'d,g){g-'d,g){g-'dxg), (6.2) 

where k is an integer called the level of the WZW model. 
The definition (|6.2p of the WZ term requires an extension 
of the cr-model field ^(r) = g{x, y) to the third dimension, 
< s < 1, such that g(r, 0) = 1 and g(r, 1) — g{r). Such 
an extension is always possible, since the second homo- 
topy group is trivial, Tr2{H) — 0, for all the three classes. 
Further, the value of the WZ term does not depend on 
the particular way the extension to the third dimension 
is performed. (This becomes explicit when one calculates 
the variaton of the WZ term: it is expressed in terms of 
g{r) only.) More precisely, there is the following topo- 
logical ambiguity in the definition of Sy/zig)- Since the 
third homotopy group is non-trivial, 7r3(if) = Z, 5'wz(ff) 
is defined up to an arbitrary additive integer n times 27rfc. 
This, however, does not affect any observables, since sim- 
ply adds the phase nk x 27rz to the action. 

The WZ term arises when one bos onizes certain mod- 
els of Dirac fermions llWittenl . 1X9841) and is a manifes- 
tation of the chiral anomaly. In particular, a cr-model 
for a system of the AIII (chiral unitary) class with the 
WZ term describes Dirac fermions in a random vector 
potential. In this case the u-model coupling constant 
is truly marginal (as is typical for chiral classes) and 
one finds a line of fixed points. On the other hand, for 
the class CI there is a single fixed point. The WZW 
models of these classes were encountered in the course 



of study of dirty d-wave supe rconductors (jAltland et all . 
I2QO2I : iNersesvan et all . Il995l ) and, most recently, in the 
context of disordered graphene. We will discuss critical 
properties of these models in Sec. IVI.G.3I 



B. Symplectic Wigner-Dyson class (All) 

In metals with spin-orbit coupling the spin of a parti- 
cle is no longer conserved. The spin-up and spin-down 
channels are coupled, and an electron needs to be repre- 
sented as a two component spinor. If the time-reversal 
symmetry is preserved, the system belongs to the sym- 
plectic Wigner-Dyson class AIL The one-loop quantum 
correction at large conductance g takes then the form of 
weak antilocalization, see Sec. III. B. 21 IVI.A.ll At lower 
g the one-loop /3-function is not sufficient anymore, and 
higher-order terms lead to localization, with the Ander- 
son transition at some . While the /3-function has been 
calculated up to the four-loop order, see Eq. (|2.25p . this 
does not help to get quantitative predictions for critical 
properties. In particular, an attempt to use the four-loop 
/3-function to extract the localization length exponent 



(IWegned . \l9M) yields (|C(3))^^^ wO.193, which is 
an order of magnitude smaller than the numerical re- 
sult (see below) and even violates the Harris criterion 



V >2/d (jChaves et al\ . Il986l ). This is not very surpris- 
ing: the considered Anderson transition takes place at 
strong coupling, ^ 1, so that keeping just the first few 
terms of the perturbative expansion is an uncontrolled 
procedure. In this situation, numerical simulations are 
particularly important. On their basis, a detailed quan- 
titative picture of the transition has been developed; the 
key findings are summarized below. 



1. Microscopic models 

Most numerical studies employed a tight binding 
Hamiltonian. It is defined on a two dimensional square 
lattice with nearest neighbor coupling 



E 



,C- C ■ . 



(6.3) 



^^v.^^, ^ a) denote creation (annihilation) operators 
of an electron with spin cr on site i. The on-site energies 
are taken to be random numbers drawn from the interval 
[— W/2, VF/2] with a homogeneous distribution. There 
exist various versions of the model, that differ by the 

choice of the hopping matrix Vi^o-y.o-' ■ 

Most studies employ the Ando model ( Andol . Il989l ) 
characterized by non-random hopping between next 
neighbors only. 



yi,a-iJrk,a' = (Vb exp(l6'feO-/c))c 



k^x,y, (6.4) 



where ax,(Jy are the Pauli matrices. Convention- 
ally, the spin-orbit energy scale is set to unity, 
Vq—1, and the mixing angles take cons tant values 
9k=TT/Q. In the Evangelou-Zim an model (Evangelou|, 
ll995tlEvangelou and ZirrLaiil . ll987| ). the components of V 
that are proportional to <Jx,y,z^ are chosen to be random 
with prefactors drawn independently from a box distri- 
bution of a width Vn . Recently, a third variant was intro- 
duced (jAsada et all l2002l |2004 - the SU(2) model, in 
which the random matrix exp(i0fe(Tfe) is chosen to be uni- 
formly distributed on the SU(2) group. This last model 
was found particularly suitable for numerics, since finite 
size corrections appear to be very small. 



2. Localization length exponent 

The small magnitude o f finite-size c o rrections in the 
SU(2) model has allowed lAsada et a/] (|2002l l2004l ) to 
determine the localization length exponent with a high 
precision, z/=2.746 ± 0.009. Very recently, a similar re- 
sult (th ough with somewhat higher un certainty) was ob- 



^ r- W K 

tained (jMarkos and Schweitzer! . 12006!) in th e cont ext of 



the Ando model, v = 2.8 ±0.04. lAsada et a/] (|2004l ) eval- 
uated numerically the whole /3-function and found that 
it has the expected shape, with a single zero determining 
the transition point. 
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3. Critical conductance 

Since the transition takes place in the strong cou- 
pling regime, the mesoscopic conductance fluctuations 
at criticallity are comparable to the mean value, so 
that an ensemble of macroscopically identical coherent 
samples is to be characterized by the whole distribu- 
tion function Vig). Quite generally speaking, critical 
conductance distributions are scale-invariant but depend 
on the shape of the sample (similar to the IPR dis- 
tribution function. Sec. III. C. 51 and the level statistics, 
Se c. III.Ep; a revie w of numerical results has been given 
bv iMarkod (|2006l) . The distribution Vig) for a square 
sample a t the symplecti c Ande rson transition was deter- 
mined in lOhtsuki et al\ t004 for the SU(2) m odel, see 
Fig. 



and 



Markos and Schweitzeil ( 20061 ) for the 
Ando model, with essentially identical results. The av- 
erage value {g) was found to be (q)=l. 42 ± 0.005 with a 
variance v&rg— 0.36 (|Markos and Schweitzeil . I2006D . 



4. Multifractal spectrum 

The spectrum Tq can be calculated for the 2D sym- 
plectic transition with very good accuracy, because 
corrections to scaling turn out to be extreme ly sma ll 
(|Asada et al\ . 120041 : iMildenberger and Everl [2003). 
Therefore, many generic features of the critical wavefunc- 
tion statistics can be studied in great detail. 

The results of a high-precision study of wave func- 
tion multifractality at the symplect ic Anderson transi- 
tion ( Mildenberger and Eversl l2007l ) are summarized in 
Fig. [25l In order to highlight the non-trivial features, 
the reduced anomalous dimensions S, 
depicted. The key observations are as follows 

(i) One finds that 6o = uq - 2=0.172 ± 0.002. This re- 
sult may be used as a check on the conformal invariance 
which imposes the exact condition (|janle3,[i99l[l993) 



q=Aq/q{l-q) are 




0.17 - 



FIG. 25 Multifractal spectrum 5q for the Ando model 
(dashed, 14/^=5.84) and the SU(2) model (solid, 17^=5.953). 
To highlight deviations from parabolicity, reduced anomalous 
dimensions 5q = Aq/q{l — q) are plotted. Anomalous di- 
mensions S^q"^ obtained from typical IPR are also shown (o). 
Dashed lines indicate the estimated error {2a) in Sq. Inset: 
blow up of the solid line behavior near q=\ is represented 
by empty circles (o). Filled symbols (•) show original trace 
after reflection at q=\. Dot-dashed line is a fit (offset: 10"'^) 
5, = 0.1705 + 0.0043(g-i)^ (|Mildenberger and Eversl. |2007^ 



7r(5nAn = 1. 



(6.5) 



uncertainty in determination of Tq also increases, explain- 
ing some deviation between the data for both models. 

Very similar r esults for the multif ractality spectrum 
were obtained in ( Obuse et aZ.I . [20"07b[ ) . In this work, the 
multifractality was also studied at the boundary and at 
the corner of the system. It was found that the multifrac- 
tal exponents fulfill the relation (|2.59p . thus providing a 
further strong evidence for the conformal invariance at 
this critical point. 



Here Kc—^m/M, where is the localization length in a 
quasi- ID strip of width L at criticality. ) W ith the above 
value for 5q and with Ac=1.844 ± 0.002 (jAsada et all 
|2004| ). the left-hand side of Eq. (|6.5p becomes indeed very 
close to unity: 0.996 ± 0.013. 

(ii) The function 5q fulfills the symmetry relation 
Ea. (|2.38p . as demonstrated in the inset of Fig. [55] where 
the data for the range —0.5 < g < 1.5 are displayed. 

(iii) 5q has a small but non-zero curvature, implying 
that the the multifractal spectrum is not parabolic. 

(iv) The results are essentially identical for the SU(2) 
and Ando models, confirming the universality of the tran- 
sition. An abrupt change in the behavior of the data for 
the typical IPR in Fig. [25] at q ~ 2.5 is related to the fact 
that at g > when the average IPR probes the tail of 
its distribution function, the exponents Tq and t]^^ start 
to differ, see Eq. (I2.42p . In this range of q the statistical 



5. Symplectic-class theories with Z2 topology. 



As was explained in Sec. IVI.A.4I the cr-model of the 
symplectic class allows for an inclusion of a topologi- 
cal term with 9 = n. A microscopic realizations of 
such a n on-trivial topology was f or the first time iden- 
tified in ' Ostrovskv et aZ.I ( 2007al) where the model of 
Dirac fermions in disordered graphene was studied (see 
Sec. IVI.G)) . It was found that for the case of long-range 
impurities, when two valleys in the spectrum are decou- 
pled and the problem reduces to that of a single species 
of Dirac fermions in random potential, the field theory is 
the class- All cr-model with 9 ~ tt topological term. 

The fermionic sector of the corresponding cr-mode 
manifold is Mf — 0(4n)/0(2n) x 0(2n). In fact, for 
the "minimal" supersymmetric cr-model (n — 1) the sec- 
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ond homotopy group is richer: 



Tr2[Mi 



Tr2[Mi 



(6.6) 



For n — 1 the compact sector of the model is the 
manifold (S^ x S^)/Z2 (product of the "diffuson" and 
"Cooperon" 2-spheres divided by Z2). Thus two topo- 
logical invariants, Ni^2[Q], counting the covering of each 
sphere, emerge in accordance with Eq. (j6.6p . The most 
general topological term is iS'top = iOiNi + 16*2 A^2. How- 
ever, the time-reversal symmetry requires that the ac- 
tion is invariant under interchanging the diffuson and 
Cooperon spheres, which yields 9i = 02 = where 
is either or tt. Hence only a Z2 subgroup of the whole 
Z X Z comes into play as expected: the phase diagram 
of the theory sh ould not depend on n. For the Dirac 
fermion problem ( Ostrovskv et a/.l . l2007a[ ). an exphcit ex- 
pression for n = 1 topological term can be written using 
u = TVT- ^ (where the a-model field Q = T'^AT, see 
Sec. lH.B.ll) . 



tS2[Q]^'-^Str[{A±l)T2U^up] 



iTT{Ni[Q]+N2m, 



yielding = n. The non-trivial value of the topological 
angle (0 — tt) holds for higher n as well. It implies that all 
configurations Q{r) are subdivided into two topologically 
distinct classes ( "even" and "odd" ) ; the former give pos- 
itive and the latter negative contribution to the cr-model 
partition function . This was confirmed numerically in 
iRvu et all (|2007bl ). 

At large conductance g the contribution of topologi- 
cally non-trivial configurations is exponentially small and 
can not affect the metallic phase in any essential way. 
On the other hand, at small g the topological term is 
expected to suppress localization, similarly to its role in 
the quantum Hall effect fSec. IVI.C]) and in the quasi-lD 
symplectic model (Sec. IV. D[) . This leaves room to two 
possibilities: 

(i) The /3-function changes sign twice. This would 
mean that, in addition to the conventional repulsive fixed 
point of the symplectic class, a new attr active fixed point 
arises. This scenario was proposed in lOstrovskv et al\ 
(l2007aD . Then, if the RG flow starts with a sufficiently 
low conductivity, it ends up in this new critical point with 
a universal conductivity of order unity. 

(ii) The (3- function remains positive everywhere. The 
RG flow then necessarily leads the system into the ideal- 
metal fixed point with infinite conductivity. 

In view of the strong-coupling nature of the prob- 
lem, numerical simulations are needed to resolve 

this dile mma. Recent simulations of disordered 

graphenc ("Bardarson et al!, '2007*; 'Nomura et al!, '2007|; 
[Nomura and MacDonald, 2007; Rycerz et al, 2006) do 
confirm the suppression of localization in the sym- 
pl ectic class with Z2 topology. While the results 
of iNomura and MacDoimldl (|2007f ): iRvcerz et al\ (|2006f ) 
were consistent with the scenario (i), w ith a critical con- 
ductivity ^ e^ /h, most recent works ( Bardarson et al\ . 



l2007t iNomura et al\ . l2007t ) appear to favor the second 
scenario. 

It is worth reminding the reader of a different type of 
Z2 topology in 2D systems of the symplectic symmetry 
class (AH). It arises in the context of the quantum spin 
Hall (QSH) effect and is related to a non-trivial first ho- 
motopy group, T:i{Mp) — Z2. This enables, in full simi- 
larity to the 2D situation, a 0-term with equal to or 
TT also in ID case, inducing a Z2 topological classification 
of edge states in QSH systems, see Sec. IV.DI There- 
fore, these systems possess in addition to the metallic 
phase, two distinct insulating phases, with different edge 
properties (normal insulator and QSH insulator). An im- 
portant question is whether there is a direct, quantum- 
Hall- type tran sition between these two phases. Recent 
numerics dEssin and Moord . l2007l : lObuse et all . l2007al: 
lOnoda et al\ . 2007 ) on some models of QSH systems gives 
a negative answer: the insulating phases are found to be 
everywhere separated by the metallic phase. It would be 
interesting to find out whether such a direct transition is 
generically prohibited, independently of the microscopic 
model. More activity in this direction may be expected 
in near future. 



C. The integer quantum Hall effect 

Our presentation in this section compleme nts the re- 



Uur presentation m tnis section compieme i 
views ( Huckesteinl . Il995l : iKramer et al. I. (20051) . 



1. Pruisken's cr-model 



As was discovered in lv. Klitzing et all ( IQSOl ). the Hall 
conductivity Oxy of a 2D electron gas in a strong trans- 
verse magnetic field develops plateaus at values quantized 
in units of e^ /h. While the physics of the Hall plateau 
is fairly well understood by now, the theory of the quan- 
tum critical points separating the plateaus - the quantum 
Hall transition - remains a challenging issue. 

From the field-theoretical point of view, the IQHE 
is described by the cr-model (j2.17p with a topological 
term (j6.ip . It w as first deriv e d by Pruisken i n the 
replic a formahsm ( Levine et ali Il983l : iPruiskeiJ . 1 1984 
19871): a supersymm etric generalization was obtained in 
Weidenmiilleij ( 19871 ). The action of the model reads 



S[Q] = istrhc7,,(VQ)2 + 2axyQVxQVyQ], (6.7) 

where a^x, f^xy are dimensionless conductivities. The 
Hall conductivity (Jxy is related to the topological an- 
gle as axy = 0/2tt. There is strong evidence that 
the corresponding two-parameter fiow dia gram has the 
form s hown in F ig. as proposed in iKhmelnitskiil 
(|l984l ): iPruiskenl (|l985rfl987l ). The fixed point at 6* = 
(2n -|- l)7r describes the QH transition. While the the- 
ory (j6.7p is highly important for understanding of qual- 
itative features of the problem, it allows to make only 
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(e2/h) 



FIG. 26 Two-parameter fl ow diagram of the Pruisken a- 
model, as first proposed in (jKhmelnitskii ITqS^ I . 



rou gh predictions for parame t ers of the critical behav- 
ior ( Pruisken and Burmistrovl l2005t ). This is because a 
controllable calculation in this framework can only be 
performed at weak coupling, <Txx S> 1, while the fixed 
points are at strong coupling, Gxx ~ 1. In this situation, 
numerical simulations are particularly important; their 
results will be reviewed in Sec. IVI.CTsUVI.C.Tl 



2. Further analytical approaches 

A great deal of effort has been invested by many re- 
searchers in order to attack the problem of the QH transi- 
tion from the analytical side. In addition to the Pruisken 
model, Sec. IVI.C.li several other analytical frameworks 
have been used. While this activity has not led to an ul- 
timate success in the quantitative description of critical 
behavior, a variety of important connections between the 
models has been established. In particular, it has been 
shown that the a-model (|6.7|) is also obtained as a contin- 
uum limit of the Chalker-Coddington network described 
in Scc lVI.C!4l ( Zirnbauer., .1997, ). Further, either of these 
two models can be map ped onto a quantum antiferromag 
netic superspin chain ([Kon dev and Marston'. '1997|; iLee 



11994 iMarston and Tsail . Il999: Zirnbauer, 1994. Tl 997 
Unfortunately, attempts to find an integrable deforma- 
tion of this spin chain have failed. A further approach to 
QH criticality is based on the model of Dirac fermions; 
it will be reviewed in Sec. IVI.G.2I 



is promoted to the group SU(2) (isomorphic to the 3- 
sphere S^) at criticality. The idea is thus to identify 
the corresponding critical theory for the QH problem, 
with a hope that it is of the WZW type and is solv- 
able by the methods of t he conformal field theory. Such 
a proposal was made in IZirnbaue"il ( 1999t) . along with 
a detailed analysis of constraints on the sought fixed- 
point theory. The target space of the theory conjec- 
tured by Zirnbauer is a real form of the complex super- 
group PSL(2|2). Its base A^_f x Mb is a product of 
the 3-sphere Mf = SU(2) = S^ and the 3-hyperboloid 
SL(2, C)/SU( 2) = H^. A model of the same type was also 
proposed in (iBhaseen et a/] . |2000| ) and most recently in 
(|Tsvelikl . l2007i r The proposed theories have the form of 
the WZW model, see Sec. IVI.A.5I 



(6.8) 



where i^wz is the WZ term (|6.2p . The peculiarity of the 
WZW models on the considered manifold is that they 
are critical at any val ue of the couplin g constant t and 
lev el fc € N. While (IZirnbaueii Il999l) argues for k = 
1, (jBhaseen et al considers the model with Kac- 

Moody symmetry, k — 1/t. The later condition restricts 
1 /t to be integer but facilitates t he analy s is of t he model. 
Very recently, it was proposed ( Tsvelikl . l2007f ) that the 
later model with fc = 8 may be the required fixed-point 
theory. 

Both variants of the theory make a prediction for 
the statistics of critical eigenfunctions. Specifically, it 
is found that the multifractality spectrum is exactly 
parabolic, 



A, = 79(1 - q) ; 

/(a)-2-(a-ao)V4(«o-2); 



(6.9) 

ao = 2 + 7, (6.10) 



with ao — 2 = At in the case of (IZirnbaiieil Il999[ ) and 
ao — 2 = 2i for ( Bhaseen et all . I2OOOI ) . This prediction of 
parabolicity of A, and /(a) indeed agrees with numerical 
simulations, Sec. IVI.C.7[ supporting this type of models. 
However, many questions related to the above conjec- 
tures remain open. In particular, if there is a whole line 
of fixed points (parametrized by t) , then is there univer- 
sality at the QH transition? If yes, how is it established? 
From the numerical point of view, there is no indication 
of non- universality at present. 



3. Quest for conformal field theory 



4. Chalker-Coddington network 



Another line of activity is the search for a conformal 
field theory of the QH transition. The guiding principle is 
related to the fact that a relative of the Pruisken's model, 
the 0(3) (T-model with 9 — tt topological term, describ- 
ing a ID quantum antiferromagnet with half-integer spin, 
flows under renormalization to a SU(2) WZW model. 
This means that the target space - which is the 2-sphere 
0(3)/0(2) = SU(2)/U(1)=S2 for the 0(3) cr-model - 



The Chalke r-Coddington network (CCN) m odel was 
introduced in I Chalker and Coddingtonl (|1988[ ) as an ef- 
fective description of the IQHE in a smooth random po- 
tential. In brief, the model is motivated in the following 
way. One considers electrons in a Landau level broad- 
ened by a potential with large correlation length. The 
electrons then drift along equipotential lines and tun- 
nel between the lines near saddle points of the random 
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FIG. 2 7 Chalker-Coddington ne t work model of the IQH tran- 
sition (|Chalker and Coddingtonl . ll988l ). The symbols ±c, ±s 
denote the components ±cos0, ±sm9 of scattering matrices 
at the nodes of the network. 



potential. When the energy is sufficiently close to the 
band center (classical percolation threshold) , the tunnel- 
ing probability becomes ~ 1, and a random network with 
directed links is formed. At each node of the network 
two incoming and two ongoing links meet. In the CCN 
model, this geometrically random structure is replaced 
by a regular square network, as shown in Fig. 1271 the 
disorder is accounted for by random phases associated 
with all links. A state ^ of the network is defined by 
i ts amplitudes on the edg es of the network. Originally 
(jChalker and Coddingtonl [l988) the network was char- 
acterized by a transfer matrix, as appropriate for finite- 
size scaling analysis of the localization length. Later 
( [Klesse a nd Metzlcr, 1995) an equivalent description in 
terms of a scattering matrix was introduced. Each real- 
ization of the network is determined by a unitary oper- 
ator U = UnUe acting on states ^ and modeUing the 
evolution of the state in a time step. Here Wat is an oper- 
ator describing the unitary scattering at nodes with am- 
plitudes ± cos6', ± sin0, as shown in Fig.[271 The second 
factor, lAE^ is a diagonal operator with random elements 



e^'*'^ on all edges e of the network. In the simplest for- 
mulation of the model, the angle 9 is the same for all 
nodes, and the phases cj),, are independent random vari- 
ables distributed uniformly over [0; 27r]. Changing the 
parameter 6 allows one to drive the system through the 
IQH transition, with the critical point at cos^ = 1/2. 

The CCN model has been extensively used for numeri- 
cal simulations of the IQH transition point; it turned out 
to be particularly well suited for the analysis of statisti- 
cal properties of energy levels and wave functions. Key 
results of computer simulations are reviewed below. 

The model has been generalized to other symmetry 
classes. This is most naturally done for the counter- 
parts of the IQHE in the superconducting classes C 
and D, namely, SQHE and TQHE. In these classes the 



symmetric spaces of the Hamiltonian (Table |T| are the 
groups Sp(N) and 0(N), and the required modification 
amounts to a replacement of the factors e*"^" G U{\) 
by the elements of Sp(2)=SU(2) for the SQHE (in this 
case the amplitudes are spin doublets) and of 0(1) for 
the TQHE, see Sec. IVLDl IVLEl for more detail. For 
several other symmetry classes non-directed generaliza- 
tions of the CCN have been constructed and used to 
study the corresponding critical behavior; specifically, 
this has been done for the ch iral classes (AIII, BDI, CII) 
(|Bocauet and Chalkeiil2003D. Sec. lVLFl and for the sym - 
plectic class AH (|Merkt et a/.l . [19981 : lObuse et a/] . l2007al ) 
considered in Sec. IVI.BI 

A further important aspect of the CCN model and 
its generalizations is that they can serve as a start- 
ing point for analytical work. We have already men- 
tioned established equivalences between the CCN and 
other IQH models (Pruisken model and superspin chain) 
in Sec. IVI.C.2I Further, a connection with the mod- 
els of disordered Dirac fermions has been established 
(|Ho and Chalkerj . I1996D . The network model of the 
SQHE has led to a number of exact analytical results. 

Sec. Eld] 



5. Localization length exponent 

Several microscopic models have been used to study 
numerically the critical properties at the IQH transi- 
tion. This includes tight-binding models. Landau-space 
models where the problem is projected on one or sev- 
eral Landau levels, and the CCN models. Sec. IVI.C.4I 
For a more detailed review o f the models the reader is 
referred to iHuckesteinI ( 19951 ). The first high-precision 
determination of the localization length exponent v has 
been a chieved in t he framework of the Landau-space 
model ( Huckesteirj . 19951 : [Huckestein and Krameil Il990l : 
iHuckestein et aLl . Il992l) ~with the result v = 2.35 ± 0.03. 
Results of later simulations on different models are all 
in agreement with this value, thus favoring the univer- 
sality of the IQH critical behavior. At present, the pre- 
cise value of the leading irrelevant scaling index, y, is 
known with much less accuracy. While several authors 
find values close t o y=0.4, e. g. (Ever s and Br enig, 19981 : 
lHuckesteinlll995D. in some ca s es val ues up to y=0.6 have 
been reported ( Kramer et a/.l . [20051 ) . 



6. Critical conductivity and conductance distribution 

In a number of works, the critical conductivity was 
found numerically in the range 0.5 - 0.6. More specif- 
ically, the results are: from 0.50 ifc 0.03 to 0.55 ± 
0.05 for diffe rent types of disorde r ( Huo et oil . Il993l ): 
0.50 ± 0.02 dC ammel and Bremgl, Il994[ ): 0.58 ± 0.03 
( Schweitzer and Markos, 20(5^^ Results in higher Lan- 
day levels are consistent with these values to the extent 
that the critical regime (of system sizes) could be reached 
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FIG. 28 Distribution function of the conductance at the IQH 
(solid) and the 2d symplectic (dashed) critical point using 

eriodic boundary co nditions and a square sample geometry 

Ohtsuki et"ttHl2004D 



(|Gammel and Everd . fl998l ) . 

Several authors also studied the conductance distri- 
bution 'P(g) of a square sample with periodic boundary 
conditions in the transverse direction; see Sec. IVI.B.3I 
for a quaUtative discussion of Vi^g) at criticality. As ex- 
pected, a scale invariant distribution was found, see Fig. 
E8l with the average (g) = .58 ib 0.03 and th e vari- 



ance var(5) = 0.081 ± 0.005 (IWang et alV Il996l); sim - 
il ar results were obta. i ned in ICho and Fished (llQQTbl ). 
in lKramer et (|2005l ): IOhtsuki et al\ (|2004l ) where the 
average found to be (q ) = .57 ± 0.02, as well as in 
ISchweitzer and MarkosI (|2005t ): the latter work yields 
(g) = 0.60 ±0.02. 



7. Wave function multifractality 

A high-precision evaluation of the multif ractal spec- 
trum at the IQH transition was carried out in lEvers tt al\ 
|2001) for the CCN model of a size L x L with L 
ranging from 16 to 1280. Figure [29] shows results for 
the /(a) spectrum. It is seen that after extrapolation 
to the thermodynamic limit, L — > oo, the /(a) spec- 
trum is well described by the parabolic form (|6.10p with 
Qfo — 2 = 0.262 ± 0.003. One observes that deviations 
from parabolicity - if they exist - are too small to be 
resolved in this plot. In Fig. [501 reduced anomalous di- 
mensions Ag/(7(1 — q) are plotted. This quantity is con- 
stant (equal to ao — 2) for an exactly parabolic spectrum, 
Eq. (|6.9p . The observed deviations from the constant are 
very small (~ 1%) and within the error bars. Since there 
is no small parameter in the problem, it seems that an 
accidental closeness of the spectrum to a parabola with 
such a high accuracy is very improbable. So, the nu- 
merical results are in favor of exact parabolicity of the 
multifractal spectrum, thus supporting the possibility of 



FIG. 29 Multifractal spectrum at the IQH transition. Fat 
solid line: numerical results for j(a) obtained from scaling of 
average IPR after extrapolation to L ^ oo. Datas for sev- 
eral finite values of L = 16, 128, 1024 are also shown. Dashed 
line: parabolic approximation, Eq. (|6.9|) . Inset: data points 
from typical IPR from systems with sizes L = 16, 128, 1024 
(open symbols) and extrapolation to infinite system size (full 
circles); data from average IPR are shown by solid line. 
([Evers et ali |2001D 



a conformal theory of the IQH critical point of the type 
discussed in Sec. IVLC.31 

It is worth mentioning that earlier studies of the 
IQH multifractalit y [the references can be found in 
([Evers et al\ . [200l[ )] gave considerably different results, 
showing, in particul ar, strong deviation s from parabolic- 
ity. As was shown in [Evers et al\ ([200l[ ). earlier numerics 
suffered strongly from the absence of ensemble averaging 
(its role was explained in Sec. lII.(J?5[) and from finite-size 
effects. Importance of a careful analysis of the latter is 
illustrated in Fig. [2S1 where also data for finite sizes of a 
system are included. 

A high-precision evaluation of multifractality allows 
to test the conformal invariance of the problem via 
Eq. (|6.5p . The pa rameter Ac = jmlM was found to be 
A^ = 1.22 ± 0.01 ([Evers and Brenisi [l998t [Evers et ali 



[200l[ ). implying, in combination with the above value 
of ao — 2, that (|6.5p is perfectly fulfilled and thus con- 
firming the expectation that the IQH critical theory is 
conformally invariant. 



8. Statistics of the two-point conductance 

Closely related to wavefunction multifractality is the 
statistics of two-point con ductances, Sec. III.C.9I In 
[Klesse and Zirnbaued ( 200l[ ) a relation between statisti- 
cal properties of the wave functions and two-point con- 
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FIG. 30 Anomalous multifractal dimensions Aq [divided by 
q{l — g)] at tlie IQHE critical point, as obtained from ex- 
trapolation of the average inverse participation ratio (IPR). 
In the case of exact parabolicity, the plotted quantity should 
be constant. The dotted lines indicate the error bars ob- 
tained for Qo — 2. The circles give the A*^P/q(l — q) as 
obtained from the typical IPR. As explained in Sec. III. C. 51 
Ag^'P = A, for q < q+; for IQHE with parabolic spectrum 



(tO)) . ((6T0)| we find q+ = [2/{ao - 2)Y^'^ ~ 2.76. Thus, the 
last two data points for A^^^ /q{l — q) are in the range q > q+, 
which explains the ir downward deviations from Ag/g(l — q). 
(|Evers et all l200ll ) 



ductances was derived, 



'2T^p{ymf{yrn/yi)) = {F{Ti,, 

F{T)= r '^f(T-^\l-e'^VT^f 
Jo Svr V 



(6.11) 



Here f{x) is an arbitrary function, yi = t/m — 

l^'mp are wave function intensities at two links I and 
m for an eigenstate \1/ of a closed network, and Ti^m is 
the two-point conductance defined for an open network 
with the edges I and m cut and attached t o two ter- 
minals. This result was used in lEvers et all ( 200ll ) to 
derive the the relation (|2.63p between the correspond- 
ing critical exponents. The parabolic spectrum for wave 
function multifractality, Eq. (|6.9p . found for IQHE trans- 
lates thus into parabolic spectrum of exponents Xq for 
the two-point conductance, 



X. 



q<l/2 



Xtq{l - q) ; 



X, 



q>l/2 



Xt/A , (6.12) 



with Xt = 2(ao - 2) = 0.524 ± 0.006. In an earlier 
work ( JanBen et a/.l . [l999f ) an explicit expression for the 
distribution of the two-point conductance on the CCN 
was derived under the assumption of the parabolic law 
(pT^ for Xq. In Klessc and Zirnbaucr (2001) this result 
was used to test the conformal invariance of the theory, 
utilizing a numerical analysis of the moments of the two- 
point conduc tance in the quasi-l D (cylinder) geometry. 
The result of lJanfien efall (|l999l for the distribution of 
the conductance T between the points (0, 0) and {x, y) in 



V{T) = 



2^-l/2^-X,/4 



T2(XanC)3/2 Vcosh^t-T-i 
C = (VF/7ra)| sinh[7r(a; + iy)/W]\ , (6.13) 

where a is the non-universal microscopic scale that sets 
the length unit. The numerically determined moments of 
T were in perfect agreement with this formula, thus sup- 
porting the conformal invariance. The best fits yielded 
the values Xt — 0.54 ± 0.01 from the analysis of the mo- 
ments (T"+i/2) and Xt = 0.57 ± 0.05 from the analysis 
of the typical conductance, (InT). These values are con- 
sistent with the above result Xt — 0.524 ± 0.006. 



9. Classical percolation vs. quantum Hall effect 

If the disorder correlation length is large, there is 
an intermediate, parametrically broad, range of energies 
where the physics is dominated by classical percolation. 
On corresponding length scales, tunneling between the 
percolating contours is ex ponentially small. This range 
of energies was studied in iMil'nikov and Sokolovl ( 19881 ) 



where the scaling of the localization length with an ex- 
ponent i> — Vperc + 1 — 7/3 was found (I'porc — 4/3 is the 
correlation length index of the 2D percolation problem). 
In this regime quantum interference effects play no role. 
When energy approaches still closer the critical point, the 
probability of tunneling between the contours ceases to 
be small and the quantum interference becomes crucially 
important - the system enters the true critical regime of 
the quantum Hall transition. It is this latter regime that 
is described by the CCN model. While the "quasiclassi- 
cal QHE" exponent = 7/3 is remarkably close to the 
numerical value of the true QH exponent v = 2.35 ±0.03, 
this coincidence is apparently fully accidental, as the 
physics of the true QH crit ical regime and the intermedi - 
ate quasiclassical regime of lMiPnikov and Sokolovl (I1988D 
is completely different. 

It is worth mentioning that the physics of the in- 
termediate, quasiclassical regime, where the physics is 



(lEvers and Brenid. 1994, 1 


19981: ICammel and Brenis 


1996 


; Klesse and Metzler. 


Il995t Kratzer and Brenig 



is interesting in its 
own right. In particular, quasiclassical time evolution 
generates a long time tail in the velocity correlation func- 
tion {vx{t)vx{0)) ^ i"^, which leaves crxx(^) with a non- 
analytic cj-dependence inside an intermediate (quasiclas- 
sical) frequency window. The corresponding results may 
be relevant to experiments if the latter are performed on 
structures with smooth disorder, in which case the true 
QH criticality may in fact be totally unobservable for 
realistic temperatures. 
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10. Experiment vs. theory. Interaction effects 

a. Experimental results. Tlic IQH transition and the as- 
sociated critical properties have been studied in numer- 
ous experiments. All the basic features — the existence 



of phase transitions with critical values axy 



1/2 



and axx of the order of unity, as well as the power-law 
scaling behavior — are in agreement with the theoreti- 
cal expectations. The situation with the values of critical 
exponents is not so simple, as we are going to discuss. To 
do this, we will have to touch the question of interaction 
effects, which is left out in the rest of this review, except 
for Sec. IVII.AI 

Experiments yield the following results for the critical 
exponents. First, the index v o f the localization length 
is found to he V = 2.3 ± 0.1 llKoch et all Il99lh: this 
value w as confirmed more recentlv in lHohls et al\ (|200lL 
l2002bD . Second, the width of the critical region (peak in 
axx and plateau transition in cixy) scales wit h the tem- 
perat ure T as AS cx T« with k = 0.42 ± 0.04 ()Wei et all 
While different values of k were obtained in other 
works, it was emphasized in Ivan Schaiik et all (l2000D 
that this results from macroscopic inhomogeneities that 
complicate observation of the true IQH critical behav- 
i or with k 0-42. More recent work of the same group 
( Pruisken et al liooiidi; Visser et ^,12006') favors again 
K=0.56 ± 0.02, however. On the other hand, th e impact 
of den sity inhomogeneities was reconsidered bv lLi et al\ 
(|2005D . [t was found there, that for short-range dis- 
order, when the true IQH criticality can be achieved, 
K = 0.42 ±0.01, whereas the larger value k=0.58 was as- 
cribed to impurity clustering. Finally, the frequency scal- 
ing of the transition w idth was found to be AB ^ uj'^ , 
with C = 0.41 ± 0.04 (lEngel all Il993l ). A more re- 
cent work ( Hohls et all 2002a ) yields a result consistent 
with this value, but with somewhat larger uncertainty, 
(^ = 0.5±0.1. To summarize, the experiments yield v 
that agrees with its numerical value (2.35 ± 0.03), as 
well as the dynamical exponents zt = l/nv — 1 and 
z = l/C,!^ ~ 1. The remarkable agreement in the value 
of v is in fact surprising, in view of the electron-electron 
interaction. 



b. Finite-range interaction. We consider fi rst the case of a 
finite- r ange interact i on ■^(r — r' ) , following iLee and Wand 
(ll996D : IWang" et al. I (|2000D . ;n this case the interaction is 
irrelevant. Indeed consider the Hartree-Fock interaction 
between the (close in energy) states a and /3, normalized 
to the level spacing. 



'M2 



\ = pL' j dWr'([|^„(r)|^|VMr')l 

-^„(r)C(r')V'/3(r')V'a(r')]«(r - r')). (6.14) 

In the (T-model language, see Sec. IILB.II the scaling of 
(I6.14P with the system size L is governed by the scaling 
dimension of the operator 

Qn (r)Q^2'(r) - Q\'^{v)QI\{v) . (6.15) 



While each of the two terms in (|6.15p is relevant in the RG 
sense, having a dominant negative scaling dimension A2 
(which governs the IPR scaling, Sec JII.C.llVI.C.Tp . the 
difference (pTTSl) is RG-irrelevant (fWegneil . 119801) . with 
a scaling dimension X2 > 0. The numerical value of 
X2 at the I QH critical point wa s estimated to be X2 = 
0.66 ±0.04 ("Le e and Waji3 . ll996t) . With increasing L the 
interaction (|6.14p scales as A cx so that the fixed 

point is unaffected by it. This implies that the critical 
index v of the localization length and the multifractality 
spectrum Ag remain the same as in the non-interacting 
problem. This is also true for the dynamical exponent 
z governing the destruction of localization by finite fre- 
quency: the corresponding scaling variable is with 
z = 2. (In general, for a non-interacting transition with 
finite DOS in d dimensions the scaling variable is uj/q'^ 
( Wegneii Il976t) ). The interaction cannot be fully dis- 
carded, however (it is said to be dangerously irrelevant), 
as the conductivity at finite T would be zero without the 
interaction-induced dephasing. The dephasing rate scale s 



T, CX with p 



2x2/ z ()Wang et all l2000t) . 



1/2 



Thus, the dephasing length is oc 
zt = 2/p — 2z/{2x2 + z), yielding zt 



1.2 with the 

above estimate for X2- Therefore, for a system with a 
metallic gate (which screens the interaction) one expects 
the exponents v ~ 2.35, z — 2, and zt — 1.2. 



c. Coulomb interaction. The situation with l/r Coulomb 
interaction (as in typical experiments) is much less clear. 
In this case the interaction is R G-relevant and drives 
the system to a nove l fixed point (jBaranov et al\ . 120021 : 
iLee and Wand . Il996[ ). While the conductivity and the 
screen ed compressibility c^n/9u rema i n finite at the tran - 
sition (jBelitz and Kirkpatriclj . 11994 iFinkelsteinl Il990[ ). 
much less is known theoretic ally about other critical 
prop e rties. Several authors ( Polvakov and Samokhinl . 
119981: IWang and Xio^ . l2002h argued that charging ef- 
fects analogous to those responsible for the linear 
Coulomb gap in the tunneling DOS of the insulator will 
lead to z — Zt — \ (which is the natural scaling dimen- 
sion of the 1 /r interaction) . The status of this argumen- 
tation is unclear, however, for the following reasons: 

(i) These arguments essentially identify z and 
Zt with a dynamical expo nent (23 in notations of 
iBelitz and Kirkpatriclj ( 1994t) ) governing the scaling of 
the density response function. It is known, however, 
that in metals the plasmon pole governed by this expo- 
nent, Z3 = 1, determines the interaction- induced quan- 
tum correction to tunneling D OS but not to conductivity 
( Altshuler and Aronovl Il984} ) , in view of gauge invari- 



ance ( Finkelsteinl . 119941 ) . The conductivity correction is 



governed by the conventional diffusion pole in the irre- 
ducible density-density response function, which corre- 
sponds to the dynamical exponent z"^ = 2. 

(ii) Each Goldstone mode (or, equivalently, conserved 
quantity) is in general characterized by a dynamical ex- 
ponent. For the QH transition this implies that, in 
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addition to the exponents related to the particle num- 
ber conservation - = 1 and its irreducible counter- 
part, z^" = 2, - there is another exponent associated 
with the energy conservation. This latter exponent, 
controlling the renormalization of the frequency term 
in th e a-model, is denoted by C in IFinkelsteinI (1993), 



Zi m 



Belitz and Kirk patricd (|1994[ ). and by 2 



+ 7* 



iBaranov al. (2002). It is believed to govern the fre- 
quency scaling of the conductivity at the critical point of 
the Anderson transition in a system with Coulomb inter- 
action and broken spin-rotation invariance in 2-1- e dimen- 
sions (jBelitz and Kirkpatricd Il994l; IFinkelsteinI . Il990[ ). 
One may thus expect that this dynamical exponent plays 
a central role at the quantum Hall tr ansition as well 
()Baranov et all l2002t iBurmistrovl . l2006t) . 

The problem of the index of the localization length 
is also far from being solved. While it was found that ly 
is equal to it s non-intera. c ting v alue within the Hartree- 
Fock theory (jYang et all Il995l ). it is not clear whether 
this should be applicable to the true fixed point in the 
problem with Coulomb interaction. To summarize, in 
our view, the theoretical problem of the critical behavior 
in the presence of Coulomb interaction remains open. In 
particular, it remains to be seen whether the remarkable 
agreement of v with its non-interacting value as well as 
z = zt = i — as suggested by experiments — are indeed 
exact properties of the interacting problem. 



In the IQHE the step Aa^y between the quantized values 
of the Hall conductivity is e'^/h per spin orientation. In 
the case of the S QHE, the elem e ntary charge e is replaced 
by h/2, yielding (jSenthil et ail . Il999l ) 

Aaly = 2n ■ h/Sir , (6.17) 

where the factor 2 accounts for the spin and n for the 
valley degeneracy. In particular, for the case of Dirac 
fermions in a d 4- id superconductor (n = 2), one finds 
two SQH phases with quantized values 



±ft/47r 



(6.18) 



In the presence of disorder, these two phases become An- 
derson insulators separated by the SQH transition. The 
field theory of this problem is analogous to Pruisken's 
theory of the IQHE (Sec. IVI.C.ll) — it is a cr-model of 
the class C with the topological t erm (16.11). pos s essing 
a critical point at 9 = (2n -|- l)7r dA ltland et al!, l2002t 



iRead and GreCTl l2000l : FSenthil et aq ,1998., 1999. ). The 



corresponding flow diagram is expected to have qualita- 
tively the same form as for the IQHE, Fig. |26l The criti- 
cal behavior at the SQH transition is analyzed below. In 
the presence of a Zeeman term, the spin-rotation sym- 
metry is broken, and the system crosses over to the sym- 
metry class A of the IQHE transition (ICui et al , |2004 
iKagalovskv et al\ . 119991 : ISenthil et aLl . ll998lll999D . 



2. Mapping to percolation 



D. Spin Quantum Hall Effect (Class C) 

1. Physical realization 

The SQHE is a counterpart of the IQHE in 
superconductors with broken t ime-reversal bu t pre- 
served spin-rotat i on invariance (ICruzberg et all Il999t 
IKagalovskv et all \l99^ . ISenthil et aLl . ll999l). The class- 
C Hamiltonian satisfies the symmetry (|4.12l) (with Ty 
the Pauli matrix in the particle-hole space) and has the 
block structure (|4.1ip . Several possible physical real- 
ization of SQHE systems have been proposed: (i) a d- 
wave superconductor with complex d^2 ,,2 + id xv pairing 
(IKagalovskv et all Il999t ISenthil et all . 1l999D that was 
conjectur e d for high- superconductors ( Balatskvl . ll998l : 
iLaughlirJ . Il998f ): (ii) granular supercond ucting film in 
a magnetic field (iKagalovskv et all Il999f) : (iii) a state 
of composite fe rmions at filling f ractio n u = 5/2 with 
d-wave pairing (iRead and Greenl [20001 ) . as proposed in 
iHaldane and Rezavil 1^198^. 

Similar to the case of the IQHE, the key signature of 
the SQHE is the quantization of the appropriate Hall con- 
ductance. Specifically, while the quasiparticle number is 
not conserved, the spin is, so that the relevant quantity 
is the spin Hall conductivity a^y. It describes the trans- 
verse spin current induced in response to a gradient of 
the Zeeman field. 



Th e network model for the SQHE ( Kagalovskv et all 
I1999I ) is the SU(2) versi on of th e Chalker-Coddington 
IQHE network, see Sec. IVI.C.4I and Fig. ^ The di- 
rected links of the network carry doublets of complex 
fluxes representing propagation of spin-1/2 particle. The 
scattering at each node is spin-independent and defined 
in the same way as for the IQHE. Each realization of 
the network is characterized by a set of random 2x2 
spin matrices Ue associated with all edges e of the net- 
work. In view of (|4.12p . the evolution operator of the 
network U satisfies the symmetry U = TyU*Ty, implying 
that Ue e SU(2). Starting with this network model, it 
turns out to be possible to establish a remarkable prop- 
erty of the SQH transition: some physical observables 
and critical indices can be calculated exactly via map- 
ping to classical percol ation. This was shown f or the DOS 
and the conductance in iGruzberg efai] (|l999l ) via the su- 
persymmetry; an alternative der i vati on of these results 
was p resented in lBeamond et al. I (12002) . In Mirlin et al\ 
(I2OO3D the mapping was extended to all two- and three- 
point correlation functions describing, in particular, the 
wave function statistics. It was also shown there that 
the mapping breaks down for generic rt-point correlation 
functions with n > 3. 

We briefly sketch the idea of the approach of 
iBeamond 6^0/] (|2002l ): [Mirlin a/.l(|2003D . The primary 
objects are Green functions on the network. 



3. 



-dB^{y)/dy] 



(6.16) 



G(e',e;z) = {e'\{l ~ zU)-^\e) 



(6.19) 
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i) 



ii) 

s c 



iii) 

s c 



FIG. 31 Possible configurations of paths passing four times 
through a network node. The symbols c and ±s denote the 
elements cos 6, ± sin 9 of the S-matrix at the node. The state- 
ment 2 in Sec. IVI.D.2I allows one to get rid of the quantum 
interference contribution (iii) and to associate with the contri- 
butions (i) and (ii) the weights cos'^ and sin^ 9, respectively, 
yielding a mapping to the classical percolation. 



for z = expi(e±H7) they have a meaning of retarded {Gr, 
\z\ < 1) and advanced {Ga, \z\ > 1) Green functions at 
energy e. The Green function is straightforwardly repre- 
sented in the form of a sum over paths 



G(e, e'; z) 



E 

paths e'- 



■zUe^Sj-zUe^^j^Sj+i- 



(6.20) 



where Sj is the corresponding matrix element (cos 9, sin 0, 
or — sinO) of the S'-matrix between the edges Cj and ej+i. 
Equation (|6.20p generates a convergent expansion in pow- 
ers of z when \z\ < 1; otherwise the identity 



Gt(e,e';z)-l-,5ee'-G(e',e; (z*)-^) 



(6.21) 



is to be used. As shown below, each of the sums over 
paths obtained by substituting (|6.20p , (|6.2ip in products 
of n < 3 Green functions can be reduced after disorder 
averaging to a sum over classical paths (hulls) in the per- 
colation problem. This remarkable reduction crucially 
relies on the following two statements: 

1. Only paths visiting each edge of the network either 
or 2 times are to be taken into account; contributions 
of all the remaining paths sum up to zero. 

2. Using the statement 1, it is easy to see that each 
node may be visited 0, 2, or 4 times. The second state- 
ment concerns the nodes visited four times. As illustrated 
in Fig. [nil there are three possibilities how this may 
happen; the corresponding contributions have weights (i) 
cos^ 9, (ii) sin^ 9, and (iii) — sin^ 9 cos^ 9 from the scatter- 
ing matrix at this node. The statement is that one can 
equivalently take into account only the contributions (i) 
and (ii) with the weights cos^ 9 and sin^ 9, respectively. 

After application of the statement 2 to all nodes, the 
network is reduced to a weighted sum over all its possible 
decompositions in a set of closed loops (such that each 
edge belongs to exactly one loop). These loops can be 
viewed as hulls of the bond percolation problem. As a 
result, the correlation functions (averaged products of 
Green functions) are expressed in terms of classical sums 
over the percolation hulls. The results obtained in this 
way are listed below. 



3. Density of states and localization length 

The result fo r the average of a sin gle Green function 
reads (|z| < 1) (|Beamond et aLl . l2002[ ): 



(TrG(e,e,0)) = 2-^P(e;iV) 

N>0 

(TrG(e,e,z-i)) = ^P(e;7V)2 



^2N 



2N 



(6.22) 
(6.23) 



N>0 



where P(e, N) is the probability that the edge e belongs 
to a loop of the length TV. (In the bulk of a large system, 
L — > oo, or for periodic boundary conditions, this prob- 
ability does not depend on e, P{e,N) = P{N).) This 
yields the DOS 



piE) = (1/2^)[1- ^ P{N)cos{2NE)] 



(6.24) 



N>0 



In the insulating phases {t = cos^ 9^1/2) this yields 

p{E)-^Ti~^{N'^)E'^, (6.25) 

which is the expected behavior of DOS in the class C. On 
approaching the percolation transition point, t = 1/2, the 
characteristic diameter of largest loops diverge. 



4/3. 



e^l^-l/2r^ 

At the critical point (t = 1/2) 

P{N) - iV-2M. = 



(6.26) 



(6.27) 



where dh =7/A is the fractal dimension of the percola - 



p ercola - 

tion hull (llsichenkol. [l992HSaleur and Duplantierl. Il9 87^. 
yielding (jBeamond et all l2002t iGruzberg et all Il999ll 



p{E) ^ \E 



1/7 



(6.28) 



The characteristic length of loops contributing to (|6.28p is 
Ne ^ E~^, yielding their characteristic "diameter" ^ 
l^j-i/dh _ \E\~'^^'^ , which is the locahzation length of 
the states with energy E. The percolation hull exponent 
dh plays therefore a role of the dynamical exponent for 
the SQH transition. 



4. Conductance 

To define the two-terminal conductance g, one opens 
the system by cutting two subsets of links and attach- 
ing them to two reservoirs. The average dimensionless 
conductance is given by 



id) 



(6.29) 



eGlout; e'e2i, 



where P{e, e') is the probability that a path incident from 
the second reservoir on the link e' escapes to the first 
contact via the link e, and the sum goes over all such 
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links. Equation (|6.29p was used in (|Cardvl . [2000h to cal- 
culate the critical conductivity by determining the con- 
ductance of a wide sample {W > L), {g)L/W — V3/2. 
A very close result was earlier obtained by numerical 
evaluation of the Kubo formula, ac — 2 x (0.45 ± 0.01) 
( Eversl IT997I ). It is worth stressing that, despite of the 
vanishing density of states, the critical conductance is 
finite. From the point of view of the Einstein rela- 
tion, a = hpD, this results fro m a mutual cance l lation 
of the percolat ion exponents (|Evers and Breni j . 1 1994 
IZiff et al\ . \im\) . p{E) - ~ £;2/dh-i = £;i/7 ^nd 

For the averag e two-point conductan ce (Sec. III.C.9|) . 
Eq. (jOg)) yields (|Gruzberg et all |1999| ) 



(5(e,e')) =2P(e,e')~r-i/2 , 
where r ^ 1 is the distance between e and e'. 



(6.30) 



5. Higher correlation functions and multifractality 

The results present ed in this subsection were obtained 
in lMirlin et al ] (|2003). To determine the fractal dimen- 
sion A2 governing the scaling of two-point correlations of 
wave functions (|2.30p . (|2.32p . one considers the correla- 
tion functions 

V{e',e;j) = (2^)-2(Tr[G(e', e; z) - G(e', e; z'^)] 

X [G(e,e';z)-G(e,e';z~i)]), (6.31) 

D(e',e;7) = (2^)-2(Tr[G(e, e; z) - G(e, e; z"!)] 

X Tr[G(e',e';z)-G(e',e';z-i)]), (6.32) 

with a real z — e ' < 1 and 7 <C 1 playing a role of the 
level broadening. The mapping to percolation yields for 
the averaged products of two Green functions entering 

(IS3I1), 

(TrG(e',e;z)G(e,e';z)) = (TrG(e', e; z-i)G(e, e'; z"!)) 
= -2^P(e',e;iV)z2^,(6.33) 



N 



(TrG(e',e;z)G(e,e';z-i)) 



-2 5]Pi(e',e;iV)z2^,(6.34) 



N 



where P{e',e;N) and Pi{e',e;N) are probabilities that 
the edges e and e' belong to the same loop of the length 
N (resp. with the length TV of the part corresponding 
to the motion from e to e'). According to the classical 
percolation theory, P and Pi scale as 

P(e', e, N), Pi(e', e, N) - A^-s/V-i/-^ , r < N^''^ 

(6.35) 

and fall off exponentially fast at r ^ N^^'^ , where r is the 
distance between e and e'. This yields for the correlation 
functions in (|6.33p and (|6.34[) (which we abbreviate as 

(GflGfl), {GaGa), (GrGa)) 



(GrGr) = (GaGa) ^ (GrGa) - r-i/^ 



-4/7 



(6.36) 



i n full agreement w ith the scaling argument of 
( Gruzberg et all Il999t ). However, these leading order 
terms cancel in (j6.3ip . and the result is non-zero due 
to the factors z^^ only, implying that relevant N are 
N ^ 7"-^. As a consequence, {{Gr — Ga){Gr — Ga)) 
scales differently compared to (j6.36p . 



V{e', e; 7) = ^ Y.'^Pir, N) ~ Pi{r, N)]{1 



-2Ny 



N 



<C7- (6-37) 



The analysis of the correlation function (|6.32p yields sim- 
ilar results. One thus finds the scaling of two-point wave 
function correlations for r < ^e, 

i'(V'*Je)^,o(e)V.Me')V'*^(e')> 1 ( ^A'^ 33. 

with a, P labeling the spin indices. This implies that the 
fractal exponent A2 = —1] is 



-1/4, 



(6.39) 



at variance with what one might naively expect from 
the r^^/^ scaling of the diffusion propagator (GrGa), 
Eq. (|6.36p . An analogous calculation for three-point cor- 
relation functions yields A3 = —3/4. For correlation 
functions of higher orders (determining, in particular, the 
exponents Aq with q > 3) the mapping to percolation 
breaks down. 

The point q = 3 deserves special attention. It satisfies 
the relation 



T{q) = qxp + Ag = , 



(6.40) 



where Xp — 1/4 is the scaling dimension of DOS defined 

by p{E) ^ C_E^''- It separates two regimes with differ- 
ent scaling of correlation functions. For smaller q, when 
r(q) > 0, the correlation functions 



n(«)..,^(ei,...,e,;i?i, 
= (TrG,, (ei, 62; e'^^) ■ ■ • G,, (e„ ei; e*^')),(6.41) 
where Sj = R or A, show the scaling 

ni«) (ei,...,e,;i?i,...,£;,)^r-«-'' (6.42) 



and are, to the leading approximation, independent of 
the indices Si. However, when one calculates the wave- 
function correlations, 

P('')(ei,...,e,;£;i,...,i?,)-(27r)-'' 

X {Tt[{Gr - GA)(ei, 62; e^^i)(Gfl - GA)(e2, 63; e*^^) 



X...X (Gfl,-GA)(e„ei;e*^0]>, 



(6.43) 



these leading-order terms cancel, yielding the multifrac- 
tal behavior. 



I?(9'(ei,...,e,;i;i,, 



E ' 



(6.44) 
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On the other hand, for q > 3, when r(q) is negative, 
the correlation functions 111'^,...,^, start to depend in a 
singular way on the infrared cutoff {^e) and scale in the 
same way as P*^'' , Eq. (|6.44p (with a numerical prefactor 
depending on indices Si), similarly to the conventional 
Anderson localization transition. 

It is instructive to analyze this situation within the 
field-theor etical approach to the wave - function multi- 
fractalitv (iBernard and LeClaiii 2002b Bhaseen et al. , 



20001: iDup lantier and Ludwid. Il99ll: iMudrv et all 11996 : 
Wegneij , Tl 980 , 1985). In the renormalization-group lan- 
guage, r(g) defined by Eq. (|6.40p are scaling dimensions 
of operators of the type C^'^ ~ V'si^]' ■ • • V's,''/'!' i where 

tp, ■0'!' are electronic fields. Averaged products of Green 
functions are expressed as correlation functions of the 
corresponding operators O^''''; in particular, (|6.4ip takes 
the form 

ni^..., - (TrOW^(e2)0«3(e3)...OW,(ei)). (6.45) 

To calculate the scaling behavior of such correlation 
functions, one applies the operator product expan- 
sion (OPE) toupl antier and Ludwi3 , ll99ll : lMudrv et all . 
[I996- Wcgncr, .198^ Generically, the identity operator 
will be among those generated by the OPE. Moreover, 
under the condition T{q) > it will be the most relevant 
operator and will dominate the expansion, leading to the 
gap scaling 11*^'^) ~ r^i'^W ^ in agreement with (I6.42|) . On 
the other hand, if T{q) < 0, the operator O'^'^-' will give a 
dominant contribution to OPE, leading to a multifractal 



type of scahng, 11^'^ 



oc r 



-gr(l) 



(r/^E)i^(9), asinEq. (lOl) . 



What is, however, non-trivial from this point of view, is 
that the scaling of the wave function correlator (|6.43p has 
the multifractal form (|6.44p independently of the sign of 
r(g). This means that in the regime r{q) > the leading 
(gap scaling) terms (I6.42p cancel in the particular com- 
bination of the functions n*^') corresponding to I?^'^ , and 
subleading terms determine the result (|6.44p . 

A related analysis can also be performed for the mo- 
ments of the two-point conductance, Eq. (|2.62p . The 
corresponding exponents Xq are found to be linked to 
the wavefunction multifractal indices via (|2.64p . 



6. Numerical results 

The nume rical simulations of the SQHE network 
( Mirlin et all [2003 ) have allowed to confirm the analyt- 
ical predictions (Sec. IVI.D.3ffVI.D.5p as well as to deter- 
mine some physical quantities that are not known ana- 
lytically, - most notably, the whole spectrum of multi- 
fractality. We present a brief summary of the numerical 
results. 

In Fig. [321 the numerically calculated DOS p{E) for 
different system sizes L is displayed. After a proper 
rescaling all data collapse onto a single curve. The scale 
invariance of p{E) at criticality is reminiscent of the 
analogous property of the level statistics at the conven- 
tional Anderson or QH transition. Sec. III. El At £' 3> ^ 




FIG. 32 Scaling plot of the density of states for system sizes 
L — 16(o), 32(n), 96(o). Dashed and dotted lines indicate 
power laws (dashed: E^^"^ , dotted: E'^), S = l/27rL'^/'' denotes 
the level spacing at i? = 0. Inset: s ame data on a line ar scale 
and the RMT result (solid curve). (jEvers et all l2003l ) 



the critical DOS scales as in agreement with 

the analytical prediction (|6.28p . On the other hand, 
at E ^ d one observes an oscillatory structure quali- 
tatively analogous to the RM T behavior for the class C 
(|Altland and Zirnbaueii Il997l ). 

The anomalous multifractal dimensions Ag [divided by 
q{l — q)] are shown by a solid line in the upper panel of 
Fig. [331 They have been obtained from the scaling of 
the average IPRs. According to the analytical calcula- 
tions fSec. lVI.D.5|) , Ag/(?(1 - q) is equal to 1/8 for both 
q — 2 and q — 3; this value is marked by the dashed 
line in the figure. The numerical results agree perfectly 
well with the analytical findings at g = 2 and q = 3. 
Furthermore, the parabolic dependence may serve as a 
numerically good approximation, 



Aq~q{l — q)/8. (approximate!) 



(6.46) 



Nevertheless, Eq. (|6.46p is not exact: at < g < 2 
the numerically found Aq show clear deviations from ex- 
act parabolicity (|6.46p . which are of the order of 10% 
near q — 0. In particular, the deviation of the limit- 
ing value A,/g(l - q)\q^o = 0.137 ± 0.003 from 1/8 well 
exceeds the estimated numerical uncertainty. The lower 
panel of Fig. [331 depicts the singularity spectrum f{a). 
The dashed line represents the parabolic approximation 
Eq. ([07|1 with ao-2 = 1/8, corresponding to ([OS)) . 
The deviation of ao - 2 = 0.137 ± 0.003 from 1/8 high- 
lighted in the inset corresponds to non-parabolicity of Tq 
discussed above. The numerical results for the multifrac- 
tality at the SQH transition rule out conjectures of crit- 
i cal field theories that predict exac tly parabolic spectra 
([Bernard and LeClaid , [200l l2002bD . 

In Fig. [34l the scaling of the average {g) and the typ- 
ical (7typ = exp(ln(7) values of the two-point conduc- 
tance is shown, along with analogous quantities (|Gp) 
and |G|typ = exp(ln|Gp) for a closed system, |Gp = 
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FIG. 33 Upper panel: anomalous dimensions Aq (solid line) 
and Ag^P (circles) at the SQHE critical point. Numerical re- 
sults agree very well with analytical findings A2 = — 1/4 and 
A3 = —3/4. Data are presented in the form Aq/g(l — q) em- 
phasizing deviations from exact parabolicity (dashed line). 
Lower panel: singularity spectrum /(a) (solid) and the 
parabolic approxi mation (dashed). Inset: magnified view of 
the apex region. (jMirlin et all |2003| ) 
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FIG. 34 Scaling of the two-point conductance with distance 
r between the contacts: average value (empty symbols), (g), 
and typical value (filled symbols), gtyp = exp(ln g), in systems 
of sizes L — 128(n) and L — 196(o). Also shown is scaling of 
the two-point Green function, (|Gp) and |G|typ = exp(ln |Gp) 
(L = 128(A), L = 196(o)). The lines correspond to the r"^/^ 
(dotted) and r"'^^* (dashed) power laws. Deviations from 
power-law sc aling at large valu es of r are due to the finite 
system size. (jMirlin et all 120031 ') 



— TrG(e',e; 1)0(6,6'; 1). For the average values, {g) and 
(|Gp), the numerics fully confirm the theoretical results 
(I6.30p . (|6.34p predicting that both quantities scale as 
r^-^/'^ and, moreover, are equal to each other. A non- 
trivial character of the equality (g) = (|Gp) is well il- 
lustrated by the data for typical quantities: gtyp and 
|G|typ are not equal. Nevertheless, they are found to 
share a common scaling: gtyp, \G\typ ~ r~'^*, confirming 
the analytical expectations. Furthermore, the numeri- 
cally obtained value of the exponent, Xt ~ 3/4, is in 
agreement vifith the theoretical prediction based on the 
relation (PTM)) . Xt = 2xp + 2(ao - 2) ~ 0.774. 



E. Thermal quantum Hall effect (class D) 

1. Physical realizations and general considerations 

Systems belonging to the class D are disordered super- 
conductors where both time reversal and spin rotation 
symmetries are broken. The corresponding Hamiltonian 
has the structure described in Sec. lIV.Dl see Eqs. (|4.8p . 
(14. 9p and text below them. Possible physical realizations 
of this symmetry class include: (i) d-wave superconduc- 
tors with strong spin-orbit scattering, (ii) p-wave paired 
states of spinless or spin-polari zed fermions, e.g. paire d 
states of composite fermions ( Read and Greenl l200Cl() : 
(iii) triplet od d-parity (p- o r /-wa ve) superconductors, 
like SrRu304 (|Nelson et all I2004D : (iv) type-II super- 
conductors in a strong magnetic fi eld in the presence 
of spin-orbit scattering impurities ( Senthil and Fished . 



|200C|). While neither the quasiparticle number nor the 
spin are conserved for this symmetry class, one still can 
speak about thermal transport. The TQHE corresponds 
to the quantization of the ratio K^y/T of the thermal 
Hall conductance to the te mperature in units of Tr^k^/6h 
(|Senthil and Fished . l2000l) . 



For a combination of reasons, class-D systems show 
particularly rich behavior from the point of view of lo- 
calization, quantum phases, and phase transition. First 
of all, class D allows for two mechanisms of 2D critical- 
ity, Sec. IVI.Al (i) a topological 0-term associated with a 
quantum-Hall-type transition and (ii) a metallic phase, 
in view of broken spin-rotation invariance. Thus, gener- 
ically, three phases are possible: metal, insulator, and 
quantized Hall conductor. A further striking feature of 
class D is that the type of disorder affects crucially the 
phase diagram. At the level of the cr-model, the reason 
is believed to be that the relevant target space has two 
disconnected pieces, and that, depending on the choice of 
the underlying microscopic model, it may or may not be 
necessary to consider configurations containing domain 
walls on whic h the g-model f ield jumps between the two 
components ("Bocquet et al. Chalker et all 120021: 

[Gruzbcrg et al, 2005; Read and Ludwig, 20oi^In the 
following we mai nly concentrate o n the C ho-Fisher (CF) 
network model (|Cho and Fisheil . Il997al ) of the TQHE, 
which is generic in the sense that it displays all three 
possible phases. Other models of disorder will be briefly 
discussed in the end of Sec. IVLE.2I 
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Metal 



sin2(e) 

FI G. 35 Phase diaRram of the Cho-Fisher model as obtained 
m l|Chalker et ali |2002| ) from transfer matrix calculations. 
The plane is spanned by the parameters sin^ 9, the interpla- 
quette tunneling probability, and p, the concentration of vor- 
tex disorder: these control the short-distance values of the 
conductivity components a^y and (Txx respectively. 



2. Network model and phase diagram 

To obtain a disordered network model of class D, one 
can start from the ordered network, Fig. 1271 and then 
allow for independent fluctuations of the node parameters 
0i with some distribution function V{9). The CF model 
corresponds to the choice 



5+^o)+|<5(e+0o-7r), (6.47) 



implying that disorder is introduced as isolated defects 
by making the change 6*^— 6'or6'^7r — 6', fora sub- 
set of nodes randomly distributed with a concentration 
p. This amounts to flipping signs of either both sinfl or 
both cos 9 associated with such a node. This procedure 
can be viewed as the insertion of two additional half-flux 
lines into two plaquettes adjacent to the node and be- 
longing to the same sublattice. Note that the vortex pair 
appears with equal probability on the C- or S-sublattice. 
It is this feature that distinguishes the CF model from 
the rand om bond Ising m odel (RBIM) fCho and Fishei^, 
Il997a. : Merz and Chalkerll2002 ) , which is obtained if all 
the additional vortices are placed on the same sublattice. 

The phase d iagram of the CF model was established 
in IChalker et al . (.2002), and it was found that all the 
three expected phases are indeed present, Fig. [35l The 
DOS in these phases and at transitions between them was 
studied in Mildenbcrgcr et al. (2007a), the results will be 
presented in Sec. | VI.E.3I IVI.E.4I It was also checked in 
iMildenberger et aU (i2007a) that the CF model is indeed 
generic: the same behavior is obtained for a model with 
a Gaussian distribution 'P{0). 

We briefly discuss now two other disorder models, with 
properties qualitatively different from the CF model: 

(i) A fermionic version of the zt J RBIM is described 
by a disordered network model with I Chalker et aLl . [2002t 



iMerz and Chalked . [200I iRead and Ludwid . l200lh 

V{0) = (1 - p)S{9 - 60) + pS{0 + 60). (6.48) 



This implies that all pairs of vortices are inserted 
in th e same sublattice. It has been shown analyti- 
callv (iRead arid Ludwi3 . l200lh and verified numerically 
(IChalker et a^.l ■ l2002^ that the metallic phase is absent in 
the RBIM. Two phases with localized states (separated 
by the TQHE transition) correspond to the paramag- 
netic and ferromagnetic phases in the Ising spin language. 
The self-dual state of the disorder-free network {p = 0, 
sin^ 9 = 1/2) maps onto the critical point of the clean 
Ising model. 

(ii) The 0(1) model is obtained if one includes in the 
regular network factors (-1) for propagation along some 
links, randomly selected with a concentration p. The cru- 
cial feature of such a defect (that distinguishes it from 
the randomness in the nodal parameter such as the vor- 
tex pairs in the CF model and the RBIM) is that it 
cannot be "switched off" by any continuous transforma- 
tion. Therefore, such a defect has the topologica l char- 
acter of a vortex. It was found (Chal ker etdl . 120021: 
iRead and Ludwi^l200lh that such topological defects de- 
stroy completely the localization, so that the 0{1) model 
is always in the metallic phase. From the cr-model per- 
spective, it was shown that the effect of vortices is in 
suppression of t he second (disconnect ed) component of 
the target space ( Bocauet et a/.l . [20001 ). 



3. Thermal metal 

The metallic phase can be treated analytically by using 
the (T-model approach. The corresponding RG analysis 
yields 



dt/dlnL 



-t' 



(6.49) 



where t is the running coupling constant inversely pro- 
portional to the dimensionless conductivity, t = I/ttj?. 
The infrared behavior of the system is governed by the 
perfect-metal fixed point, t —> 0. Specifically, the con- 
ductance increases logarithmically with the system size, 
g{L) — go + (I/tt) lnL/£o {£0 is the mean free path), jus- 
tifying the perturbative RG. The RG equation for the 
second coupling constant e, whose bare value is given by 
the energy E, reads 



de/d\iiL = {2 + t)e 



(6.50) 



leadi r ig to a logarithmic increas e of DOS ( Bocauet et all . 
I2OOOI : ISenthil and Fisherl . |2000D . 



In 



D 



Att^D \E\£l 



(6.51) 



where D is the diffusion constant (remaining unrenor- 
malized to this order), go = 2npQD. 
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FIG. 36 Low-energy DOS in the metallic phase. Parameters 
(upper curves): p—0.5, Q=7r/4, system sizes L=128 (squares) 
and L—256 (full circles). The straight dashed line represents 
the logarithmic asymptotics. For lowest energies the RMT 
oscillations are clearly visible; they can be collapsed on a sin- 
gle curve as shown in Fig. 1371 For comparison, the results for 
p — 0.15 and p — 0.1 are also shown; the latter point is close 
to the expected boundary of the metallic phase, see Fig. 1351 
It is seen that when the system approaches the phase bound- 
ary, the logarithmic increase of the DOS disappear s and t he 
RMT oscillations get damped. (|Mildenberger ef aZ.I. l2007al '). 



Numerical results for the DOS in the metallic phase 
are shown in Fig. 1361 The data exhibit a logarithmic in- 
crease of the DOS over almost three decades in E for the 
larger system size, L=256. It is worth stressing that the 
increase continues to be of logarithmic form even though 
the renormalized DOS at small energies becomes much 
larger than its bare (large-£^) value po— O-l- This is a 
signature of the fact that the RG flow is towards weak 
coupling, so that the one-loop result (|6.5ip is valid down 
to arbitrarily low energies in the thermodynamic limit. 

At the smallest energies, pronounced oscillations in the 
DOS are observed. These are RMT oscillations due to 
finite system size and serve as another indication of the 
fact that we are dealing with a metallic phase. The RMT 
origin of these oscillations is demonstrated in Fig. [571 
where these parts of DOS curves are replotted, with the 
energy rescaled to the mean level spacing 5l at lowest en- 
ergy for the corresponding system size. The data collapse 
on a single curve, which shows that the renormalized level 
spacing 



1 



1 



(6.52) 



with E^h being the Thouless energy, is indeed the only 
relevant energy scale in the regime E < E^h where the 
RMT is applicable. As further seen in Fig. [371 the ob- 
tained curve agrees with the RMT prediction, 



PiE) = 



1 



L^Sr 



sm{27rE/SL) 
2ttE/Sl 



(6.53) 




0.75 - 



FIG. 37 Renormalized DOS at maximal disorder p=0.5 and 
on the symmetry line sin^a=l/2 for different system sizes 
vs. the energy measured in units of the level spacing 5l- 
The RMT result, Eq. ([6]53}, is plotted as a dashed line 
for comparison. Inset: logarithmic dependence of 1/L^5l 
on the system size L, con sistent with the data of Fig. 1361 
([Mildenberger et allWOTei ) . 



up to _E/(5l^ 1.5-2; for larger energies the oscillations are 
strongly suppressed. This is fully consistent with the ex- 
ponential vanishing of the RM T oscillations b eyond the 
Thouless energy, see the review (lMirlinl . l2000bl) . With in- 
creasing system size, the ratio Etii/Sl increases (though 
only logarithmically), so that the RMT range includes 
progressively more oscillation periods. This tendency is 
clearly seen in Fig. [571 



4. Localized phases and TQH transition 

An analytical approach to the problem alternative to 
the cr-model is based on the model of Dirac fermions with 
random mass, in the spirit of th e analy sis of the Ising 



ranaom mass, m tne spirit oi tn e analy sis or t 
model in iDotsenko and Dotsenkol ( 1983( ). Being pertur- 
bative in the disorder strength, this approach is appro- 
priate for the description of the localized phases and the 
transition between them. The disorder-free system has a 
transition, driven by tuning a uniform mass through zero, 
which in the CF model lies at p=0, sin^ 0=1/2 and cor- 
responds to the clean Ising transition. In the vicinity of 
the clean fixed point representing this transition the dis- 
order strength gm is marginally irrelevant. This implies 
for th e critical DOS a loga r ithmic correction te r m of t he 
form (|Bocauet et all . l2000l : [Mildenberger et all [2007al) 



p{E) 



27r 



2.9A/ , 1 

In ; 

TT E 



(6.54) 



Slightly away from the critical value 9 = ir/i, the system 
is in a localized phase with a large localization length 
^ cx \d — 7t/A\~'^. As the RG flow is towards the clean 



48 



Ising fixed point, the corresponding index should be the 
same as in the Ising model, = 1. The behavior of 
the DOS in the localized ph ases can be also understoo d 
using the Dirac fermion RG ( Mildenberger et aZ.l . l2007ah . 
Specifically, for energies that are not too small, behavior 
will be the same as at criticality, Eq. ()6.54p . However, 
for smallest energies, it is the localization length ^ (rather 
than E) that will terminate the RG process. In this sense, 
the role of ^ is fully analogous to that of finite system size 
L at criticality. This implies that p{E) saturates at the 
value 



P(i^)-|(l + 2^1nf 



1/2 



E<E^, (6.55) 
The energy E^ at which the saturation takes place is 



E, 



£ \ -1/2 

l + 2^1n-^^ 



(6.56) 



The Anderson transition from the metallic to one of the 
localized phases should therefore be signalled by a tran- 
sition from logarithmically diverging to finite p{E — > 
0). This has been verified by numerical simulations in 
iMildenberger et all (|2007al ). Fig.[38l 

A brief comment on the regions of localized phases 
where the interplaquette coupling is very weak (sin^ 9 
close to zero or to unity), i s in o rder here. As shown 
recently ( Mildenberger et a/.l . l2006l ). in this situation the 
DOS of the RBIM acquires a non-universal power-law 
singular ity, with z > 1 associated with Griffiths 

strings (jMotrunich et all [2OO1I . ,2002). The same mech- 
anism of the formation of divergent DOS in these parts 
of the localized phases is expected to be operative in the 
CF model as well. 

We turn now to numerical results I Mildenberger et all . 
l2007al ) on the DOS at the line sin^ 6* = 1/2, where the 
TQH phase boundary is located. They are shown in 
Fig. [311 While at sufficiently large p (p > 0.1) the 
DOS shows a logarithmic increase characteristic for the 
metallic phase, for lower p the DOS behavior agrees with 
Eq. (|6.54p . as expected for the TQH transition. This is 
demonstrated in the inset of Fig. [311 where p{E)/\E\ as 
a function of log \E\ is plotted for p = 0.05. 

While at moderately low E the DOS at the TQH tran- 
sition line is in good agreement with the Dirac-fermion 
RG, Eq. (16.541). the results for the lowest energies ob- 
tained in IMildenberger et all { 2007al ) constitute a puzzle. 
Specifically, it was found that the DOS saturates at a very 
low energy scale and even shows a weak upturn. The rea- 
son for this behavior is not understoo d at present; several 
possibl e scenarios were proposed in [Mildenberger et all 
( 2007a[ ): (i) the position pT of the tricritical point T is 
in fact not pr — 0.1 as in Fig. [35lbut rather considerably 
smaller, pT < 0.05; (ii) in addition to the tricritical point 
Pt there is a second, repulsive fixed point on the TQH 
transition line svc? 9 = 1/2 , at some pn < pt- This 
point would then act as a "flow splitter" which is similar 
to the role of the Nishimori point in the RBIM; (iii) the 
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FIG. 38 DOS near £ = for disorder values p = 0.13, 0.08, 
0.04, and 0.02 and for two system sizes L = 64, 128 at fixed 
interplaquette coupling sm^{9) = 0.579. The DOS diverges 
logarithmically as i5 ^ in the metallic phase {p = 0.13) 
and remains finite in the localized phase (other values of p). 
The results for the lowest impurity concentration, p = 0.02, 
show an oscillatory feature induced by the band structure of 
the clean system, as well strong scatter in the data at the 
lowest energies, which is due to insufficient ensemble averag- 
ing. Upper inset: Low-energy peak at p = 0.13; its amplitude 
increases with L, in agreement with Sec. IIII.AI Lower inset: 
Low-energy DOS at p = 0.08. No peak at i? ^ is detected; 
p{E 0) is a constant independent on L, indicating that 
the system is in the insulating phase. Statistical noise in the 
lower inset is more pronou nced than in the upper one due to 
the smallness of the DOS. ([Mildenberger et ad[2007al ). 



RG treatment of the theory of Dirac fermions with Gaus- 
sian random mass is in fact insufficient, and some effects 
- possibly of non-perturbative origin - eventually drive 
the system away from the clean Ising fixed point. The 
clarification of this important issue remains a subject for 
future research. 



F. Chiral classes (AMI, Cll, BDI) 



We have seen in Sec. IV. El that quasi-lD systems of 
chiral symmetry show a criticality accompanied by a very 
slow (L~^/^) decay of the average conductance and by 
a Dyson- type singularity (l/|£;in^ £;|) in the DOS. As 
we discuss in the present subsection, a similar type of 
criticality takes place also in two dimensions. 



1. Gade-Wegner cr-model 



In their 



[Gade and Wegnerl 
tems of the chiral 



pio neering 
19911) 



works, [Gadd ([1993[) : 
derived cr-models for sys- 
classes and performed their RG 
analysis at and near two dimensions. They used the 
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FIG. 39 DOS at low energy on the self-dual line sin^(Q) = 0.5 
for disorder concentrations p — 0.2 (o, o), 0.1 (□, x), and 0.05 
(A, where in each case the first symbol is for L — 128 
and the second is for L — 256. Inset: p{E)/\E\ at p=0.05 
on a log-linear scale. The logarithmic correction is clearly 
observ ed, in agreement with Eq. (|6.54p . (jMildenberger et all 
l2007al '). 



fermionic replicas, so that the models are defined on the 
spaces U(7V), U(2iV)/Sp(2iV), and U(A^)/0(iV) for the 
chiral unitary, orthogonal and symplectic classes, respec- 
tively. As usual, the supersymmetric generalization (see 
Table U for the symmetry classes of the corresponding 
models) is equivalent to the replica version on the level 
of the perturbation theory. For definiteness, we consider 
the chiral unitary class (AIII), the results for the other 
two are very similar. The action of the cr-model has the 
form 



S[Q] = / d-r 



^ TrVQ-^VQ - :7^[Tr(Q-iVQ)]2 



+i 



16nt 6Attc 
;^eTr(Q + 0-i)}, 



(6.57) 



where Q G U(iV), t and c are two coupling constants, 
the first of which is related to the conductivity in the 
usual way, t^^ = 2Traxx, e is the running coupling 
whose bare value is the energy E (which breaks the 
chiral symmetry), and po the bare density of states. A 
special feature of the tr-model for chiral classes is the 
existence of the second term (known as Gade term), 
which governs fluctuations of det Q. This cr-model 
was later obtained and analyzed in a number of works 



was later optamea a. na anaiy ze a m a number oi worJts 
(lAltla.nd and S imonj. 19991: Fabrizio and Gastellanil . 



I2OOOI : lFukuil ."fr999: Gu ruswamv et aLl . I2OOOI ) 



The one-loop RG equations in 2 + e dimensions read 

-dt/\iiL = et (exact!); (6.58) 
-dc/d\nL = €0+20^ ■ (6.59) 
dlne/dlnL = (2 + e) +tV2c . (6.60) 

It is of central importance that in 2D (e = 0) the 
/?-function (|6.58p is identically zero ( Gade and Wegnerl . 



Il99ll ). i.e. the coupling constant t (and thus, the con- 
ductivity) is not renormalized. If one considers a system 
at non-zero (but small) energy, the chiral symmetry is 
broken down to that of the unitary Wigner-Dyson class 
(A), and the localization should happen at a sufflciently 
large scale. To find the scaling of the localization length 
with energy, one performs the RG transformation until 
the running energy e{L) ceases to be small (i.e. reaches 
a characteristic scale A of the problem; typically, the 
bandwidth). Beyond this scale, the RG takes the form 
characteristic for class A, driving the system towards the 
localized regime. The localization length £,{E) is thus 
given by the crossover length Lc{E), up to an energy- 
independent factor ~ exp(— l/4i^). Further, the DOS is 
given by 

p{E) ^ l/ELliE). (6.61) 
Integration of Eqs. (HSH), ([OP)) (with e = 0) yields 

\nL,{E) ^ + 4i2| ln(£;/A)| - b) , (6.62) 

where B ^ 2 + t^/co and co is the bare value of the cou- 
pling c. For asymptotically low energies, |ln(_E/A)| ^ 
B^/4t'^, this reduces to 



^{E)ocexp[t-^\\n{E/A)\^/^^ , (6.63) 
p{E) cx -B-i exp (-2t-^\\n{E/A)\^/^) . (6.64) 



The last formula is the 2D counterpart of the Dyson sin- 
gularity. On the other hand, at intermediate energies, 
|ln(i?/A)| ^ B^/4t^, the behavior is of the power-law 
type with a non- universal exponent. 



^E) oc ; p{E) oc E-'+^Z'' . 



(6.65) 



It is worth emphasizing that the bare coupling cq is 
of the order of unity even for a system with large con- 
ductance {t ^ 1). This means that B ~ 2 and the 
asymptotic behavior ()6.63|) . (I6.64p establishes at expo- 
nentially low energies only, ln(A/£') ^ 1/t'^. For this 
reason it is very difflcult to reach the true asymptotics in 
numerical simulations, whereras the intermediate power 
law regime (16.651) can be st udied very well, see e.g. 



iMarkos and Schweitzel ([2003). 

Like in the ID geometry, one can drive the system 
away from criticality by introducing a staggering in hop- 
ping strength ( "dimerization" ) in a lattice model. The 
DOS in the localized phase was predicted to show a 
power-law behavior, p{E) ~ governed by the 

Griffiths mechan ism, with a non-univers al dynamical ex- 
pone nt z > (Motrunich et_ Ml, I2002D . In the work 
( Bocquet and Ch alkcr. 2003) network models for all the 
chiral classes were constructed. A numerical study of 
the class- AIII network model confirmed the existence of 
the critical (Gade-Wegner) and the localized (Griffiths- 
type) phases. It was found that the DOS exhibits a non- 
universal power-law behavior in the localized 
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phase and that z asymptotically tends to infinity with 
decreasing energy in the critical phase, in agreement with 
analytical predictions. 



2. Dirac fermions approach. Strong-coupling effects 

The Gade-Wegner prediction of a critical state with 
a diverging DOS of the form (|6.64p can also be re- 
produced b};_sta£yiig_Jom of disordered Dirac 
fermions ( Guruswamv et a/] . l2000t) . Specifically, a model 
of fermions o n a bipartite lattice w ith tt flux and random 
real hopping f|Hatsugai tt aLl . [l997l ). which belongs to the 
class BDI, is described by a two-flavor model of Dirac 
fermions subject to a random vector potential (coupling 
qa) and a chiral random mass (gm)- This is a partic- 
ular type of two-flavor disordered Dirac fermion models 
considered in Sec. IVI.G.41 in the notations used there 
gA = 7_l/2, Sm = /3?;/2 . For this model, the ex act /3- 
functions can be found ( Guruswamv et oil . l2000f ) to all 
orders in Qm'- 



dgA ^ 

dlni ~ (1 + 25,„)2 



dXnL 







(6.66) 



Further, the dynamical scaling function was found to be 



dine 
ci In L 



1 + 



(1 + 25^)2 



2.g„ + 0{gl 



(6.67) 



According to Eqs. (j6.66p . gm is not renormalized, while 
gA grows logarithmically with L. Assuming for simplicity 
weak disorder, gm 1, we have 



(6.68) 



For sufficiently large L the InL term dominates. Substi- 
tuting it in Eq. (|6.67p . one gets 



Ine(L) ~ 2gl^\n^L 



(6.69) 



This reproduces Eqs. (|6.63p . (|6.64p for the localization 

length and the DO S with t qm^/^ - 

As was shown in lMotrunich et cd. 1^2002'), this result is, 
however, not fully correct, for the following reason. It is 
known that the multifractal spectrum of Dirac fermions 
in a random vector potential undergoes a transition from 
a weak-disorder to strong-disorder phase at (/^i = 1 
("freezing transition"), see Se c . IVI.G.31 As shown i n 
iHorovitz and LeDoussal (|2002h : [Motrunich et aTl (|2002h . 
this transition is accompanied by a change of the behav- 
ior of the dynamical exponent 




5A < 1 ; 
5A > 1 ■ 



(6.70) 



In the presence of the second coupling {gm), 9a{L) will 
flow according to ()6.68p . and z will develop following 
Eq. (p7D)) . While the first line of (p7D)) agrees with 



Eq. (|6.67p . it is only valid at short distances, gA{L) < 1 
(which corresponds to the ballistic regime). At suffi- 
ciently long (diffusive) scale, where the cr-model is ap- 
plicable, gA will become larger than unity. Using the 
second fine of Eq. (|6.70p instead of Eq. (|6.67p yields 

p(i?)~i?-iexp|-i(3ff^i|ln(i?/A)|)2/3| . (6.71) 



Therefore, the exponent 1/2 in (|6.63p and in the second 
(subleading) factor in (|6.64p is replaced by 2/3. 

This result was also rederived in the framework of the 
RG approach. The key point is that, when the running 
coupling gA ceases to be weak, it is not sufficient any- 
more to characterize the disorder by the lowest-order cu- 
mulants. Instead, one should take into account an infi- 
nite number of couplin gs. Such a functional RG method 
was developed in .Garpentier and Le Doussal (1999) in the 
context of a related random XY model. For the problem 
of Dirac fermions with random ma ss and vect or poten- 
tial this program was carried out in lMudrv et a l. (20031): 
the re sults confirm the conclusion of iMotrunich et al\ 
( 20021). Eg. (16.641). Analogou s results were also obtained 
in lYamada and Fukuil (|2004f) . 



G. Disordered Dirac Hamiltonians 

Localization and criticality in models of 2D Dirac 
fermions subjected to various types of disorder have 
been studied in a large number of papers and in a va- 
riety of contexts, including the random bond Ising model 
(Dotsenko and Dot senkol . ll983( ). the quantum Hall effect 
(jLudwig et al, 199 ^1 dirty super c onduc t ors with uncon- 
ventional pairing TAltland et all l2002t iBocauet et all . 
I2OOOI: iNersesvan all Il995t). and some lattice mode s 
with chiral symmetry ( Guruswamv et a/1 l2000f ). Re- 
cently, this class of problems has attracted a great 
deal of attention in connect ion with its application to 
graphen e, see, in particular, Aleiner and E fctov (2006}; 
LAltiandi (!200 7l): iKhvesh chcnko (2006[): iMcCann et M 



f|2006l) : lOstrovskv et al\ (,2006„ ^2007anbl ). 

In the presence of different types of randomness, Dirac 
Hamiltonians realize all ten symmetry cl a sses of disor- 
dered systems; see ( Bernard and LeClaiil . l2002al) for a 
detailed symmetry classification. Furthermore, in many 
cases the Dirac character of fermions induces non-trivial 
topological properties (0-term or WZ term) of the cor- 
responding field theory (ct- model). In Sec. IVI.G.ll we 
review the classification of disorder in a two-fiavor model 
of Dirac fermions describing the low-energy physics of 
graphene, the RG treatment of the problem, and types 
of criticality. The emergent critical theories will be dis- 
cussed in Sec. IVLG.2ffVLG:4l In Sec. IVI.G.5I we briefly 
discuss a four-node Dirac fermion model appropriate for 
a description of dirty d-wave superconductors. 
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1. Disordered two-node Dirac Hamiltonians: Symmetries of 
disorder, renormalization group, and types of criticality. 

The presentat ion below largely follows the papers 
(|Ostrovskv et all [2006, 2007b). We concentrate on 
a two-flavor model, which is in particular relevant to 
the description of electronic properties of graphene. 
Graphene is a semimetal; its valence and conduction 
bands touch each other in two conical points K and K' 
of the Brillouin zone. In the vicinity of these points the 
electrons behave as massless relativistic (Dirac-like) par- 
ticles. Therefore, the effective tight-binding low-energy 
Hamiltonian of clean graphene is a 4 x 4 matrix operating 
in the AB space of the two sublattices and in the K-K' 
space of the valleys. We introduce the four-component 
wave function 

* = {(t)AK,(t)BK, (jiBK'AAK'}'^ , (6.72) 

In this representation the Hamiltonian has the form 



H = woTacrk. 



(6.73) 



Here T3 is the third Pauli matrix in the K-K' space, or — 
{ci, (J2} the two-dimensional vector of Pauli matrices in 
the AB space, and vq the velocity (uq — 10^ cm/s in 
graphene). It is worth emphasizing that the Dirac form of 
the Hamiltonian (|6.73p does not rely on the tight-binding 
approximation but is protected by the symmetry of the 
honeycomb lattice which has two atoms in a unit cell. 

Let us analyze the symmetries of the clean Hamiltonian 
(I6.73P in the AB and KK' spaces. First, there exists an 
SU(2) symmetry grou p in the s pace of the valleys , with 
the generators (jMcCann et al lllooi) 



^X—CTZTI, = 0-3X2, = CToTa 



(6.74) 



all of which commute with the Hamiltonian. Second, 
there are two more relevant symmetries of the clean 
Hamiltonian, namely, time inversion operation (Tq) and 
chiral symmetry (Co). Combining Tq, Cq, and isospin ro- 
tations Ao,a;,j,,2, one can construct twelve symmetry op- 
erations, out of which four (denoted as T^) are of time- 
reversal type, four (C^) of chiral type, and four {CT^) of 
Bogoliubov-de Gennes type: 



To : Ai-^aiTiA ctiTi, Cq : A 1-^ -a^ToAa^To, 

CTq : A ^ -cF2TiA^ a2Ti, 
T^: A^ a2ToA^a2To, C^: A^ -cron^o'oTi, 

CT^ : A ^ -aiToA^aiTo, 
Ty : Ai-^ a2T3A^a2T3, Cy : At-^ -aoT2AaoT2, 

-asT^AasTs, 



: A 



CTy : A ^ -aiT^A^aiTs, 
aiT2A^ (J1T2, Cz'- Ai- 
CTz : A ^ -a2T2A^ a2T2- 



It is worth recalling that the C and CT symmetries apply 
to the Dirac point {E = 0), i.e. to undoped graphene, 
and get broken by a non-zero energy E. We will assume 



the average isotropy of the disordered graphene, which 
implies that A^, and Ay symmetries of the Hamiltonian 
are present or absent simultaneously. They are thus com- 
bined into a single notation A^; the same applies to Tl 
and C^_. In Table |V] all possible matrix structures of 
disorder along with their symmetries are listed. 

To derive the field theory of the problem, one intro- 
duces a superfield '(/'(r), see Sec. lII.B.TK or. alternatively, 
its replica counterpart). After the disorder averaging, 
the action containing all possible disorder structures from 
Table El reads [ip = i/'U) 



Pa - 

+ az{ipcr3Totp)'^ + Y [(V'CTanV')^ + {tpcr3T2tp)'^] 

+ Y [(V'O-inV')^ + (?/'0-lT2'(/')^ + {'llja2Ti1py + {■lpa2T2lpy 
+ Y [(V'O-on'0)^ + ('0croT2'(/')^] + loi^CroTjIp)'^ 



(6.75) 



A complete set of one-loop perturbative RG equa- 
tions for nine disorder amplitudes has the form 
(|Ostrovskv et all l2006t) 



d\nL 



d\nL 

duz 
d\nL 



dpQ 
d\nL 



dP± 
d\nL 

dPz 
d\nL 



d'jo 
dlnL 



dj± 
d\nL 

d-fz 
d\nL 



= 2ao(Q!o + /3o + 70 + a_L + /3i + 7± 

+ az + + jz)'2a±az + (}±f3z + "^li^lz] 

= 2(2aoaz + I3qI3z + 2joJz); 

= -2az{ao + /5o + 70 - a± - - 7-L 

+ az + l3z+ Iz) + 2aQa± -f Po[3± -f 2707^; 

= 2[/3o(ao - 70 + a_L + - 7^) 

+ a^Pz + azl3±_ -f /3i7o]; 

-iiaoPz + azl3o -\- /3o7o + /3_l7-L + Pzlz)\ 

^ 2[~(3z{ao - 70 - a_L + ttz - Iz) 

+ ao/3_L + a_L/3o + /3_l72]; 

= 27o(ao - /3o + 70 + a± - /3i + 7± 

+ otz- Pz+ Iz) + ^o-i^lz + 20^7^ -(- /?o/3±; 

= 4ao7z + 4a^7o + /3o + /?! + /^f; 

= -27^(qo - a± + - /3o + /3± - /3z + 70 



- 7-L+7^) + 2ao7± + 2ai7o + /?J-/3;.- (6.76) 
The RG equation for the running energy e (whose bare 



52 



structure 


coupling 


AE 


GLL 


Ax 


A. 


To 






Co Cx 


Cz CTo CT± 


CTz 




QO 


7o/27ri;2 




+ 


+ 


+ 


+ 


+ 








0"{1,2}''"{1,2} 


/3x 


2/3_L/7rt;2 




- 


- 


+ 


- 


- 


+ 


- + - 


- 


O"l,2T0 


7x 


'y±/Tvv'^ 




- 


+ 


+ 


- 


+ 


+ 


+ + - 


+ 


fori, 2 


P. 




2gm 


- 


- 


+ 


- 


- 


- 


+ - - 


+ 


0"3T3 


Iz 


7z/27ru^ 




- 


+ 


+ 


- 


+ 


+ 


- - + 


- 


0"3Tl,2 


l3o 




23m 










+ 




+ + - 




o'oTa 


70 








+ 
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TABLE V Disorder symmetries in graphene. The first five rows of tlie table represent disorders preserving tfie time reversal 
symmetry To; the last four — violating Tq. First column: stru cture of disorder in the sublattice (cr^) and valley (ti,) spaces. 
Second column: disorder strength according to the notations of (lOstroyskv et a;.ll2006l. i2007b.l used i n this review. Third (AE) 
and fourth (GLL) columns: disorders considered in l|Aleiner and Efetovl . l2006l ) and (fCuruswamv et all\2QQ(i ) a nd notations used 
there. The remaining columns indicate which symmetries of the clean Hamiltonian are preserved by disorder. l|Ostrovskv et alX 
120061 '). 



value is E) reads 
dine ^ 



dlnL 



(6.77) 



If all types of disorder are present (i.e. no symmetries 
is preserved), the RG flow is towards the conventional 
localization fixed point (unitary Wigner-Dyson class A). 
If the only preserved symmetry is the time reversal (Tq), 
again the conventional locali zation (orthogonal Wigner - 
Dyson class AI) takes place ( Aleiner and Efetovl l2006l ). 
A non-trivial situation occurs if either (i) one of chiral 
symmetries is preserved or (ii) the valleys remain decou- 
pled. In the Table IVII we list those situations when the 
symmetry prevents the localization and leads to critical- 
ity and non-zero conductivity at i? = (in the case of 
decoupled nodes - also at nonzero E). Models with de- 
coupled nodes will be analyzed in Sec. IVI.G.2i and mod- 
els with a chiral symmetry in Sec. IVI.G.3I (Co-chirality) 
and Sec. IVI.G.4I (C^-chirality). 



2. Decoupled nodes: Disordered single-flavor Dirac fermions 
and quantum-Hall-type criticality 

If the disorder is of long-range character, the valley 
mixing is absent due to the lack of scattering with large 
momentum transfer. The couplings that do not mix the 
valleys are ao, q;_l, az, 7o, 7±, and 7^. For each of the 
nodes, the system can then be described in terms of a 
single-flavor Dirac Hamiltonian, 

H = vo[(T'k + aM^)\- (6-78) 

Here disorder includes random scalar (Vo) and vector 
(^1.2) potentials and random mass (V3). The corre- 
sponding couplings are 5y = ao -|- 70, 2^^ = ax + 7±, 



and Qm = ctz + Iz ■ The clean single- valley Hamilto- 
nian (|6.78p obeys the effective time-reversal invariance 
H = (T2-ff"^cr2- This symmetry (Tx) is not the physical 
time-reversal symmetry (Tq): the latter interchanges the 
nodes and is of no significance in the absence of inter- 
node scattering. The RG equations have the form 

dgA/dlnL = 2gvgm ; (6.79) 

dgv/dlnL = 2{2gA + gv){gv + g,n) ] (6.80) 

dgm/dliiL = 2{2gA- gm)igv + gm) ; (6.81) 

d\ne/d\nL = 1 + 2gA + (jm + gv ■ (6.82) 

Remarkably, single-flavor Dirac fermions are never in 
the conventional localized phase! More specifically, de- 
pending on which of the disorders are present, four dif- 
ferent types of criticality take place: 

(i) The only disorder is the random vector potential 
(g^). This is a special case of the symmetry class AIII. 
This problem is exactly solvable and has been studied in 
a great detail; we will discuss it in Sec. IVI.G.3I 

(ii) Only random mass {cjm) is present. The system 
belongs then to class D. The random-mass disorder is 
marginally irrelevant, and the system flows under RG 
towards the clean fixed point. This problem has been 
analyzed in Sec. IVI.E.3I 

(iii) The only disorder is random scalar potential (gv)- 
The system is then in the Wigner-D yson symplectic 
(All) s ymmetry class. As was found in lOstrovskv et al\ 
( 2007a( ). the corresponding cr-model contains a topologi- 
cal term with 9 = ir, which leads to delocalization. This 
model has been discussed in Sec. I VLB. 51 



" (i) The tilde 

in 5m serves to distinguish it from the "chiral ran- 
dom mass" coupling gm, see Table FVl and Sec. [VLF.21 IVI . G.4I 
(ii) our couplings are related to those of lLudwig et al\ ||1994^ via 
gA = Aa/Ztt, gv = Av/27r, g,„ = Am/2n. 
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Symmetries 


Class 


Criticality 


Conductivity 


Cz, To 


BDI 


Cade 






AIII 


Cade 






CII 


Cade 


PS 4e^/-Kh 


Co, To 


CI 


WZW 




Co 


AIII 


WZW 




Az, Co 


2 X AIII 


WZW 


4e^/Trh 


Az, Tx 


2xAII 


e = TT 


4crJ* or oo (?) 


A., CTx 


2xD 


e = TT 




Az 


2xA 


e = TT 





TABLE VI Possible types of disorder in graphene leading 
to criticality." First column: preserved symmetries. Second 
column: symmetry class. Third column: type of criticality. 
The first three row correspond to Cz chiral symmetry lead- 
ing to Gade-Wegner-type criticality, Sec. IVI.G.4I The next 
three rows contain models with Co chiral symmetry (ran- 
dom gauge fields), inducing a WZ term in the a-model ac- 
tion, Sec. IVI.G.3I The last four rows correspond to the case 
of decoupled valleys (long-range disorder), see Sec. IVI.G.2I 
from top to bottom: random vector potential, scalar poten- 
tial, mass, and any of their combinations. For these models 
(except for the random vector potential, which is Co chiral at 
the same time) the a- model acquires a topolog ical term with 
6 ^Tv. Adapted from (|Ostrovskv et gZ.I . l2007bf i. 

"A further possible mechanism of delocahzation is generation of 
the metallic phase due to broken spin- (or isospin-) rotation in- 
variance, Sec. I VI. A. II If the only preserved symmetry is Tz , the 
system is in class All, while if only CTo or CTz invariance is pre- 
served the system belongs to class D. (No topological term arises in 
these situations.) In both these classes, the system flows towards 
a perfect-metal fixed point if the bare conductivity is above the 
localization threshold. 



(iv) At least two types of randomness are present. This 
is the same as to say that all of them are present, as 
the third one will be generated by RG, see Eqs. (|6.79p - 
(|6.8ip . The model belongs to the Wigner-Dyson un itary 
class A, and it was argued in lLudwig et all i 1994J) that 



it flows into the IQH transition fixed point. This is con- 
firmed by the deriva t ion of the correspondin g g-model 
(jAltland all . l2002t lOstrovskv et all . l2007allbl ). which 
contains a topological term with 9 = tt, i.e. is noth- 
ing but the Pruisken u-model at criticality. A particular 
consequence of this is that the conductivity of graphene 
with this type of disorder is equal to the value a^ of the 
longitudinal conductivity axx at the critical point of the 
IQH transition multiplied by four (because of spin and 
valleys). Note that the conclusion of IQH criticality for- 
mally holds for arbitrary energy e, although in reality it 
only works near half-filling; for other s the critical point 
would only be approached for unrealistic system sizes and 
temperatures. 

If a uniform transverse magnetic field is applied, the 
topological angle 9 becomes energy-dependent. However, 
at the Dirac point {E = 0), where axy = 0, its value 
remains unchanged, 9 = tt. This implies the emergence 



of the half-integer quantum Hall effect, with a plateau 
transition point at E = 0. 



3. Preserved Co chirality: Random abelian and non-abelian 
vector potentials 

a. Generalities. Let us consider a type of disorder which 
preserves the Co-chirality, H = —a^Hu^] according to 
the Table |Vl the corresponding couplings are a±_ , (3± , 
and 7_L. The one-loop RG equations then read 



da±/dlogL — 0, 
9/3i/91ogL = 4/3i7i, 
97i/91ogL = /3i . 



(6.83) 
(6.84) 
(6.85) 



The specifics of the Co-chirality is that it corresponds 
to the more common meaning of the term "chirality" , 
which refers to a distinction between "left" and "right" 
particles. The model takes then the form of the euclidean 
version of 1 -I- 1 theory of massless Dirac fermions coupled 
to a gauge field. Indeed, according to (|6.73p the matrices 
aiT3 and CT2T3 play the role of Dirac 7-matrices, and the 
matrix CT3 determining the Co-chirality is nothing but 
75. Depending on further symmetries, three different Co- 
chiral models arise: 

(i) The only coupling present is a± , which corresponds 
to the random abelian vector potential. In this case the 
nodes are decoupled, and the Hamiltonian decomposes in 
two copies of a model of the class AIII. This model has 
already been mentioned in Sec lVI.GT2] 

(ii) If the time-reversal symmetry Tq is preserved, only 
the couplings P± and j± are allowed, and the problem 
is in the symmetry class CI. The model describes then 
fermions coupled to a SU(2) non-abelian gauge field, and 
is a particular case of analogous SU(N) models. 

(iii) All three couplings are present. This describes 
Dirac fermions coupled to both abelian U(l) and non- 
abelian SU(2) gauge fields. This model is in the AIII 
symmetry class. 

Remarkably, all these critical Co-chiral models are ex- 
actly solvable: one can calculate exactly the spectra of 
multifractal exponents, the critical index of the DOS, and 
the value of conductivity. They have been studied in a 
large number of works, and we summarize the key ideas 
and results. 



b. Abelian vector potential. The model of 2D Dirac 
fermions moving in a random vector potential, 



(6.86) 



is exactly solvable at zero energy and ch aracterized by a 
line of fixed points ( Ludwig et all . ITOQ^ ). Decomposing 
the vector potential in the longitudinal (pure gauge) and 
transverse parts, 



A. 



(6.87) 
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and assuming that the total magnetic flux is zero^'^, 
one can expHc i tly w r ite the zero-energy wave functions 
(|Castillo et all Il997l : iLudwig et a"H ■ ll994^ : 



*+(r) 



V'±(r) = S±'/'e-''^("")e±^("") , (6.89) 

where B± — f d^re^^'^*^''^ is the normahzation factor. We 
consider the first of the functions (|6.88p for definiteness 
and analyze the statistical properties of |'(/'+(r)p. They 
are governed by fluctuations of (/'(r). Assuming a white- 
noise distribution of the random vector potential with 
{Af^{r))A„{r')) = 27rg^(5p^(5(r— r'), one gets the following 
statistics of 4>{r), 



V[(l)] cx exp 



-1 



(6.90) 



Equation (|6.90p implies that (j>{r) is a free massless 
bosonic field characterized by the correlation function 



(0(r)0(r')) = .gAln(L/|r-r'|) 



(6.91) 



c. Multifractality. To find the multifractal spectrum, one 
considers the moments (|i/)+(r)|)^'^, see Sec. III.CTl For 
not too strong disorder {qa < 1), it turns out to be suf- 
ficient to average separately the exponential e^"^"^^""^ and 
each of the q nor malization factors [ c Pre^^'^^'^K The 
resulting spectrum ( Ludwig et a/.l . [l9M ) 



= 2iq ^ - gAq) 



(6.92) 



has an exactly parabolic form. The corresponding f{a) 
spectrum is given by Eq. (|2.37p with d — 2, ^ = 2gA- 



fia) = 2 - (a - 2 - 2gA)V8g, 



(6.93) 



In view of the exact parabolicity, the spectrum necessar- 
ily contains a termination point at a = 0, see Sec. III.C77l 
The corresponding value of q is q^ = (1 + gA)/'^9A- 
Thus, Eq. (|6.92p is only valid for q < qc, for larger q 
the exponent Tg saturates at a constant value, Tq>g^ = 
— /(O) = (1 — gA)^ /2gA- It is worth recalling that we 
consider the spectrum describing the ensemble-averaged 
IPRs, {Pq), see Sec. III.C.5I If one looks at the scal- 
ing of the typical IPR, Pg^^, the information about rare 
events encoded in the part of the spectrum with /(a) < 
gets lost, and only the range < a < «_ is probed, 
where a± — 2{g]^^ =F 1)^. The corresponding r^^P spec- 
trum reproduces Tq in the range q^ < q < q^, where 



For a nonzero total flux there wiU be additional zero 
modes, in view of th e Atiy ah-Singer index theorem; see, e.g. 
lAh aronov and Cashed lll979l) . 



q± 

Eq, 



= gA^^'^i a-nd becomes linear outside this range, see 
p.42p . This behavior was found for the random 



vector potential problem inlde C. Chamon et al\ 
ICarpentier and LePoussall (|200l[ ): ICastillo et all J l99 
where the spectrum r*^P was studied (in these works the 
notation qc is used for our q^). 



d. Fre e zing. As was found in Ide C. Chamon et al. 
1996 ): ICarpentier and LePoussall (|2001 ): ICastillo et^. 



1997D . with increasing disorder the system undergoes a 
phase transition ("freezing") a,t gA — I- In the strong 
disorder phase, gA > the spectrum takes the form 



riq) = -2{l-g'fqr , 
/(a) = a(8g]l^ - a)/8gA 



-1/2 

q<qc = 9 a 



(6.94) 
(6.95) 



In this phase /(O) = 0, which implies that there is a prob- 
ability of order unity to find a point in the sample where 
the wave function is of order unity. Correspondingly, the 
saturation value of Tq for q > qc is Tq-^q^ = 0. At first 
sight, this may seem to imply that the wave functions are 
localized. The situation is, however, not so simple: a non- 
trivial multifractal spectrum shows a complex structure 
of the wave functions. Further, it can be shown that the 
probability to find a secondary spike in the wave function 
of approximately the same magnitude as the main one 
and separated by a di stance comparable to the system 
size L is of o rder unity ( Carpentier and LeDoussaj l200ll 
lFukuil . l2003h . The nature of these "quasilocalized" wave 
functions is therefore similar to that of critical states in 
ID systems of chiral symmetry, see Sec. IV. El 



e. Density of states. The scaling of the DOS is governed 
by the dynamical exponent z via 



p{E) cx £;(2-^)A 



(6.96) 



In the weak-disorder phase {gA < 1) th e value of z is 
straig htforwardly found to he z = l-l-2q4 (lAltland et~al\ . 
2002.; iLudwig et al\ . \l994 iNersesvan a^.l■ll995^ . How- 
ever, as was already discussed in Sec. lVI.FT2| the freezing 
has also impact on the dynamical exponent; the value of 
z in the strong-disorder phase is given by the second line 
of Eq. (j670)) . 



f. Non-abelian random gauge Held. The problem remains 
exactly solvable in the case of a SU(N) non-abelian gauge 
filed. However, in contrast to the abelian case, where 
one finds a line of fixed points, the theory flows now 
into an isolated fixed p oint, which is a WZ W t heor; 
on the level k = - 2N (ICauxl. Il998l: ICaux et all . Il99; 
iMudrv et all . Il996l: iNersesvan et all \199$- The spec 



trum of multifractality is parabolic with 7 = {N — 1)/N'^ 
and is given by Eqs. (I6.92p . (|6.93p with a replacement of 
2gA by {N-l)/N^. The DOS is given by Eq. 1^:^^ with 
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a dynamical exponent z = 2 — In the case N = 2 

that arises in the two-node model one thus finds z ~ 7/4, 
yielding the DOS scaling p{E) cx E^^'^. 

When both abelian and non-abelian random gauge 
fields are present, they contribute additively to the ex- 
ponents (in the non- frozen regime), as the two s ectors 
of the theory decouple. This yields fMu drv et al. I. I1996D 
the multifractal scaling (|6.92p . (|6.93p with 7 = 2gA 
2gA + {N - 1)/N^ and the DOS scaling ((Oel) with 
z = 1 + 2.g^ + (iV2 - 1)/N\ 



g. Conductivity. A further remarkable feature of the 
models with Co-chiral randomness is the independence 
of conductivity on the disorder strength. We sketch th e 
proof of this statement given in lOstrovskv et al\ (|2006t) . 
The conductivity is given by the Kubo formula 



2T:h 



■Tr 



(6.97) 



where 'Tr' includes both the matrix trace and the spatial 
integration. The chiral symmetry implies the identity 



CT3 



(6.98) 



allowing one to trade all advanced Green functions in Eq. 
()6.97p for retarded ones and thus to present the conduc- 
tivity at zero energy in terms of retarded Green functions. 
Further, in view of the Dirac character of the spectrum, 
the components of the current operator are related via 



'J crs = I] 



y 



(6.99) 



By making use of Eqs. (|6.98p and (|6.99p the Kubo formula 
can be cast in the following form: 



a{E = {]) = -^ V Tr j^'G'^fG 



■a r~iR -a. /-iR 



(6.100) 



a—x,y 



At first glance, this expression may be thought to be zero 
due to gauge invariance. Indeed, the right-hand side of 
Eq. (16. loop is proportional to the second derivative of the 
partition function Z[A] — TrlogG'^[A] (or, equivalently, 
first derivative of the current Tr j"G'^[A]) with respect 
to the constant vector potential A. The gauge invariance 
implies that a constant vector potential does not affect 
gauge-invariant quantities like the partition function or 
the current, so that the derivative is zero. This argu- 
ment is, however, not fully correct, in view of a quantum 
anomaly. The elimination of A amounts technically to a 
shift in the momentum space k — > k — eA, which naively 
does not change the momentum integral. If we consider a 
formal expansion in the disorder strength, this argument 
will indeed hold for all terms involving disorder but not 
for the zero-order contribution. The momentum integral 
/ (fkTT J'^Gq (k) is ultraviolet-divergent and the shift of 
variable is illegitimate. This anomaly was first identi- 
fied in ISchwingeil (|1963 ) for 1 -|- 1-dimensional massless 
Dirac fermions. In the Schwinger model, the polarization 



operator is not affected by an arbitrary external vector 
potential A(a;, t) and is given by the anomalous contribu- 
tion, yielding a photon mass in the 1-1-1 electro dynamics 
( Peskin and Schroedeii Il995l : ISchwinged. I1962D . In the 
present context, the role of A{x,t) is played by the chi- 
ral disorder. The explicit calculation of the zero-order 
(disorder-free) diagram yields (including a factor of two 
accounting for spin) 



=2, ,2 



(27r)2 (t;2fc2 + ^2)2 



4e^ 
Trh 



(6.101) 



Therefore, the dimensionless conductivity acquires the 
universal value 4/7r, with no corrections at any order in 
the disorder s t rength . This result was earlier obt ained in 
iLudwig et all (|l994f) for the abehan case and in ITsvelikl 
(|l995l) for a certain model of non-abelian gauge field. 



4. Disorders preserving Cz chirality: Gade-Wegner criticality 

Let us now turn to the disorder which preserves the 
Gz-chirality, H = —a^T^HasTs; according to Table IVl 
the corresponding coupling constants are /3o, ct±, 7_l and 
(3z- If no time- reversal symmetries are preserved, the 
system belongs to the chiral unitary (AIII) class. The 
combination of Cz-chirality and the time reversal invari- 
ance Tq (only couplings j± and are present) corre- 
sponds the chiral orthogonal symmetry class BDI; this 
model has already been discussed in Sec. lVI.I\2] Finally, 
the combination of C^-chirality and T^-symmetry (cou- 
plings 7_L and (3o) falls into the chiral symplectic symme- 
try class GIL The RG flow and DOS i n thes e models have 
been analyzed in iGuruswamv et al] (|2000l ) in notations 
9A = 7J-/2, gm = /3z/2, Of, = /3o/2, gA' = a±/2. . In 
all the cases, the resulting theory is of the Gade-Wegner 
type, see Sec. IVI.Fl 

In contrast to random gauge field models (Go chi- 
rality), the proof of the universality of the conductiv- 
ity based on gauge-invariance argument does not ap- 
ply to models with Cz chirality. Nevertheless, the zero- 
energy conductivity is found to have the same univcr- 
sal value, a = Ae^/Trh, in the li r nit of weak disorder 
(|Ostrovskv et all . 120061 : iRvu et a/.l . [2007a[ ) . In this case, 
however, there are corrections to this value perturbative 
in the disorder strength. In particular, the leading cor- 
rectio n is found in the second order ( Ostrovskv et all . 
M^,Sa(^) ^{eynh)iPo-P.?. 

The following comment (applicable both to Go and Gz 
chiral models) on the infrared regularization of the prob- 
lem is in order here. This role is played in Eq. (j6.10ip by 
(5, which is an infinitesimal imaginary part in the denomi- 
nator of the Green function. Physically, it has a meaning 
of the inverse electron lifetime or, alternatively, a dephas- 
ing rate, and can be thought of as modelling processes 
of escape of electrons in some r eservoir or some d ephas - 
ing mechanism. As discussed in lOstrovskv et all (|2006h . 
for a different type of the infrared regularizations (i.e. 
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by finite frequency) the critical conductivity will take a 
different value, while still being of order /h. 



5. Dirac Hamiltonians for dirty d-wave superconductors. 

We briefly discuss the application of Dirac fermion 
theory to disordered d-wave spin-singlet superconduc- 
tors with dx2-y2 symmetry. In such systems the su- 
perconducting gap vanishes at four points in the Bril- 
louin zone, and the dispersion relation near these points 
of the Dirac type. One is thus led to analyze the low- 
energy physics of the problem in terms of a four-flavor 
Dirac Hamiltonian. The investigation of symmetries of 
this problem in the presence of different types of impuri- 
tie s and the corresponding RG treatment were pioneered 
in iNersesvan et al\ { 1995( ): for a recent detailed analy- 
sis the reader is referred to lAltland et all (|2002) . The 
main physical quantity of interest for this problem is the 

low-energy be havior of DOS . 

Following lAltland et~ai\ (|2002l ): INersesvan et all 



nodes coupled T Class Criticality 



DOS 



( 19951 ). we specialize on potential disorder, assume that 
the spin-rotation invariance is preserved, but allow for a 
possibility of time-reversal symmetry breaking (i.e. by 
magnetic field in vortices in type-II superconductors). 
In full analogy with the two-flavor model. Sec. IVI.G.li it 
is crucially important, whether the disorder couples the 
nodes or not. Depending on the range of disorder and 
the interval of energies considered, one can distinguish 
three different situations: 

(i) short-range disorder: all four nodes are coupled; 

(ii) long-range disorder: inter-node scattering is negli- 
gible, and the nodes are decoupled. The problem acquires 
then a single-node character; 

(iii) each node is coupled to the opposite one but not to 
the other two. The system then decomposes in two parts 
- each of them describin g a pair of the op p ositel y located 
nodes. It was shown in INersesvan et all (1995) that, in 
view of strong anisotropy of the nodal Dirac Hamilto- 
nians, there is an intermediate energy range where this 
model becomes physically justified. 

Combining these three types of disorder with possi- 
bilities of preserved or broken time-reversal invariance, 
one gets six distinct models. Their symmetries, emerg- 
ing types of criticality, and the corresponding behavior 
of DOS are summarized in Table IVIII When all four 
nodes are coupled, the system is in the conventional lo- 
calized regime of symmetry cl asses CI and C des c ribing 
spin-singlet superconductors ( Senthil and Fish"erl . Il999l : 
ISenthil et aL . ,1998 ) . When only the opposite nodes are 
coupled, the T-invariant problem becomes a model of 
non-abelian random gauge field, and the theory acquires 
the WZ term, leading to the E^^'^ scaling of the DOS, 
see Sec. I VI. G. 31 For broken T-invariance, the two-node 
problem is described by the class-C cr-model with the 
6 = IT topological term, i.e. it is in the SQH transition 
universality class. Sec. IVI.D] The DOS scales as i?^/^ in 
this case as well. Finally, when the nodes are completely 



+ AIII WZW |£;|(i-2sa)/(i+29a) 



+ CI 

+ CI 

- A 

- C 

- C 



WZW 

^-term 
S-term 



non-critical 



TABLE VII Criticality and DOS in d-wave superconductors 
with preserved spin rotation invariance and different types of 
randomness. First column: number of nodes coupled. Sec- 
ond column: presence or absence of time-reversal invariance. 
Third column: symmetry class. Fourth column: type of crit- 
icality. Fifth column: low-energy scaling of DOS. 



decoupled, the Hamiltonian for each of them takes the 
form analyzed in Sec. IVI.G.2I More specifically, when 
T-invariance is preserved, the problem reduces to the 
class- AIII model of random (abelian) vector potential. 
Sec. IVI.G.3[ with the DOS following the non-universal 
power law (jOe]) . p{E) - |£;|(i-29a)/(i+2s^), the case 
of broken T-invariance, the problem belongs to the class 
A and is described by the Pruisken cr-model with topo- 
logical term, i.e. it is in the universality class of the IQH 
critical point. In the latter situation, the DOS is uncrit- 
ical. 

We also mention two further types of disorder that 
have been studied in detail in the literature. In the case 
of magnetic impurities botn T-invariance and spin rota- 
tion sym metry are broken and thus the system belongs t o 
class D (Bocauet et all . I2OOOI : ISenthil and Fisheii l2000f ). 
see Sec. I VI. El For the case of perfect nesting, where the 
chemical potential /i=0, and very strong potential scat- 
terers one encou nters the symm etry class AIII leading to 
the Gade phase (|Altlandl . l2003 ). 



VII. CONCLUDING REMARKS 

This concludes our review of physics of Anderson tran- 
sitions between localized and metallic phases in disor- 
dered electronic systems and of associated critical phe- 
nomena. In our opinion, during the recent years a quite 
detailed understanding of these transtions has emerged, 
including such salient features as symmetry and topol- 
ogy classification, mechanisms of criticality in quasi-lD 
and 2D systems, and wavefunction multifractality. We 
have attempted to give a systematic exposition of these 
issues and hope that this will help interested researchers 
to navigate in this remarkably rich field. 

For several reasons (including limits in space, time, 
and - last but not least - our expertise), we were not 
able to discuss all aspects of the problem. The selection 
of material in the review certainly reflects the authors' 
personal perspective of the field. We apologize to all 
researchers whose work is not sufficiently discussed in the 
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review. Probably, the most important issue out of those 
that we largely left apart in the review is the effects of 
electron-electron interactions. (We only briefly duscussed 
it in Sec. IVI.CI in the context of the IQH transition.) 
This by itself is a very rich and complicated problem; we 
restrict ourselves to just a few remarks and references. 



A. A few words about interaction effects 

Physically, the impact of interaction effects onto low- 
temperature transport and localization in disordered 
electronic systems can be subdivided into two distinct 
effects: (i) renormalization and (ii) dephasing. 



a. Renormalization. The renormalization effects, which 
are governed by virtual processes, become increasingly 
more pronounced with lowering temperature. The im- 
portance of such effects in diffusive low-dimensional sys- 
tems was d emonstrated in by Altshuler and Aronov, see 
the review (|Altshuler and AronovL 11984 ). To resum the 
arising singular contributions, Finkelstein developed the 
RG approach based o n the g-model for an interacting sys- 
tem, see iFinkelstei"3 (|l990f ) for a review. This made pos- 
sible an analysis of the critical behavior at the localiza- 
tion transition in 2 -|- e dimensions in the situations when 
spin-rotation invariance is broken (by spin-orbit scatter- 
ing, magnetic field, or magnetic impurities). However, 
in the case of preserved spin-rotation symmetry it was 
found that the strength of the interaction in spin-triplet 
channel scales to infinity at certain RG scale. This was 
interpreted as some kind of magnetic instability of the 
system; for a detailed exposition of propose d scenarios 
see the review (|Belitz and KirkpatricR Il994l ). 

Recently, the problem has attracted a great deal of at- 
tention in connection with experiments on high-mobility 
low-density 2D electron structures (Si MOSFETs) giv- 
i ng an evidence in favo r of a metal-insulator transition 
( Abrahams et a/.l . l200lh . Whether these results are due 
to a true metallic phase existing in these systems or, else, 
are explained by interaction effects at intermediate ( "bal- 
listic") ter nperatures remains a debated iss ue. In a re- 
cent work ( Punnoose and Finkel'steinl . l2005l ) the RG for 
CT-model for interacting 2D electrons with a number of 
valleys > 1 was analyzed on the two-loop level. It was 
shown that in the limit of large number of valleys N (in 
practice, = 2 as in Si is already sufficient) the tem- 
perature of magnetic instability is suppressed down to 
unrealistically low temperatures, and a metal-insulator 
transition emerges. 

The interaction-induced renormalization effects be- 
come extremely strong for correlated ID systems (Lut- 
tinger liquids). While ID systems provide a paradigmatic 
example of strong Anderson localization, a sufficiently 
strong attractive interaction can lead to delocalization 
in such systems. An RG treatment of the corresponding 
localization transition in a disordered interacting ID sys- 



tems was deve loped in Giamarchi and Schuld ( 19881 ). see 
also the book ( Giamarchil 2004 ). 



b. Depliasing. We turn now to effects of dephasing 
governed by inelastic processes of electron-electron 
scattering at finite temperature T. The dephasing 
has been studied in great detail for metallic systems 
where it provides a cutoff fo r weak localization effects 
( Altshuler and Aronovl [r984j ) . As to the Anderson tran- 
sitions, they are quantum (zero-T) phase transitions, and 
dephasing contributes to their smearing at finite T. This 
has been discussed in Sec. IVI.C.lOl in the context of dy- 
namical scaling at the IQH transition. There is, how- 
ever, an interesting situation when dephasing processes 
can create a localization transition. We mean the sys- 
tems where all states are localized in the absence of in- 
teraction, such as wires or 2D systems. At high tem- 
peratures, when the dephasing is strong, so that the 
dephasing rate t^^{T) is larger than mean level spac- 
ing in the localization volume, the system is a good 
metal and its conductivity is given by the quasiclassical 
Drude conducti vity with relatively small we ak localiza- 
tion correction ( Altshuler and Aronovl [19841 ). With low- 
ering temperature the dephasing gets progressively less 
efficient, the localization effects proliferate, and eventu- 
ally the system becomes an Anderson insulator. What 
is the nature of this state? A natural question is 
whether the interaction of an electron with other elec- 
trons will be sufficient to provide a kind of thermal bath 
that would assist the v ariable-range hopping transport 
( Fleishman et aLl . ll978t ). as it happens in the presence of 
a phonon bath. The answer to thi s question was given 
bv [Fleishman and AndersonI ( 1980f ). and it is negative. 
Fleishman and Anderson found that at low T the inter- 
action of a "short-range class" (which includes a finite- 
tange interaction in any dimensionality d and Coulomb 
interaction in d < 3) is not sufficient to delocalize other- 
wise localized electrons, so that the conductivity remains 
strictly zero. In combination with the Drude conductiv- 
ity at high-T this implies the existence of transition at 
some temperature Tc- 

This conclusion was recently corroborate d by an analy- 
sis ()Basko et aLl . l2006HGori^ et al. in the frame- 
work of the idea of An derson localization in Fock space 
( Altshuler et all . ll99l\) . In these works the temperature 
dependence of conductivity a{T) in systems with local- 
ized states and weak electron-electron interaction was 
studied. It was found that with decreasing T the sys- 
tem first shows a crossover from the weak-localization 
regime into that of "power-law hopping" over localized 
states (where cr is a power-law function of T), and 
then undergoes a localization transition. The transition 
is obtain ed both within a s elf-consistent Born approx- 
imation ( Basko et all . |2006| ) and an appr oximate map- 
ping onto a model on the Bethe lattice ( Gornvi et all 
|2005[ ). The latter yields also a critical behavior of a{T) 
above Tn, which has a characteristic for the Bethe lat- 
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tice non-power-law form lna{T) ~ (T — Tc)~^/^, see 
Sec. III. C. 91 Up to now, this transition has not been ob- 
served in experiments^'^, which indicate instead a smooth 
crossover fro m the metallic to the insulating phase with 
lowering T (Hsu and Vallesl Il995t iKhavin eToII . Il998l: 



iMinkov all 120071 : IVan Keuls et al\ . Il997t ). The rea- 
son for this discrepancy remains unclear. A recent at- 
tempt to detect the transition in numerical simulations 
also did not give a clear confirmation of the theory 
( Qganesyan and Husd . l2007[ l. possibly because of strong 
restrictions on the size of an interacting system that can 
be numerically diagonalized. 



B. Experimental studies of localization transitions 

Of all the localization transitions, the best studied ex- 
perimentally is the IQH transition. We have discussed 
the corresponding experimental findings in Sec. IVI.C.lOl 
Its superconducting counterparts (SQH and TQH transi- 
tions) have not been observed yet, although several phys- 
ical realizations of them were proposed. Sec. IVI.D.ll and 
IVI.E.ll Much effort has been invested in research on 
strongly interacting disordered 2D systems in zero mag- 
netic field but it remains controversial whether what is 
observed there is a true metal-insulator transition, see 
Sec. IVII.AI 

Below we briefly review the situation with the experi- 
mental observation of the Anderson transition in 3D elec- 
tronic and optical systems. For electronic systems (doped 
semiconductors. Sec. IVII.B.ll) the localization transition 
has been observed unambiguously. However, a theoreti- 
cal analysis of the critical behavior is complicated by the 
presence of Coulomb interaction, which modifies the crit- 
ical behavior. As a result, one can not expect that the 
experimentally extracted critical exponents agree with 
numerical values obtained from computer simulations 
on non-interacting systems. The localization of light 
fSec. lVH.B.2[) has an advantage in this respect, since the 
photon-photon interaction is negligibly small. However, 
it turns out that implementation of sufficiently strong dis- 
order so as to reach the Anderson transition and strong 
localization of light in 3D is a remarkably complicated 
endeavor. The second major difficulty is posed by the 
absorption. 



1. Anderson transition in doped semiconductors 



Si:As, Ge:Sb. In most of the works, samples with a sub- 
stantial degree of compensation [i.e. acceptors in addi- 
tion to donors, e.g. Si:(P,B)] were used, which allows one 
to vary the amount of disorder and the electron concen- 
tration independently. On these samples, values of the 
conduct ivity exponent s in the vicin i ty of s » 1 were re- 
ported (Field and Rosenbaum', '1985'; Hirsch et ai, 'l988|; 
Thomas et ai, .198 2; Zabrodskii and Zinov'eva, 198^ 
with scattering of values and the uncertainties of the or- 
der of 10%. A similar result was obtained for an amor- 
phous material Nb^^Sii-^:. (Recall that in 3D s is ex- 
pected to be equal in 3D to the localization length expo- 
nent v according to the scaling relation s = ^{d — 2)) 

On the other hand, the early study of the transition in 
undq p ed Si:P (iPaalanen et all . Il982t iRosenbaum et all 
119831: iThomas et all Il983l ) gave an essentially differ- 
ent result, s ~ 0.5. A reso l ution of th is discrepancy 
was proposed in IStupp et all (119931 [ 1994 ) where it was 
found that the actual critical region in an uncompen- 
sated Si:P is rather narrow and that the scaling anal- 
ysis restricted to this range yields s » 1.3. A more 
recent study ( Waffenschmidt et all Il999l ). which em- 
ployed uniaxial stress to tune through the transi t ion (a s 
used in Paalanen et al. I (|1982| ): IRosenbaum et al. 1 (I1983D : 
Irh omas et all (I1983D T has essentially confirmed these 
conclusions, yielding s — l.OibO.l, in agreement with the 
values obtained for samples with compensation. Further, 
in this work dynamical scaling near the transition with 
varying temperature was demonstrated; the correspond- 
ing dynamical exponent was found to be z = 2.94 ± 0.3. 
Good scaling was also o bserved in a similar expe riment 
on uncompensated Si:B ( Bogdanovich et aZl . ll999f) . how- 
ever with somewhat different critical exponents (s « 1.6, 
z « 2). A possible explanation for this discrepancy is 
that the temperatures reached in this work were not suf- 
ficiently low. Another possibility is that Si:B belongs to a 
different universality class, in view of stronger spin-orbit 
scattering. 

The fact that the experimental value s w 1 found in 
the majority of works differs from what one would expect 
based on numerical studies for non-interacting systems 
is not surprising, since the Coulomb interaction affects 
the critical exponents. For a more detailed discussion of 
the experimental data and their comparison with theo- 
retical expectations for the Anderson transitions in the 
presence of Coulomb interact ion the reader is refered to 
iBehtz and Kirkpatrickl (|l994 ). 



The 3D localization transition was extensively stud- 
ied on doped semiconductor systems, such as Si:P, Si:B, 



Of course, in a real system, phonons are always present and 
provide a bath necessary to support the hopping conductivity at 
low T, so that there is no true transition. However, when the 
coupling to phonons is weak, this hopping conductivity will have 
a small prcfactor, yielding a "quasi-transition" . 



2. Anderson localization of light 

Experimental efforts to observe the Anderson localiza- 
tion with light turned out to be very challenging. The 
main difhculty is, that the characteristic signature of 
wave localization, the exponential decrease e~^l^ of the 
transmission with the system size L, is often hard to 
disentangle from another exponential e~^"/^ originating 
from photon absorption. 
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iGarcia and GenacS ( 199ll ): [Genack and Garcial ( 199ll ) 
studied the transmission of microwave radiation through 
a random mixture of aluminum and Teflon spheres inside 
a copper tube. In the vicinity of the localization transi- 
tion, they report power law scaling of the effective diffu- 
sion constant D{L) ^ L^^. The presence of strong ab- 
sorption in the microwave system, makes it complicated 
to interpret the data that were reported as evidence for 
the localization of light. 

To diminish photon losses. IWiersma et all (|l997l ) used 
a powder from /im gallium arsenide crystals employing 
near infrared radiation. However, still the exponential 
signal of localization reported in this work is not undis- 
puted, again due to th e presence of residual absorption 
( Scheffold et all . Il999l ) . A similar study in silicon pow- 
ders did not show any hint of localization ( Rivas et all . 
I1999D . 

In order to overcome the notorious problem of separat- 
ing localization fro m absorption, an intere sting proposal 
has been made bv lChabanov et 'al. 1 (|2000D . Localization 
is accompanied by very large mesoscopic fluctuations in 
the spectral function which carry over also to the fre- 
quency dependent transmission function Tiuj). Hence, 
the measurement of the statistical properties of T[uj), 
e.g. its variance, could provide a sensitive means to un- 
cover quantum interference effects. Using this method 
the authors were able to confirm localization in quasi- 
one-dimensional waveguides; an application to a three- 
dimensional system, a mixture of aluminum spheres, did 
not show signatures of strong localization. 

To summarize, up to now there does not seem to be 
a convincing demonstration for the Anderson localiza- 
tion and t ransition in 3D syst ems. However, in very re- 
cent work lStorzer et all ( 20061 ) report the observation of 
an anomalous time dependency of the light diffusion in 
a Ti02 powder indicating the vicinity of the Anderson 
critical point. Another very promising research direc- 
tion has app e ared with the advent of photonic crystals 
(|Busch et aLl . l2007t) . where in the presence of disorder a 
localization transitio n should take pla.ce for states near a 
photo nic band edge ( Busch and Johnl Il999t iJohnl . 1 1984 
I1987D . 
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